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ABSTRACT 


The  two-dimensional  ablative  heat,  transfer  computer 
program  generated  under  contract  AK>4(6ll;-9?14  was 
refined  and  extended  to  handle  anisotropic  materials 
more  than  one  charring  material  and  reradiation  at 
the  heated  surface.  Results  cf  the  two-dimensional 
program  were  compared  to  experimental  data  to  deter¬ 
mine  effective  values  of  material  properties  used  in 
the  analysis  to  simulate  gas  generation  and  cracking 
r  ’actions,.  The  resulting  effective  properties  for 
two  ablative  materials  were  used  in  a  parametric 
study  generating  basic  information  for  the  design 
of  ablative  systems  in  liquid  rocket  thrust  chambers 
Graphs  anti  charts  showing  the  variation  of  thermal 
penetration,  char  depth  and  surface  erosion  are 
included. 


CONTENTS 


Review  of  2D-AB1ATE  Program  . .  3 

General  . . 3 

Mathematical  Model  .........  5 

Numerical  Solution  .......  a  . . 9 

Results  . . 13 

Program  Structure . „ . 13 

Task  1,  Extension  of  the  Two-Dimensional  Ablation  Program  ....  17 

General  Modifications . 18 

Task  l.A,  Radiative  Exchange  at  the  Exposed  Inside  Surface  ...  23 

Task  l.B,  Thermal  Analysis  of  Anisotropic  Materials  .  31 

Task  l.C,  More  Than  One  Charring  Material . 48 

Task  2,  Effective  Ablation  Properties  . . 51 

Introduction  .  ...........  51 

Properties  . . 53 

Estimation  of  Effective  Thermal  Conductivity  .......  64 

Task  3j  Designer's  Guide . 81 

Significant  Parameters . 81 

Nozzle  Throat  Erosion . 103 

Two-Dimensional  Ablation  .........  108 

Summary  of  Design  Procedures  . . 131 

Recommendations  for  Ihture  Efforts  .  135 

Nomenclature  . . 137 

References . 141 

Distribution . 169 

Appendix  A 

Calculation  of  View  Factors  by  the  Disk  Method  .......  145 

Appendix  B 

Development  of  the  Anisotropic  Energy  Equation  .......  151 


v 


C  X 

Finite  Difference  Analogs  of 


sendix 


Simplified  Calculation  of  Nozzle  Heat  Transfer  Coefficients 


and  Radiative  View  Factors 


Simplified  Calculation  of  Char  Hate  in  the  Overwrap  of  an 


Ablative  Wall 


ILLUSTRATIONS 


1.  Typical  Thrust  Ch saber  Configuration  and 

Boundary  Conditions  ••••»....  *  ...  4 

2.  Typical  Mesh  Configuration  and  Categorization 

of  Be3ulting  Mesa  Points.  .  ^ . .  10 

3.  Thrust  Chamber  Wall  Configuration  Used  For 

EAFB  Firing  . .  14 

4.  Structure  of  the  Belinked  Program .  19 

5.  Point  Arrangement  for  Discretizing  the 

Continuity  Equation  . .  21 

6.  Subareas,  B ^ 10,  Used  to  Calculate  Input  Viev 

Factors  G.  ,  for  the  Checkout  of  Task  l.A .  28 


7-  Comparison  of  Surface  Temperature  Histories  Computed  in 
Checkout  of  Task  l.A  for  the  Chamber  Zone  and 

Throat  Region .  29 

8.  Comparison  of  Surface  Temperature  Histories  Computed 
in  Checkout  of  Task  l.A  for  the  Throat  Region 

and  Exit  Nozzle  .  . .  JO 

9.  Various  Hesh-Boundary  Configurations  Affecting 

2 

Formulation  of  Difference  .Analog  6  T  .......  Jk 

xy  o 


10.  Anisotropic  Conductivity  for  Checkout .  40 

11.  Anisotropic  Checkout  Cases  Comparison '  of  Temperature 

After  12  Seconds  of  Firing  .  42 


12.  Checkout  Cases  for  Charring  in  More  Than  One  Material  ...  49 

13-  Density-Specific  Heat  Product  for  Ablative 


Chamber  Materials  .  Jk 

14.  fraction  Resin  Pyrolyzed . .  Jk 


rii 


15-  Variation  of  Pyrolysis-Gae  Bathaipy 

With  Temper store . . .  58 

16.  dermal  Conductivity  of  Phenolic  Refrasil  .......  58 

1?.  OFg/Amine  Test  Motor  (Typical),  Model  Used 

in  2B-Afclate  . .  65 

18.  Typical  Variation  of  the  Heat  Transfer  Coefficient 
With  Longitudinal  and  Circumferential  Location  for 
Self-Impinging  Doublet  Injector  Employed  in  Ablative 

Hotor  Firings . 67 

19.  Axial  Distribution  of  the  Adiabatic  Wsll  Temperature 
and  the  Gas-Side  Heat  Transfer  Coefficient  as  Derived 

From  an  ATJ-Graphite  Reference  Test  Cnee .  72 

20.  Computed  and  Experimental  Temperature  Histories  For 

an  ATJ-Graphite  Case  (Test  l)  .  73 

21.  Temperature  Profiles  in  Carbon  Phenolic  Ablators  at 
280  Seconds  for  an  Adiabatic  Wall  Temperature  of 
2800  F  and  a  Heat  Transfer  Coefficient  of 

O.OQ&3  Btu/in^-sec-F  . . 74 

22.  Apollo  Qualification  Test  Model  For  2D-ABLATE .  76 

23.  Throat-Station  Backwall  Temperature  .  78 

24.  Nozzle-Station  Backwall  Temperature  ........  79/80 

25.  Effect  of  Thermal  Conductivity  on  Char  Front  Advance 

For  an  Adiabatic  Wall  Temperature  of  3500  F .  86 

* 

26.  Effect  of  Thermal  Conductivity  on  Char  Front  Advance 

For  an  Adiabatic  Wall  Temperature  of  5000  F .  86 

27.  Effect  of  Thermal  Conductivi+y  on  Char  Front  Advance 

for  an  Adiabatic  Wall  Tempentu.  e  of  7000  F .  87 

28.  Effective  of  Thermal  Conductivity  at  High  and  Very 

Law  Values  of  Heat  Transfer  Coefficient .  87 


viii 


29.  Effect  of  Thermal  Cond  xtivity  on  Thermal 

Penetration  Through  Ablative  Walls . ,89 

30.  Effect  of  Char  Thermal  Conductivity  on  the  Gaa-Side 

Surface  Temperature  of  Ablative  Walls  ........  89 

31.  Effect  of  Adiabatic  Wall  Temperature  on  Char  Front 

Advance . 91 

32.  Enthalpy  Parameter  ij;  as  a  Function  of  Adiabatic 

Wall  Temperature . 91 

33*  Effect  of  Adiabatic  Wall  Temperature  on  Thermal 

Penetration  Through  Ablative  Walls  .  93 

34.  Effect  of  Adiabatic  Wall  Temperature  on  the  Surface 

Temperature  of  Ablative  Walls  .  93 

35«  Eelatioa  Between  Adiabatic  Wall  Temperature  and  Surface 
Temperatures  Expressed  as  Thu..  ^omplished 

Temperature  Difference  .  94 

36.  Effect  of  Resin  Fraction  on  Char  Front  Advance  For  High 

Thermal  Conductivity  and  High  Heat  Transfer  Coefficient  ,  .  96 

37-  Effect  of  Resin  Fraction  on  Char  Front  Advance  For  Low 

Thermal  Conductivity  and  Low  Heat  Transfer  Coefficient  .  .  96 

38.  Effect  of  Resin  Fraction  on  Thermal  Penetration  ....  97 

39«  Effect  of  Resin  Fraction  on  Surface  Temperature  of 

Ablative  Walls  . 97 

40.  Effect  of  Heat  Transfer  Coefficient  on  Char  Front 

Advance  Through  a  High  Conductivity  Ablative  .  99 

41.  Effect  of  Heat  Transfer  Coefficient  on  Char  Front 

Advance  Through  Low  Conductivity  Ablative  .  99 

42.  Effect  of  Heat  Transfer  Coefficient  on  Surface 

Temperature  of  Low-Conductivity  Ablative  .  1Q0 


lx 


45-  Effect  of  Seat  Transfer  Coefficient  e=a  Thermal 
Penetraiien  Through  Low-Cccdactivify 
Ablative  Walls  .  -  -  ,  . . . 

44  r  Effect  of  Charier  Badius  on  Char  Front  Advance  . 

45-  Recession  of  AxJ  Graphite  Throat  Inserts  in  a 

Water  Vapor  Environment  at  500-psi  Charter 
Pressure  ..  . 

46.  Secession  of  ArJ  Graphite  Throat  Inserts  in  a 

Hydrogen  Ehrironseni  at  500  psia  . 

47-  Recession  of  Silicon  Carbide  Throat  Inserts  . 

48.  Thrust  Chamber  Configuration  for  Computer  Ban  1  - 

49-  Thrust  Chrdier  Configuration  for  Computer  Bans  2  and  5- 

50.  Thrust  Chamber  Configuration  for  Computer  Buns  4  and  7- 

51.  Thrust  Chamber  Configuration  for  Computer  Buna  5  and  6, 

52.  Char  Fronts ,  Gas-Side  Temperatures,  and  Outer  Skin 

Temperatures  for  Phenolic  Carton  Cloth  Ablatives 
After  500-Second  Firing  With  the  Propellant 
Combination  of  CTF/NgH^  .  . 

53-  Char  Fronts,  Gas-Side  Temperatures,  an?  Outer  Skin 
Temperatures  for  Phenol ic-Eefrasil  Ablatives  After 
300-Seeond  Firing  With  the  Propellant  Combination 
of  N204/50<£  -  5#  TOffl . 

54.  Comparison  of  Combustion-Zone  Char  Advance  By 
Two-  and  One-Dimensional  Computer  Programs 
—  High  Chamber  Pressure —  . 

55-  Comparison  of  Char  Advance  in  Nozzle  at  an  Area  Ratio 
of  2-9  by  Two-  and  One-Dimensional  Computer  Programs 
—  High  Chamber  Pressure  —  . 


100 


105 

106 
107 
111 
112 
11> 
114 


117 


118 

119 

119 


x 


Corasriscn  of  Chmr  idr-nce  fhresgh  Phenolic  Hefrssil 
Ceshssti. Chamber  ¥felis  by  Tvo-  rad  One-DiaensicnaI 
Cecsoier  Progress  at  7-3  Inches  Pros  Injector  ....  121 

Comparison  of  Char  Advance  Tcrtugh  Phenolic  Ha fr as  11 
Ccrbis tiers  Chamber  Vails  Hr  Tvo-  end  On  e-Hi  szensiosal 
Cossuter  Progress  at  8-5  Inches  Pros  Injector  ....  121 

Comparison  of  Char  Advance  in  Combustion  Chamber 
By  fvo-  and  One-Pi  sessional  Co  ranter  Pro/jxams 

at  8-5  Inches  Pros  Injector  . .  122 

Comparison  of  Char  Advance  in  Combnstiim  Chamber 
by  Two—  and  One-Dimensional  Coranter  Progress 

at  7-0  Inches  Pres  Injector  .  199 

Comparison  of  Char  Advance  in  Combustion  Chamber 
by  Ivo-  and  One— Dimensional  Computer  Prograss  at 
5-6  Inches  Pron  Injector  —  Ben  2,  I cv  Chamber 

Pressure  —  125 

Comparison  of  Char  Advenes  in  Combustion  Chamber 

by  Ivo—  and  One-Dinensional  Technique  at  5-6  Inches 

proa  Injector  —  Bun  3 ,  Lav  Chamber  Pressure  —  ...  123 

Comparison  of  Char  Advance  in  a  Phenolic-Befrasil  Soft 

Throat  by  Two-  and  One-Dimenaional  Computer  Programs  .  .  125 

Comparison  of  Char  Advance  in  a  Phenolic-Befrasil 

Soft  Throat  by  Two-  and  One-Dimensional  Computer 

Programs  . . 125 

Approach  of  Surface  Temperature  to  Adiabatic  Vail 

Temperature  for  Relatively  High  Chamber  Pressure  .  .  .  126 

Approach  of  Surface  Temperature  to  Adiabatic  Wall 

Temperature  for  Belatively  Lav  Chamber  Pressure  .  .  .  126 

Design  Procedure  for  Ablative  Eoeket  Engines  ....  134 


ri/iii 


-5? 


Ssba r«as  3.  2=3  Ticv  Eect«rs  E.  .  Used  a* 
3 

Inset  For  Task  i  .  J  Sfr^rksnt  .... 


Recogga»adea  Dens  ity— Specif  i  c  Beat  Profcrfc  Icsot 
ho  2B-Jblzie  Program  For  Ablative  ¥&I1  and  Hard 
Throat  Xateriaxs  . . 


Effective  Ssthal^y  of  Phenolic  Resin  Pyrolysis  Tspors 

Recesssended  Tbaraal  Conductivity  Input  for  2B-A2LS2S 

Erosion  Input  Data  for  Carbon-Cloth  He isf creed 
Phenolic  Ablators  end  A3SJ  Graphite . 

Erosion  Inpai  Data  for  Phenol ic-Sefrasil 
Mstire  Vails  ............ 


Erosion  Input  Bats  lor  Silicon  Carbon  . 

Son  Conditions  for  G?,-/Firirgs  (3ef .  22)  Selected  for 

Evaluation  of  Effective  Properties  for 

rh^fiuxic/ CcTinje  Cloth  .......... 

Location  of  Thermocouples  Used  to  Heasure  Experimental 
Temperature  Eistories  Used  for  Evaluation  of 
Effective  Properties  for  Phenolic/Carben  Cloth  . 

fheroacouple  Bata,  GFg/isine  Tests  . 

Eon  Conditions  lor  2D-JLBLATE  Paranetrie  Cospatationg  . 

Adiabatic  Vail  Temperatures  end  Throat  Heat  Transfer 
Coefficients  for  Tvo-Dinens  icrnal  Cocpater  Buns  . 

Char  Penotraticn  Through  Conhustion  Chancer  Wall 
of  Apollo  Kozzle  . 


xiii/riv 


-dr 


27 


yy 

57 

60 

62 


65 

65 


54 


6S 

69 

109 


115 


130 


XSJEODEiCilGS 


Effective  design  of  rocket  e^ime  thrrysi  chambers,  which  are  passively 
cooled  by  beat  ®imk  cr  ablative  tes’aimes ,  requires  predictions  of  the 
effect  of  high- v^oeraior*  e&ninsiion  products  os  candidate  tfcrnst  ebrabgr 
•wall  materials.  Gdsr  Contract  AJ&4{6ll)-9714,  Bockatdyne  conducted  a 
12-manth  program  estixled  "Effect  of  Socket  Engine  Combustion  cm  Chamber 
Materials,”  la  which  cmsaier  pntgrsms  vere  developed  far  the  numerical 
sc  let  ion  of  one-  and  oro-dimems  laral  models  of  transient  ablation  a  rd  heat 
conduction.  3e tailed  reports  of  this  program  can  be  found  in  Eei.  i  and 
2  .  including  discussions  of  the  physical  and  cathsnatical  node  Is  employed, 

the  numerical  procedures  developed  for  tms  solution,  and  a  comparison  of 
cosnseied  results  to  test  data.  Opera  ting  instructions  for  the  corstrter 
codes  are  given  in  Bef.  3  and  4. 

In  April  196b,  a  12-month  sequel  program  vss  initiated  at  Ro  ckeidyne  under 
Contract  A5U4(6ll)-il415,  entitled  " designer's  Guide  and  Computer  Program 
for  Ablative  Materials  in  liquid  Socket  Thrust  Chambers.”  The  objectives 
of  the  program  vers  to  extend  the  scope  of  the  two-dimensional  program 
/‘ir^iaijTsj  and  casdsci  a  ccssrehessiTe  parametric  stud”.  Tbs  results  of 
this  study  vere  to  be  presented  in  a  form  convenient  for  design  of  ablative 
rocket  engine  thrust  chambers .  The  vork  vas  performed  in  four  phases  Gr 
tasks : 


1.  The  tvo-dinensional  program  vas  extended  to  handle  anisotropic 
materials,  multiple  charring  ablators,  and  the  effects  of 
reradiaiion  at  the  heated  surface. 

2.  Cemputi  d  results  vere  compared  to  existing  experimental  data  to 
determine  values  of  effective  thermal  material  properties  used 

in  the  2D-ABL4TB  program  to  simulate  the  chemical  reactions  vhich 
occur  in  the  interior  of  the  charring  materials. 

3.  A  parametric  study  vas  performed  to  determine  thermal  and  ablative 
effects  as  functions  of  parameters  describing  thrust  chamber  opera¬ 
ting  conditions,  geometry,  and  materials  in  liquid  rocket  engines. 
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4.  G£?ts  asae  care?  vsrs  *3Hi3t«a  cn  tie  lasis  of  tbs  resnits  of 
3  ahs=e  for  sis  in  tbs  design  ct  ihrssi  cbsnibers.  Is  a  separate 
report,  an  operators  rareas!  -as  vritien  for  the  firs  I  version  g£ 
tbs  2WSii3  pic^ras,  ircisaisg  ibe  capabilities  added  in  1  store. 


Is  tbs  sabseqntat  sections  ei  ibis  report,  a  review  of  tea  origins!  2B- 
ABIAI3  asagras  is  first  gives.  lids  ia  foi loved  ay  separate  preemixi i esr 
azss  -Lseassicss  of  tee  v&riGsr  rbss&s  of  verk  accacnlished. 
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SSBsi&i: 


!!«:««>  it  i*  assssad  ice  rteier  is  farilifi?  =c:t»s  ih»  asierial  pre¬ 

sented  in  5 si.  2,  £  brief  review  of  the  2D~AELi52  cserpoier  pragxsa  is 
gives  here  for  ccnvenierce . 

f»t  £D45aB  progress  vas  developed  to  help  fill  -ate  void  caused  by  the 
scarcity  of  c  sspreh  ess  ire  ccesxrier  -codes  for  analyzing  tvo-feeanaosl 
transient  cblatissi  end  best  ctardcsiien  probless  is  roebsi  engine  i crosi 
chasber  vails,  fee  analysis  vss  performed  using  cylindrical  coordinates, 
axial  sad  radial.  iGr  the  will  safari ais  (as  sany  63  five)  of  an  aodsya- 
setrie  thrust  charter  cf  general  configuration  (rig.  l).  Predicted  by 
the  pregren  are  tespe rainre  and  pyrolysis  gas  s asa  flux  distributions  in- 
decsd  throughout  the  vail  rsaieriale  and  recession  rates  at  the  surface 
exposed  to  tee  hot  ccsbnstiaa  gases.  Following  are  scene  of  ibe  sore 
isrxrioirt  pregrs3  features: 

1.  Ice  thrust  chaster  geometry  simulated  using  as  cany  as  40 
quadratic  segments,  inns,  materials  with  arbitrarily  curved 
boundaries  can  be  bandied. 

2.  Curvature-preserving  techniques  are  used  to  obtain  a  second  order 
accurate  approximation  of  the  normally  directed  beat  fluxes 
encountered  at  the  curved  rate  rial  boundaries. 

3.  In-depth  charring  vitbin  the  wall  materials  and  erosion  of  the 
exposed  inside  vail  surface  are  treated  in  the  program,  and 
transpiration  effects  caused  by  gases  generated  vithin  the  vail 
materials  and  at  the  exposed  surfaces  are  accounted  for. 

4.  Physical  and  chemical  properties  say  be  specified  as  functions 
of  temperature  for  each  thrust  edsEber  material. 

5.  The  numerical  procedures  employed  are  stable*  under  most  circun- 
stances  encountered  in  thrust  chamber  applications,  thus 


*Stabie  in  the  sense  that  errors  introduced  into  the  analysis  do  not  grov 
during  the  course  of  the  calculation 


Figure  1.  Typical  Thruat  Chamber  Configuration  and  Boundary  Couditiono 


permitting  the  use  of  relatively  large  time  steps.  In  particular, 
a  generalization  of  the  unconditionally  stable  Peaceman-Eachford 
alternating  direction  method  (Eef.  5)  was  employed  to  discretize 
the  energy  equation  rather  than  the  more  frequently  used,  condi¬ 
tionally  atable,  explicit  forward  difference  method  which  often 
requires  the  use  of  prohibitively  small  time  increments  to  avoid 
error  growth. 

6.  Complex  duty  cycles  of  intermittent  engine  firing  can  be  simulated, 
including  either  steady  bursts  of  firing  followed  by  soakout  or 
high-frequeney  pulsing. 

MATEEM4TI  CA.L  MODEL 

In  formulating  the  mathematical  model  and  developing  the  numerical  solution, 
the  following  assumptions  were  made  concerning  charring  materials  treated  in 
the  analysis: 

1.  In  any  thrust  chamber  configuration,  there  is  no  more  than  one 
charring  material  (this  limitation  was  removed  during  the  subse¬ 
quent  effort  presently  being  reported).  It  is  a  continuous 
material  with  continuously  changing  material  properties. 

2.  Chemical  reactions  in  the  charring  region  are  treated  in  depth; 
i.e. ,  they  are  permitted  to  take  place  throughout  a  continuous 
range  of  temperature  rather  than  being  confined  to  a  single 
"interface"  of  constant  temperature.  (By  specifying  a  very  narrow 
temperature  range  for  the  reaction,  an  interface  model  can  be 
simulated  when  required,  e.g. ,  in  the  case  of  vaporization  of  a 
metal  filling  a  solid  porous  matrix.)  These  reactions,  however, 
are  simulated  thermodynamically  (rather  than  kinetically)  in  the 
energy  and  continuity  equations.  They  include  a  maximum  of  three 
gas  generation  reactions  (pyrolysis,  a  char-reinforcement  re¬ 
action,  and  a  further  decomposition  of  the  solid  product  of  this 
reaction)  plus  cracking  of  the  generated  gases. 


3.  The  porous  char  is  cooled  ccnvectively  by  the  gases  generated  in 
the  charring  material  and  the  gas  and  char  temperatures  at  any 
point  are  identical.  Conduction  of  heat  within  the  gas  is  assumed 
to  be  negligible  in  comparison  with  that  in  the  char. 

4.  Gi<s  density  is  assumed  to  be  negligible  in  comparison  to  the  char 
density. 

5.  The  generated  gas  mass  flux  in  the  porous  char  region  is  assumed 
to  be  oriented  in  the  direction  of  the  temperature  gradient  vector. 

In  addition,  the  following  ground  rules  were  followed  in  analyzing  tie 
effects  of  surface  erosion  and  recession: 

1.  Any  of  the  thrust  chamber  wall  materials  exposed  to  the  hot  com¬ 
bustion  gases  are  subject  to  surface  erosion.  For  each  exposed 
material,  the  surface  lemoval  modes  include  melting,  vaporization, 
and  as  many  as  three  independent  chemical  reactions  with  components 
of  the  combustion  gas. 

2.  For  vaporization  and  chemical  reactions,  the  removal  rcte  is  com¬ 
puted  based  on  mass  diffusion  in  the  boundary  layer  coupled  with 
chemical  kinetics  at  the  wall  surface. 

The  energy  and  continuity  equations  used  in  the  model  for  the  charring 
material  have  the  following  forms  (see  Ref.  2,  Appendix  A,  fox  justification) 


where  r  ranges  over  the  gas  generation  reactions  occurring  within  the 
material.  The  orientation  equation  for  the  generated  gas  mass  flux  in  the 
charring  material  takes  the  following  form: 
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Equation  5  is  a  formulation  of  assumption  5  above  for  charring  materials, 
vhich  in  conjunction  vith  assumption  4,  permits  the  solution  of  Bq.  1  and 
2  without  the  need  for  a  separate  momentum  equation.  The  dependent  vari¬ 
ables  of  Bj.  1  through  3  are  T,  G^,  and  (or  IL).  All  physical  and 
chemical  properties  appearing  are  assumed  to  be  known  functions  of  tempera¬ 
ture.  Bor  the  noncharring  vail  materials  present  in  a  throat  chamber  con¬ 
figuration,  Bq.  2  and  3  are  not  applicable,  and  Sq.  1  reduces  to  the 
pure  conduction  equation,  as  follovs: 
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In  the  event  of  surface  erosion  caused  by  engine  firing,  an  additional 
parameter  must  be  predicted,  f(x,r),  the  changing  radial  position  of  the 
receding  hot-gas  boundary  (Hg.  1).  The  following  equation  is  used  to 
relate  the  radial  recession  rate  to  the  predicted  normal  recession  rates; 
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where  r  ranges  over  those  species  generated  at  the  surface  and  does  not 

include  those  generated  within  the  material.  The  v  ’s  are  in  turn  defined 

r 

by  an  appropriate  rate  equation  (Ref.  2,  Appendix  C)  for  each  surface 
removal  mode  r.  In  practice,  the  vr's  are  calculated  iteratively  so  as 
to  satisfy  both  the  rate  equation  and  the  energy  balance  at  the  exposed 
inside  wall  surface,  which  is  taken  in  the  form; 
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More  generally,  Eq.  6  can  be  used  as  s  convenient  formulation  of  ail  the 
exterior  boundary  conditions  which  describe  the  heating  and  cooling  mech¬ 
anisms  encountered  (Pig.  l),  where  hgff  may  include  such  effects  as  those 
of  convection,  radiation,  and/or  aerodynamic  heating.  Thus  we  vrite 

eif  conv  rad  env  '  * ' 

At  an  insulated  boundary,  both  h^^  and  the  v  are  set  to  zero  in  Eq,  6. 

The  v^’s,  indeed,  apply  only  at  the  exposed  hot-gas  boundary  and  are  zero 
elsewhere.  Similarly,  T  ,  in  Eq,  6,  has  physical  significance  only  at  an 
exposed  boundary,  but  is  used  elsewhere  as  a  convenient  reference  temperature 
in  the  definition  of  the  effective  quantity  h^..  In  particular,  when  radi¬ 
ative  and  environmental  heat  flux  terms  are  required,  they  are  expressed 
as  follows: 


Thus,  the  effect  of  T  is  cancelled  when  h  .  and  h  are  substituted 

aw  rad  env 

in  Eq,  6  (this  form  is  employed  as  a  linearization  device  in  the  finite 
difference  solution  procedure).  The  cancellation  does  not  occur  in  the 
convection  tpjm  at  the  exposed  boundary,  where  ^conv-  includes  modification 
of  the  basic  convective  heat  transfer  coefficient,  h  ,  by  the  blocking 
effect  of  the  gases  ejected  at  the  surface: 
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The  samaation  of  Eq.  9  is  taken  over  all  gaseous  species  ejected  at  the 

surface,  where  G.  is  defined  to  be  ~pv.  for  the  gases  generated  at  the 
J  J 

surface,  and  for  the  internally  generated  gases  we  use 
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At  an  interface  separating  two  adjacent  material  regions,  say  regions  I  and 
J,  the  continuity  of  temperature  and  outward  normal  heat  flux  are  expressed 
as  follows: 
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Initially,  the  temperatures  are  assumed  to  be  constant  for  all  materials; 
e.g. ,  ambient  temperature. 


NUMERICAL  SOLUTION 

A  time  step  procedure  was  employed  to  uncouple  and  solve  the  continuous 
equations  given  above  where  finite  difference  techniques  were  used  to 
discretize  the  continuous  parameters.  Predicted  in  each  step  are  the 
two-dimensional  (actually  three,  because  of  the  assumption  of  axial 
symmetry)  temperature  distribution  throughout  the  thrust  chamber  materials, 
the  generated  gas  mass  flux  distribution  in  the  charring  material,  and  the 
new  position  of  the  inside  surface  exposed  to  the  hot  combustion  gases 
(in  the  event  of  surface  erosion  and  recession).  No  specific  char  front 
or  depth  is  calculated  because,  as  discussed  in  assumption  2  above,  py^  lysis 
is  permitted  to  occur  over  a  temperature  range  rather  than  at  a  single  value. 
The  pyrolysis  zone,  however,  can  be  identified  at  any  time  level  by  inspec¬ 
tion  of  the  predicted  temperatures  or,  in  the  event  of  a  stop-start  engine 
firing  cicle,  by  noting  tLe  maximum  temperatures  achieved. 

Spatially,  the  discretization  is  achieved  by  imposing  a  meah  on  the  multi¬ 
material  region  o?  interest  (Fig.  2).  Mesh  points  are  located  at  the 
intersections  of  the  mesh  lines  with  each  other  and  with  the  boundaries 
and  material  interfaces  and  are  classified  as  regular,  irregular,  boundary, 
and  interface  points  /Fig,  2).  The  continuous  temperature  and  generated 
gas  mass  flux  distributions  are  approximated  by  discrete  distributions 
defined  only  at  meah  points  aid  discrete  time  levels.  These  are  obtained 
by  solving  the  difference  equations  resulting  from  replacing  the  continuous 
derivatives  in  the  equations  above  with  their  finite  difference  equivalents 
(Ref,  2), 


9 


Figure  2.  Typical  Mesh  CondCiguratioa  oud  Categorization*  of  .Resulting  Me»h  Pointo 


Of  particular  interest  is  an  essentially  seaond  order  accurate  {in  tire 
as  veil  as  distance)  generalization  of  the  csconditi-mally  stable,  implicit 
alternating  direction  method  of  Peacesan  and  Sachford,  employed  in  each 
tine  step  to  discretize  end  solve  Iq.  1  (or  Eq,  4)  for  tesperatnre  in  con¬ 
junction  viih  second  order  accurate,  curvature— preserving  techniques  de¬ 
veloped  to  express  the  normal  gradient  conditions  at  curved  boundaries. 

A  backward  tise  difference  vas  used  to  approximate  the  continuity  equation 
(Eq.  2)  in  obtaining  the  internally  generated  gas  tags  flux  distribution 
along  viih  second  order  accurate  central  differences  in  the  axial  and  racial 
directions.  Although  a  forward  tise  difference  vas  used  to  discretize  the 
recession  equation  (Eq.  5),  the  approximation  nevertheless  approaches  a 
second-order  central  difference  procedure  because  of  the  iteration,  sta¬ 
tioned  above,  which  is  required  to  obtain  values  of  the  Details  of 
the  numerical  procedures  employed  are  given  in  Rex,  2. 

Because  of  the  use  of  essentially  central  (not  exactly  because  the  initially 
nonlinear  difference  equations  are  linearized  in  each  time  step  by  taking 
coefficients  at  the  old  time  level)  or  backward  time  differences  in  dis¬ 
cretizing  the  differential  equations  in  the  interior  of  the  material  regions, 
the  resulting  solution  procedure  approaches  unconditional  stability  for 
sufficiently  small,  externally  generated  heat  fluxes.  However,  when  high 
values  of  heat  flux  are  encountered  at  the  hot-gas  boundary  during  periods 
of  steady  firing,  limitations  may  be  required  on  the  sire  of  the  time  step, 
especially  at  the  beginning  of  such  a  period.  In  general,  such  limitations 
are  not  nearly  as  severe  as  those  required  for  explicit  forward  diffex*ence 
methods,  and  are  not  necessary  at  all  during  periods  of  sparse  intermittent 
firing  or  during  soakback.  The  favorable  stability  situation  plus  the 
generally  second  order  spatial  and  temporal  differencing  permit  the  use 
of  relatively  large  time  and  distance  increments  (compared  to  those  fre¬ 
quently  required  for  explicit  forward  difference  procedures)  and  a  sub¬ 
stantial  saving  in  computer  time. 

In  simulating  the  normal  temperature  gradient  conditions  at  curved  surfaces 
as  given  by  l’q.  6  and  11,  the  following  exact  relationships  are  utilized. 
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■roach  relate  the  norsal  gr»3  tangential  derivatives  to  these  in  the  axial 
sad  radial  directions: 


Sr  _  ±  /"  df  of?  of 
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In  Eg  12,  y  =  f  (x,  r)  is  the  equation  of  the  boundary  aegseni  under  con¬ 
sideration  and  the  plus  (sinus)  sign  is  used  at  a  lever  (upper)  boundary; 
i-e.,  where  the  outward  nornal  has  a  negative  (positive)  component  in  the 
y  direction,  uhe  direction  of  n  and  *  is  indicated  in  Fig.  1  for  a  mat- 
erial  at  a  boundary  and  an  interface.  Hie  expressions  for  df/dn  and  Bf/fcs 
can  be  combined  (in  keeping  with  the  alternating  direction  method)  bo  as 
to  eliminate  Sf/3y  (or  3f/^x)  in  an  odd  (even)  tine  step.  Finally,  on 
purely  physical  grounds,  vs  assume  that  df/d3  is  sc  dominated  by  St/Sr 
that  in  -the  resulting  expressions,  ~ie  say  neglect  terms  in  ^f/^3.  Thus 
ve  obtain,  for  odd  and  even  steps,  respectively: 


a  .  2^/2  (*3) 

!-!/[-(!  f\ 

Eq.  13  is  need  to  replace  3T/dn  in  Eq.  6  and  11  and  the  resulting  expres¬ 
sions  are  discretized  vith  second  order  spatial  accuracy. 

All  systems  of  difference  equations  generated  in  each  time  step  ere  at 
worst  tridiagonal*  and  are  easily  solved,  directly  (rather  than  iteratively) 
by  simple  recursion  formulas. 


mi  i—i— Wym—  >***  tmrmmmm* 

*A  system  of  equations  is  tridiagonal  if  all  the  nonnero  elements  of  its 
coefficient  matrix  lie  on  the  main  diagonal  and  its  ineaediately  adjacent 
diagonals. 
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RESULTS 

Evaluation  of  the  2B-ABLASE  program  included  comparisons  of  program  results 
vitfc  both  theoretical  and  test  data,  including  in-depth  checkout  of  the 
following  cases:  "• 

1.  Simulation  of  transient  one-dinensional  radial  conduction  in  a 
hollow  cylinder  heated  cn  the  inside  and  insulated  on  the  outside, 
with  no  charring  or  recession.  The  exact  series  solution  for  this 
problem  is  vell-kncvn  and  has  been  presented  in  the  fora  of  tem¬ 
perature  response  curves  in  Eef .  6. 

2.  Caparison  to  test  data  supplied  by  Edwards  AFB  from  an  engine 
firing  with  shutdown  and  soakback  in  a  two-material  thrust  cham¬ 
ber,  including  a  charring  carbon-cloth/phenolic  backed  by  a  stainless- 
steel  shell  (Fig.  3).  Here  the  full  two-dimensional  transient 

case  was  treated. 

3.  Comparison  to  measured  temperatures  obtained  from  a  complicated 
mission  duty  cycle  of  engine  firing  (including  pulsing)  and  aero¬ 
dynamic  heating  for  a  four -material  charring,  attitude  control  engine. 
Predicted  computer  results  were  m  good  agreement  with  measured 
backvali  throat- ana  exit  temperature  histories.  (This  comparison 

was  not  reported  in  Ref.  2;  details  may  be  fount?  m  Ref.  7.) 

Cases  1  and  2  above  are  fully  reported  in  Ref.  2  and  need  not  be  repeated 
here.  Some  of  the  results  of  Case  2,  however,  were  used  as  controls  for 
checking  out  the  2D-ABLATE  program  extensions  discussed  later  in  this 
report. 


PROGRAM  STRUCTURE 

The  2D-ABLATE  program  was  coded  in  Fortran  IV  for  use  on  the  IBM  7094. 
Several  links  were  employed  to  take  advantage  of  the  system  overlay  feature. 
In  addition  to  the  main  control  link,  Link  0,  two  other  links  were  provided 
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1.  Link  I,  is  viich  all  sefc-cp  is  parforoed  tor  ihs  iise-sie s  cal¬ 
culation  (Sef.  2) 

2.  Link  2.  tie  tire— step  calculation 

Because  3ie  extensions  to  tie  22MB^B  gregraa  (rsder  tie  cmeal  con¬ 
tract)  expanded  Link  2  beyond  the  7094  core  storage  capacity,  T.vrgr  2 
been  further  sublinked,  utilising  the  overlay  feature,  ghis  is  described 
in  sore  detail  in  another  section  of  this  report. 
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In  Task  *  of  the  present  contract,  the  2D-AELAIZ  caa&uter  program  (as 
developed  cnder  Contract  AHi4(6ll)-9714)  vaa  extended  by  tse  addition  of 
several  mechanisms  each  of  which,  when  applicable,  can  have  a  significant 
effect  on  the  performance  of  ablative  materials.  Included  vere  the  xol— 
loving  ncdi fixations: 

A-  The  effect  of  radiative  exchange  across  the  thrust  chamber  betveen 
opposing  surfaces  exposed  to  the  hot  combustion  gases.  During 
firing,  this  effect  can  exert  a  strong  influence  on  the  surface 
temperature  histories  attained. 

B.  Analysis  of  materials  vith  differing  temperature-dependent  con¬ 
ductivities  in  major  and  minor  directions  to  provide  adequate 
treatment  for  strongly  anisotropic  materials  oi  current  interest. 

C.  Prevision  to  permit  any  of  the  vail  materials  to  be  charring 
ablators  rather  than  being  restricted  to  just  one  such  material. 

11  ^ oKi« ^  /10+ tri  + nvilv  »  <?■»«»!  a 
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charring  material,  ®»re  than  one  is  often  required  in  the  case 
of  production  engines.  Frequently  encountered  classes  of 
multiple-ablative  construction  include  stacked  materials  (e.g. , 
the  Bockeidyne  LEH  engine  vith  a  phenolic-refrasil  laminate 
backed  by  a  lov  density  phenolic-asbestos)  and  side  by  side 
(e.g. ,  the  Apollo  attitude  control  engines  vith  a  45-degree 
phenolic-refrasil  laminate  in  the  chamber  and  the  same  material 
in  the  nozzle  extension  at  a  0-degree  orientation). 

In  addition  to  the  three  Task  1  modifications  (henceforth  referred  to  as 
Tasks  l.A,  l.B,  and  l.C),  several  other  additions,  modifications,  and  cox  - 
rectione  were  made  to  the  2D-ABLAIE  program  as  required  during  the  perform¬ 
ance  of  Tasks  1;  2,  and  3  of  the  current  study.  Because  these  changes 
do  not  fall  within  the  scope  oi  Task  1,  they  are  referred  to  collectively 
as  "general  modifications"  and  briefly  described  prior  to  the  Task  J 
discussion. 


Checkout  of  most  of  these  Bonifications  vas  performed  by  comparison  with 
results  obtained  previously  with  the  "old  version"  of  2D-AELATE.  The  test 
case  used  for  the  comparisons  vas  Case  2  discussed  in  the  Eeview  section 
above  and  described  in  detail  in  Ref.  2.  For  Task  1,  a  preliminary  check¬ 
out  step  involved  demonstration  of  the  ability  of  the  extended  program  to 
suppress  all  the  new  mechanisms  and  perform  as  originally  written. 


GENERAL-  MODIFICATIONS 

Program  Structure 

As  indicated  in  the  Review  section  above,  coding  of  the  old  version  of  the 
2D-ABLATE  program  was  performed  for  IBM  7094  application.  The  "new  version 
(including  the  Task  1  extensions)  was  to  be  nsed  on  either  the  7094  or  the 
IBM  360.  To  accommodate  the  Task  1  program  extensions,  several  methods 
were  investigated  of  increasing  available  core  storage  in  the  new  7094 
version,  (in  the  new  360  version,  no  such  problem  existed.  Indeed,  it  was 
possible  to  unlink  the  program  completely  without  exceeding  core  storage 
limitations.)  It  was  decided  to  further  subdivide  the  time  step  link 
(Link  2  in  the  old  version),  yielding  Link  2,  the  control  link,  Link  3, 
the  recession  calculation,  and  Link  4,  the  remainder  of  the  calculation 
in  each  time  step.  Thus,  Link  3  is  bypassed  in  any  time  step  for  which 
the  recession  subroutines  are  not  required  (see  Fig.  4  for  a  chart  of  the 
relinked  program  structure).  The  saving  amounted  to  over  3500  locations 
in  core  during  Link  4  (and  ever  more  during  Link  5),  which  was  more  than 
enough  for  accommodation  of  the  Task  J  extensions. 


Added  Capability 


Because  performance  of  the  Task  2  and  3  efforts  required  extensive  evalua¬ 
tion  of  calculated  temperature  histories,  a  great  deal  of  time-consuming 
manual  cross-plotting  was  avoided  by  extending  the  program  to  provide  a 


(a).  Linkage  of  the  nev  version 


Link  2 


(b)  Program  flow  in  Link  2 


Figure  4.  Structure  of  the  Relinked  Program 
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visual  display  open  inpat  signal  of  as  many  as  10  temperature- time  graphs* 
each  at  a  distinct  specified  regular  interior  point  (see  Fig.  2  for  ter®- 
inology)  within  the  vail  materials.  This  vas  accomplished  by  taking  ad¬ 
vantage  of  the  available  cathode  ray  tube  equipment  and  the  C2T  library 
subroutines . 

Improvement  in  Calculation  of  Gas  Generation 

Improvements  were  made  in  the  portion  of  the  program  dealing  with  internal 
charring.  First  of  all,  the  simulated  boundary  and  initial  conditions 
(pages  64  and  65,  Bef.  2)  for  the  solution  of  the  continuity  equation  were 
altered  to  reduce  initial  lag  in  gas  generation.  In  particular,  the  simu¬ 
lated  boundary  condition  c,  page  65,  Bef.  2,  should  now  read,  in  part, 

(see  Fig.  5  for  the  point  arrangement  used  te  discretize  the  continuity 
equation)  "If  ^  i8  less  than  T  or  if  point  4  is  a  noncharring 

material,  then  G„  is  set  equal  to  0  "  rather  than  "If  either  T  0  or 

^max  4  *s  *ess  ^an  ^py““"  Secondly,  because  the  temperature  may  jump 
significantly  into  the  pyrolysis  range  during  a  relatively  large  time 
step  at  mesh  points  close  to  an  exposed  surface,  an  upper  limit  for  gas 
mass  flux  within  a  charring  material  can  now  be  assigned  as  input  to  the 
program  to  prevent  excessive  initial  interior  gas  generation  and  possible 
temperatu  oscillation  caused  by  overcooling. 

These  program  changes  were  successfully  checked  out  using  the  test  case 
discussed  above.  As  expected,  the  result  was  a  graduated  temperature  re¬ 
duction  throughout  the  thrust  chamber,  the  most  significant  reductions 
occurring  at  the  heated  inside  surface. 


Correction  of  Errors 


In  the  cor  ae  of  performing  the  Task  2  and  3  efforts,  a  number  of  errors 
were  uncovered.  Host  were  of  a  strictly  "prograraning"  nature  such  as 
spelling  and  indexing  of  program  variables.  Correction  of  these  errors 
did  not  alter  results  of  previously  run  checkout  cases  because  of  the 
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markedly  different  material  and  msak  line  configurations.  On  the  other 
hand,  several  other  errors  were  discovered  in  the  naraerie&l  interpretation 
of  ihe  model  and  vere  subsequently  corrected,  They  vers  as  follows. 

Inspection  of  Eq.  9  above  indicates  that  tfee  values  of  G .  used  in  Es.  9 
vere  intended  to  be  negative  for  all  gaseous  species  ejected  at  the  surface, 
as  in  Eq.  10.  In  the  progress,  hcveves*,  the  pins  sign  in  Bq>  9  was  inad¬ 
vertently  written  as  a  minus.  Thus,  the  transpiration  effect  of  the  in¬ 
ternally  generated  gases  was  reversed.  For  the  checkout  ease  discussed 
in  the  fieview  aection,  the  effect  of  this  error  proved  to  be  not  overly 
significant. 


Errors  were  found  and  corrected  in  the  subroutines  which  predict  erosion 
and  recession  at  varying  axial  positions  on  the  exposed  inner  surface  and 
in  the  temperature  calculation  at  these  points  after  erosion.  In  the 
recession  calculation,  too  few  restrictions  were  placed,  on  the  number  of 
loops  of  iteration  permitted  and  on  the  number  of  circumstances  under 
which  the  iterations  would  be  permitted  to  continue,  both  in  the  calcula¬ 
tion  of  Mach  number  and  in  the  calculation  of  erosion  rate.  In  the  cal¬ 
culation  of  temperature  after  erosion  haa  occurred,  an  error  was  found  in 
the  interpretation  of  the  information  supplied  by  the  recession  subroutines. 
In  particular,  an  attempt  had  been  made  to  bypass  the  direct  simulation 
of  Eq.  6  after  erosion  by  employment  of  the  value  of  Si'/dn  as  calculated 
iteratively  in  the  recession  subroutines.  The  nature  of  the  iteration, 
however,  is  .«ueh  that  Sl/Bn  is  evaluated  at  some  intermediate  time  level 
between  the  old  and  new  levels.  The  value  required  for  the  temperature 
calculation,  on  the  other  hand,  must  be  obtained  at  the  new  level.  Con¬ 
sequently,  the  short-cat  procedure  was  abandoned  and,  Eq.  6  has  been 
simulated  directly  in  the  program. 


*  2; 


A21eviati.cn  of  a  Program  Limitation 


in  several  of  the  thrust  cheater  configurations  treated  during  Task  5,  the 
nozzle  extension  vas  taken  to  be  a  long  narrov  slanting  strip  of  vali  sa- 
tenal.  Because  of  the  ret  .resent  of  equal  spacing  of  the  horizontal 
(i.e. ,  axially  oriented)  nesh  lines,  aid  the  neper  lisit  of  15  on  the 
aroaber  of  vertical  mesh  lines,  it  vas  necessary,  in  order  to  have  a  suf¬ 
ficient  ucsther  of  horizontal  mesh  lines  for  accuracy  in  the  chzsher  and 
throat  regions,  to  force  toe  uany  into  the  nsrrcv  strip  of  extension  sa- 
terial,  The  result  vas  a  violation,  on  several  horizontal  sash  lines  in 
the  strip,  of  a  program  liniiaticr  on  the  sinirnxa  ctszher  of  points  per¬ 
mitted  per  sash  line  (four,  as  gives  in  Bale  4  on  page  79  of  Bef.  2).  As 
a  consequence,  a  duroy  subroutine  vas  written  which  can  be  used  -when  re¬ 
quired  to  replace  the  restrictive  subroutine  (called  MIH4PT)  and  thus 
reduce  the  ninisas  iron  four  to  three  (the  latter  being  at  present  an  ir¬ 
reducible  sin ;  zee  due  to  other  program  lici tat ions  as  expressed  by  Boles 
1  and  5  on  pages,  77  through  79,  Bef.  2).  lie  resale,  of  relaxing  the  four 
point  1  ini  tat  ion  vas,  as  expected,  a  worsening  of  the  ‘-condition"  of  the 
systems  of  difference  equations  generated  in  the  affected  portions  of  the 
strip  and,  in  particular,  a  certain  amount  of  error  in  the  calculated 
temperature  in  odd  time  steps  at  the  intersections  of  the  radiation-cooled 
outer  boundary  of  the  nozzle  extension  viib  some  of  the  horizontal  mesh 
lines.  Xo  transmitted  ill  effects,  hovsver,  were  observed  (again,  as 
expected)  in  the  throat  and  chamber  regions,  which  were  of  primary  interest 
in  the  parametric  study  performed  in  Task  3- 


TASK  1,  A,  RADIATIVE  EXCHANGE  AT  THE  EXPOSED  INSIDE  SURFACE 


The  effect  of  radiative  exchange  between  axially  varying  portions  of  the 
inside  vail  surface  vas  added  as  an  optional  (upon  input  signal)  term  of 
the  effective  heat  transfer  coefficient  (see  Eq.  7  and  8  above),  as  follows 

,  _  ^rerad 

rerad  T  -  T 


The  approach  taken  in  tbs  program  is  to  recalculate  values  of 
the  reradiative  heat  flirt,  after  each  even  time  step  at  the  boundary 
points  P^  forced  by  the  intersections  of  the  radially  oriented  mesh  lines f 
x  =  x. ,  i  =  n,  with  tbe-  exposed  inside  surface,  as  follows: 


Vr«d,i-Z  0  f  Fi;k  (fk  '  €>• 


In  Eq.  15,  F.  ,  ie  tbe  viev  factor*  from  an  inside  surface  area  represented 
by  P.  to  a  similar  area  containing  P^.  Include  is  the  term  with  viev 
factor  F.  .  because  tbe  inside  surface  area  is  ir  the  shape  of  a  rnng  and 
opposes  itself  as  veil  as  the  remainder  of  the  rings.  The  array  obtained 
using  Eq.  15  is  then  used  for  the  succeeding  two  time  steps  in  Eq.  14  to 
calculate  the  reradiative  component  of  the  effective  heat  transfer  coef¬ 
ficient.  In  an  odd  step,  values  of  q^  ^  at .boundary  points  lying  at 
tbe  ends  of  axially  oriented  mesh  lines  are  obtained  as  needed  by  inter¬ 
polation  from  the  ^  array. 

Values  of  the  two-dimensional  array  F.  .  are  obtained  by  interpolation 

1  ,K 

at  the  points  P..  from  a  similar  input  array  G^.  ^  defined  at  the  input 
points  Q.  at  which  the  rest  of  the  axially  varying  input  data  are  speci- 
ficd,  such  as  the  adiabatic  wall  temperature  and  the  convective  heat  trans¬ 
fer  coefficient.  This  is  performed  in  three  stages  in  Link  1  of  the  pro¬ 
gram  (the  setup  phase  prior  to  the  time  step  calculation),  as  follows; 

1.  Input  and  assignment  of  the  two-dimensional  array  G^  j,  -t  =  1, 

...,  m,  at  the  points  0.  and  of  a  one-dimensional  array  B.  of 

J  J 

associated  disjoint  surface  subareas. 


*Tb.s  use  of  viev  factors  in  Eq.  15  instead  of  overall  interchange  factors 
(P.ef.  8  and  9)>  which  would  include  the  effect  of  reflection,  follows  from 
the  analysis  in  Tasks  2  and  3  of  vail  materials  which  emit  nearly  as  black 
bodies.  The  net  reflected  heat  fluxes  are  nearly  zero  and  the  inter¬ 
change  factors  nearly  equal  to  the  viev  factors.  Modification  of  the  pro¬ 
gram  to  include  reflection  effects  would  be  straightforward. 


2k 


2V  Direct  calculation  of  a  similar  array  A.,  i  »  1,  >,.j  n,  of  snb- 
areas  associated  with  the  P^,  each  assumed  conical,  i.e, , 

Ai  =  "  Z1  +  ^  'Zi+l/2  “  Zi-l/2^  ^fl  ^Xi+l/2+Xi-l/2^ 


2(f ,-x.f!  )  J 
'  1  x  r' 


(16) 


vhere  y  =  f(x)  is  the  equation  of  the  inside  boundary  curve  and 

fi  =  f^xi)’  811,5  vhere  xi+l/2  =^xi  +  xi+l)/^  defines  the  end  points 
of  the  intervals  over  which  the  arc  calculated.  The  generated 


array  A,.  is  then  normalized  as  follows: 


m 


n 


\  *  Ai  l  V  l  V 


(17) 


j=l 


k=l 


so  as  to  satisfy  the  following  normalization  condition  (where, 
for  convenience  of  notation,  we  drop  the  bar  on  the  Aj): 


n 


m 


1  A»  -l  h- 


(18) 


i«l 


j=l 


Finally  a  two-sweep  linear  interpolation  is  performed  to  obtain 
the  F.  ,  at  the  paints  P.  from  the  G.  .  at  the  points  Q.  in  such 
a  way  as  to  satisfy  the  following  two  normalization  conditions: 


m 


l  nj,k  ■  l  V  J-1- 

4=1 


m. 


(19) 


k=l 


m 


1  AiFi,k  =  1  BjHj,k>  k=1>  ,4‘ 

j=l 


(20) 


i-1 


mis  definition  of  is  applicable  only  for  i  -  1,  2, . . . ,  n-I.  x,  /„ 

^+1/2  are  taken  to' he  the  leftmost  and  rightmost  axial  positions,  ' 
respectively,  of  the  exposed  inside  surface. 


vhere  E.  ,  ,  j=l,  ...»  ra,  k=l,  . ..,  n,  is  an  interin  two-dimensional 
array  of  view  factors,  from  the  points  Q.  to  the  points  P.,  oh- 
tained  daring  the  first  sweep  of  interpolation. 

To  make  use  of  the  added  Task  l.A  program  capability,  it  is  necessary  to 

obtain  reasonable  values  for  the  input  view  factors  G.  j  and  associated 
)  3 1 

^-^Jareas  for  each  contemplated  thrust  chamber  configuration.  For  check¬ 
out.  purposes,  an  extension  of  the  disk  method  developed  by  Br.  Simon  deSoto 
of  Rocketdyne  (Ref.  10)  was  used  to  generate  the  G^  ^  and  is  described  in 
some  detail  in  Appendix  A.  The  disk  method  is  exact  for  the  case  of  conical 
subareas,  none  of  which  are  occulted  or  shaded  from  each  other.  Because 
the  occultations  in  the  checkout  case  were  not  severe  (Fig.  3),  the  formulas 
were  used  as  if  no  occultations  existed.  No  rigorous  method  was  found  in 
the  literature  which  would  yield  view  factors  from  all  portions  of  a 
converging-diverging  thrust  chamber  to  each  other  portion,  including  the 
effects  of  occultation  and  shading.* 


The  input  subareas  B.  employed  in  the  test  case  used  for  checkout  were  ob- 

J 

tained  through  application  of  Eq.  16  at  points  Q.  rather  than  P.  .**Values 

J 

of  the  B.  and  G.  ,  arrays  used  for  the  Task  l.A  checkout  are  given  in 

J  J  J 

Table  1,  and  the  arrangement  of  the  subareas  and  their  approximation  by 
conical  rings  are  depicted  in  Fig.  6.  The  other  properties  and  heating 
conditions  used  for  the  checkout  comparison  can  he  found  in  Table  4  of 
Ref.  2,  The  checkout  run  made  with  reradiation  was  terminated  after  12 
seconds  of  steady  firing.  Comparisons  Of  computed  surface  temperature 
histories  with  the  previously  obtained  results  are  given  in  Fig. 7  and  8.  As 
expected,  temperatures  computed  with  reradiation  effects  included  were 
lower  in  the  vicinity  of  the  throat  and  higher  in  the  chamber  and  exit 
sections. 


*A  rigorous  procedure  is  presently  being  developed  at  Rocketdyne.  When 
available,  it  would  be  desirable  to  program  the  method  and  incorporate  it 
as  a  subroutine  of  Link  1  in  2D-ABLATE.  This  would  eliminate  a  great  deal 
of  arduous  computation,  necessary,  for  example j  for  application  of  the  disk 
method  using  a  desk  calculator. 

■^Selection  of  the  end  points  21+1/2  (corresponding  to  the  xj,+i/2  Eq.l6)  is 
not  restricted  to  midpoints  between  the  xj.  Instead,  because  judgement  can 
be  more  readily  exercised  in  a  band  calculation,  they  can  be  taken  at  any 
p  'nt  between  the  Xj,  Generally,  they  should  be  choaen  on  a  geometrical 
basis  so  as  to  yie’a  subsurfaces  which  are  most  nearly  conical  in  shape. 
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TASK  l.B,  THERMAL  ANALYSIS  OF  ANISOTROPIC  MATERIALS 


The  2B-ABLATE  program  was  extended  to  permit  (upon  input  signal)  any  of 
the  vail  materials  to  he  anisotropic;  i.e.,  to  possess  two  different  con¬ 
ductivities,  each  temperature  dependent,  one  in  the  major  direction  (highest 
conductivity)  and  the  other  in  the  minor  direction  (lowest  conductivity)  per¬ 
pendicular  to  the  major  direction.  Many  wall  materials  of  current  interest 
are  strongly  anisotropic  such  as  phenolic-refrasil  and  phenolic-carbon 
cloth  laminates,  which  were  analyzed  in  some  detail  in  Task  2  and  3  of 
the  present  program. 

For  treatment  of  anisotropic  materials,  modifications  were  necessary  in 
the  formulation  and  discretization  of  the  model  as  expressed  by  Eq.  1 
through  11  given  above.  The  overall  solution  procedures,  however,  were 
substantially  unchanged.  A  discussion  follows  of  the  changes  made  and 
their  effect  on  program  results. 


Model  Changes 


For  an  anisotropic  wall  material  with  conductivities  K^(T)  and  K^(T)  in 
the  major  (rj)  and  minor  (£)  directions  and  a  counterclockwise  displacement 
0* of  these  directions  from  the  radial  (y)  and  axial  (x),  the  following 
version  of  the  energy  equation  was  derived  (see  Appendix  B) ,  in  which, 
for  generality,  both  conduction  and  charring  are  accounted  for: 


c>T 


(pchff  St  =  k  (Kx  If)  +  7 1?  frKy  sf> 


dH  /„  dT  _  dTv 
dT  'Gx  Sx  +  Gy  Sy>  + 


2K 


d2T 

xy  Sx5^ 


+  2 


dT  dT 


(21) 


*It  should  be  emphasized  that  0  is  not  what  is  commonly  termed  the  orienta¬ 
tion  (or  angle  of  wrap),  the  latter  conventionally  being  taken  as  the  clock¬ 
wise  rotation  from  the  axial  to  the  major  direction. 


31 


where 


(22) 


and  where  (PC)ef£  is  nsed  to  abbreviate  the  coefficient  of  dT/d7  on  the 
left  side  of  Eq.  i. 

For  use  in  the  boundary  and  interface  heat  flux  conditions,  given  by  Eq.  6 
and  11  above,  the  following  expression  wa.B  obtained  (see  Appendix  B)  for 
the  conductivity  in  the  normal  direction  n  at  the  bounding  surface: 

.  Kn  ‘  Kr  '  2  I  V  *  K 1  +  (H>2]  (23) 

where  df/dx  is  the  slope  of  the  boundary  curve,  y=f (x , t) . 

No  direct  changes  were  required  for  the  anisotropic  analysis  in  Eq.  2, 
the  continuity  equation  for  the  mass  flux  of  generated  gases  within  a 
charring  material.  The  direction  of  flow,  however,  in  an  anisotropic 
material,  as  characterized  by  the  ratio  of  G  to  G  ,  is  now  assumed  to  be 

*  y 

oriented  with  the  heat  flux  vector,  as  follows*: 

°A  -  **  H 

which  is  seen  to  be  compatible  with  Eq.  3  above  for  an  isotropic  material; 
i.e. ,  when  Kx  =  K^.,  It  can  be  shown  that  Eq.  24  closely  approximates  the 
expression  G^/G^  =(k^  dT/d^^dl/d  Tj). 


*An  alternative  assumption,  which  was  not  employed,  wo old  have  taken  the 

direction  of  flow  to  be  in  the  major  direction  t)  of  conductivity;  i.e,, 

G  /G  =-tan  0,  It  was  felt,  however ,  ■  that  in  general  Bq.  24  would  be 
x  y 

more  appropriate. 
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Discretization  of. 


The  numerical  solution  of  the  model  equations  for  the  case  of  an  aniso¬ 
tropic  material  vas  obtained  with  a  time  step  procedure  using  finite  dif¬ 
ferences  as  in  the  isotropic  case  outlined  above  and  detailed  in  Ref.  2. 

Discretisation  of  the  anisotropic  energy  equation,  however,  required  def¬ 
inition  of  a  finite  difference  analog  of  the  aired  derivative  3^T/3xdy. 

For  the  remaining  derivatives  of  Eq.  21,  the  centered  first  and  second 
spatial  difference  operators,  6  ,  6 ^  ,  6^.,  ,  and  the  first  time  dif- 

ference  were  employed  as  defined  by  Eq,  19  through  23  of  Ref.  2, 

At  an  interior  point  0  whose  horizontally,  vertically,  and  diagonally  ad¬ 
jacent  points  are  all  regular  (see  Fig.  2  for  definition  of  regular  point) 
and  where  the  points  are  equidistant  (as  in  Fig.  9a  for  &  x.  -  Ax„  =■  Ax 
^  yjg  ~  &  y p  =  Ay),  the  following  standard  mired  difference  analog 
could  have  been  used  with  second  order  accuracy  (See  App-  .dir  C  for 
derivation): 

6^  Tq  -  (Te  -  Tp  +  Tg  -  Th)  /  4  Ax  Ay  (25) 

In  the  case  of  nonequal  spacing  (as  occurs  in  2D-ABUVTE  in  the  x-direction), 
the  following  similar  analog  could  also  have  been  applied: 

^ay  -T0  =  “  TF  +  TG  ”  ^  ZA  +  ^  XC^  ^  yB  +  ^yD^ 

but  only  with  first  order  accuracy;  i.e.,  with  0(Ax^  -  A  x^)  +  o(AyE  - 
A  y^)  accuracy  (which  improves  to  second  order  accuracy  whenihe  points  are 
equidistant  as  in  Eq.  25  and  thus  might  be  termed  accurate  of  order  3/2 
because  it  ia  effectively  between  first  and  second).  This  would  have  been 
quite  satisfactory  in  the  approximation  of  the  energy  equation  because  of 
the  similar  truncation  error  obtained  in  the  discretization  of  the  second 
spatial  derivatives.  Neither  Eq.  25  nor  26,  however,  are  employed  in  2D- 
ABLATE  because  they  too  often  are  inapplicable;  i.e.,  if,  for  a  given 
interior  point  0,  one  or  v,ore  of  the  points,  F,  F,  G,  and  H  do  not  ij.e  in 
the  same  material  because  of  interposition  of  boundaries  or  interlaces 
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{«3  jq  Fig.  9b,  for  ezaaple).  bhai  is  required  instead  in  such  a  case 
is  a  difference  analog  involving  temperatures  at  the  intercepted  irregular 
boundary  or  interface  points  (such  as  F^,  F,,,  Gv ,  and  Gg  in  Fig.  9b).  On 
the  other  bend,  to  formulate  e  separate  difference  analog  for  each  8u.»h 
combination  of  irregular  points  (Fig.  9b  depicting  just  one  possible  com¬ 
bination)  would  be  a  difficult  task.  Instead,  the  following  second  order 

accurate  sized  difference  analogs  were  derived  (see  Appendix  C) ,  one  for 
odd  and  one  for  even  time  steps,  which  treat  all  cases  of  the  type  shown 
in  Fig.  9b  because  of  their  formulation  in  terms  of  temperature  gradients 
at  the  horizontally  and  vertically  adjacent  points  instead  of  temperatures 
at  the  diagonal  points. 


m 


Odd  Step: 


4,  To  -  L4  4  Va  -  *4 > v v*  4  Vc]/ 

t— 

L*xa  ^xc  ^2a  +  Ax(p] 


Even  Step: 


&L  T 
xy  0 


*  4  6a  -  4 


-  *4  >  6,T0  -  *4  \\lf 


(27) 


^yD  ^yB  + 


-  i 


ft-  -  I 


Using  the  definitions  given  in  Ref.  2  of  the  second  order  accurate  differ¬ 
ence  operators,  6^  and  £y  it  can  be  seen  that  if,  as  in  Fig.  9a,  all  of 
the  adjacent  points  B,  C,  D,  E,  F,  G,  and  H  are  regular  and  in  the 
sane  material,  then  Eq.  27  and  28  are  identical  (hut  do  not  reduce  to  the 
less  accurate  first  order  analog  expressed  by  Eq.  26).  If,  further,  the 
spacing  is  equal,  both  reduce  to  Eq.  25. 

Equations  2?  and  28  are  emploved  in  the  2D-ABLATE  program  for  the  discre¬ 
tization  of  the  anisotropic  energy  equation  at  each  interior  point,  0  whose 
horizontal ly  and  vertically  adjacent  points  A,  B,  C,  and  B  are  also  interior 
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points.  If,  however,  any  of  these  points  are  boundary  or  interface  points 
(as  in  Pig.  9c,  for  example),  then  either  Eq  27  or  28  or  both  are  not  ap¬ 
plicable  directly.  This  is  caused  by  the  difficulty  encountered  in  expres¬ 
sing  the  first  differences  5sTB,  5^,  6yTA,  or  6yTc  when  the  points  B,  2), 
A,  or  C  are  irregular.  The  same  problem  arose  (described  in  Bef.  2,  pages 

47  through  53>  and  indicated  in  the  Review  section  above)  in  obtaining 
a  second  order  accurate  difference  analog  to  5t/5q  for  use  in  boundary  and 
interface  conditions  at  irregular'  points.  Here  we  do  not  employ  the  same 
device  used  for  maintaining  second, order  accuracy;  i.e.,  permitting  con¬ 
version  from  df/dx  to  ^T/Sy  at  points  B  and  B  or  from  &T/dy  to  St/3x  at 
A  and  C  by  essentially  letting  St/Ss  vanish  and  using  the  equivalent  con¬ 
dition,  dT/dx  =  -df/Bx  dT/3y,  for  the  conversion.  Instead  we  revert  to 
first  or  even  zero  order  accuracy,  in  the  event  of  irregular  points,  by 
replacing  Eq.  27  and  28,  as  follows: 


Odd  Time  Step: 


-  Vo 


62  J  -l 
uxy  Ao  -  \ 


6jf_  -  5L 
0  V  y 


v. 


if  A  is  regular  and  C  is  irregular 


it  0  is  regular  and  A  is  irregular  (29) 


if  A  and  C  are  irregular 


Even  Time  Step; 


if  B  is  regular  and  D  is  irregular 


if  D  is  regular  and  B  is  irregular 


V. 


0 


if  B  and  D  are  irregular 


(30) 
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Piscretizatita  of  the  imergy  Equation 

The  discretization  of  the  anisotropic  energy  equation  wag  performed  as  in 
the  isotropic  case  in  accordance  with  the  linearization  of  the  generalized 
Peaceaan-Racfcford  alternating  direction  method  as  expressed  by  Eq.  30 
through  32,  Ref.  2.  Tridiagonality  of  the  systems  of  difference  equations 
generated  was  maintained  by  evaluating  the  mixed  difference  analogs,  as 
given  by  Eq.  27  through  30,  at  the  old  time  level  (i.e..  at  time  in  the 
time  step,  Tk  j)  rather  than  the  new.  Similarly,  discretization  of 
the  term  2  K^.  (T)  St/3x  df/dy  was  performed  so  as  to  maintain  linearity 
and  tridiagonality  of  the  difference  equations  by  taking  6  T  at  the  old 

y 

time  level  and  6^T  at  the  new  in  an  odd  time  step  and  reversing  the  assign¬ 
ment  in  an  even  step.  Thus,  in  odd  and  even  time  stepa,  the  anisotropic 
difference  equations  take  the  following  form  (where,  following  the  nota¬ 
tion  cfRai.  2,  Eq.  31  and  32,  the  spatial  subscript  0  is  omitted  but  is 
understood  to  apply  to  each  parameter  appearing) : 


Odd  Step  (rk  -  rk+1):. 
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x.k  5 


eff,k  ~Et 

Kiy,k  Vk)  Vk+1  +  Ky,k  K;,k  «y\  - 
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sr  <L  +  K<  6  T,_  )  6  T,_  + 


yTk  >,k  '  “xy,k  wx‘k  ''  Vk  "  2Kxy,k  \y  Tk 


(31) 
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Even  Step  (T^  -  rk+2) 


'gff,k+l 


T  -  T 

k+2  k+1 


x,k~l  6x  Tk+1  +  kfl  Jx\~l  ~ 


6e 

Gx,k+1  *  Kxy,  k+1  5yTk+l^  6xTk+l  +  Ky,k+1  °y  Tk+2  + 


Vj'  A  y  +  1{y,  k+1  5y"k+l  “  ^'k+1  Gy,k+1  r  Ksy,  k+1  5xTk+L^  Vk 


2  Kxy,  k+1  5xy  Tk+i 


The  notation  used  in  subscripting  and  differencing  the  parameters,  which 
follows  that  of  Ref.  2,  Eq.  31  and  32,  should  be  self-evident.  In  parti¬ 
cular,  and  6  \+i/^yT{j+i  are  use^  iaatead  of  the  known  values 

H1  (l’k)  and  H1  (Tjj+1)>  as  explained  in  Ref.  2,  to  ensure  inclusion  of  all 
heat  absorption  effects  when  the  distance  increments  used  are  so  large 
that  entire  temperature  intervals  might  be  passed  over  in  which  H1  (T)  is 

relatively  large.  Also,  j  £>y  »  y.,  the  height  of  the  j—  axial  mesh  line. 
*  J 


Revised  Programming  Procedures 


T".  s  2D-ABLATE  program  was  modified  for  treatment  of  anisotropic  materials 
in  the  input  and  setup  procedures  of  Link  1  and  in  the  time  step  calculation 
performed  in  Links  2,  3,  and  4  (see  the  Review  section). 


Upon  input  signal,  any  of  the  five  wall  materials  can  be  anisotropic.  For 
each  anisotropic  material,  the  input  data  include  the. major  and  minor  con¬ 
ductivities,  K^(T)  and  K^(T),  read  in  as  piecewise  quadratic  functions  of 
temperature,  and  the  counterclockwise  angular  displacement  9  of  the 
and  £  directions  from  y  sad  x.  The  anisotropic  input  in  Link  1  is  then 
transformed  to  equivalent  conductivities  in  the  axial,  radial,  and  "mixed" 
directions  for  later  use  in  Links  2,  3,  and  4.  The  transformation  is 
accomplished  by  use  of  Eq.  22. 
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The  time  step  calculation  is  performed  in  Links  2,  3.  and  4  as  outlined  in 
lief.  2,  viih  modifications  in  the  solution  procedure  to  account  for  the 
changes  given  above  in  Iq,  20  through  24  and  27  through  32. 

In  the  solution  of  the  energy  equation,  a  separate  function  subprogram  was 
used  to  calculate  the  difference  analog  of  the  mixed  temperature  derivative, 
thus  providing  the  capability  of  testing  and  evaluating  alternative  formula¬ 
tions  to  that  given  by  Eq.  27  through  30.  The  only  alternative  subprogram 
written  thus  far  is  one  that  simply  sets  the  mixed  difference  identically 
to  zero.  This  was  used,  as  will  be  seen  in  the  checkout  section  below, 
to  help  asses?  the  significance  of  the  mixed  derivative  term. 


Checkout  of  the  Anisotropic  Capability 

Checkout  of  the  anisotropic  program  capability  was  performed  by  comparison 
to  previously  obtained  results  with  carbon  cloth-phenolic  treated  as  iso¬ 
tropic  (the  same  control  case  used  to  check  out  Task  l.A  and  discussed  in 
the  Review  section,  see  Fig.  3).  Assumed  values  were  used  for  the  aniso¬ 
tropic  input  data.  Because  checkout  was  not  performed  with  respect  to 
measured  results  but  only  by  comparison  to  computed  isotropic  results, 
designation  of  checkout  as  successful  or  unsuccessful  was  a  qualitative 
decision  based  on  the  predictability  of  the  direction  rather  than  the 
degree  of  deviation  of  the  anisotropic  results  from  the  isotropic.  Thus, 
the  objective  of  the  checkout  was  to  see  whether  the  anisotropic  results 
deviated  from  the  isotropic  in  a  predictable  manner,  over  a  representative 
tango  of  orientations,  using  assumed  anisotropic  conductivities  as  input 
which  differed  only  slightly  from  th~  isotropic  values. 

Cases  evaluated  included  orientations  of  0,  45,  90,  and  135  degrees,  in 
which  major  and  minor  conductivities  (both  temperature  dependent)  differed 
from  the  isotropic  conductivity  by  a  constant  (Fig.10).  An  orientation 
of  45  degrees  only  was  employed  for  another  case  in  which  the  difference 


dc  Conductivity  for  Checkout 


between  conductivities  varied  with  temperature  by  permitting  the  aniso¬ 
tropic  conductivities  to  vary  linearly  with  the  known  isotropic  conduc¬ 
tivity,  as  follows: 

K  (T)  =  1.1  K.  (T),  MT)  =  0*9  K.  (T). 

All  cases  were  first  run  using  the  subprogram  which  provides  a  vanishing 
mixed  difference  in  the  energy  equation  and  then  with  the  second  order 
analog  expressed  by  Eq.  2?  through  39.  Ia  this  way,  variouo  anisotropic 
effects  could  be  isolated  and  evaluated  separately. 

The  results  of  the  checkout  arp  given  in  Fig.  11  in  the  form  cf  a  comparison 
of  surface  and  interior  temperatures  obtained  for  the  various  orientations 
(cases  a  through  e)  after  12  seconds  of  firing*  The  classification  of 
cases  a  through  e  (Fig.  11  )  will  be  maintained  throughout  the  remainder 
of  the  discussion.  Temperatures  given  in  parentheses  (Fig.  11  )  are  those  -  - 
obtained  with  a  vanishing  mixed  temperature  difference  in  the  energy  equa¬ 
tion.  Cases  a,  c,  and  d  are  analyzed  below  in  detail,  and  Cdsc  e  in 
somewhat  less  detail.  Case  b  is  not  discussed  because  the  reasoning  employed 
was  similar  to  that  of  Case  a. 


Case  a,  45-degree  orientation  (0  «  45  degrees*),  K.  =  K„  =  K_.  (T) , 

"■ 1  M  x  y  iso 

Kxy=-0.2  (10"°),  1^.  =  Q.  At  the  inside  surface,  the  normal  conductivity 

is  obtained  from  Eq.  23  as  follows: 


*In  general,  to  calculate  the  equivalent  value  of  0  for  a  given  orientation 
a,  we  can  use  0  «  90  -  a  degrees  for  0  &&  s90  degrees;  6  *>  2?0  -  Ct  degrees 
for  90  <  ‘X  <  180  degrees 


Sin* 


Therefore,  to  the  left  of  the  throat,  ve  heve  K  <  K.  because  &f/Bx  <  0, 

D  10  0 

and  thus  we  should  have  T  .  >  T.  at  the  surface.  This  inequality  is 

aniao  iso  *  * 

satisfied  in  Fig.  11.  the  calculated  values  being  T  .  =  5061  E  and 

aniso 

T.  =  5053  R-  To  the  right  of  the  throat,  the  reverse  tendency  is  reflected 
(Fig.  ll);  Tan.so  is  16  R  less  than  T.gQ. 

It  might  be  argued  that  a  10  or  16  degree  difference  is  smaller  than  the 
expected  error  in  the  temperature  and  therefore  should  be  disregarded  aa 
not  being  significant.  This  reasoning  vould  be  false  because  ve  are  com¬ 
paring  two  calculated  temperatures,  not  a  calculated  temperature  with  a 
measured  temperature.  There  are  several  sources  of  error  in  this  sort  of 
calculation.  One  error  is  caused  by  the  assumption  of  a  nonrealistic 
model  and  is  a  discrepancy  between  real  temperatures  and  idealized  tem¬ 
peratures.  Another  error  is  caused  by  the  method  of  solution  and  is  a 
measure  of  how  closely  the  numerical  solution  of  the  finite  difference 
equations  approximates  the  exact  solution  of  the  continuous  model  *  3. 

Both  of  these  errors  can  he  quite  large.  However,  neither  applies  to  uhe 
present  discussion  because  we  are  comparing  two  values  which  are  calculated 
numerically  from  almost  identical  input  data.  Thus,  as  long  as  the  numerical 
procedures  are  reasonably  accurate  and  stable,  the  difference  between  the 
two  calculated  values  need  not  be  large  to  be  significant.  As  will  be 
shown,  the  consistency  with  which  the  expected  trends  are  reflected  in 
the  calculated  results  indicates  that  these  differences,  though  small-  are 
indeed  significant. 


To  evaluate  the  interior  temperatures  attained,  two  subcases  must  be 
analyzed.  First,  in  the  calculation  made  with  a  vanishing  mixed  tempera- 
ture  difference  (i.e.,  vith  =  0)  we  see  that,  because  K'^  is  zero 

as  well  and  =  K^g{),  the  anisotropic  energy  equation,  reduces  identi¬ 

cally  to  the  isotropic  case.  Thus,  any  interior  anisotropic  effect  felt 


auat  be  only  a  reflection  of  the  surface  effect.*  litis  is  borne  out  in  the 

temperatures  calculated,  where  ^an^80  =  1529°  >  1522^  =  T^so  ^faa  left 

of  the  throat  and  T  .  =  1289°  <  1293°  =  T.  to  the  right  of  the  throat, 

aniso  iso 

In  the  second  subcase,  that  calculated  with  ^  I  58  expressed  by  Eq.  27 
through  30,  we  see  that  the  only  difference  in  the  energy  equation  is 
caused  by  the  presence  of  the  mixed  derivative.  Thus,  we  can  write 


_  dr 

p°  a? 


aniso 


_  st 

7  P0  If 


-0.4  (10-5) 


Uso 


dxdy 


(33) 


Here  we  have  isolated  the  mixed  derivative  term  of  the  energy  equation  so,- 
that  the  present  subcase  serves  as  a  partial  check  of  the  aignif icancjt?  of 
the  second— order  mixed  difference  analog  by  comparison  of  results  w^tb, 
the  subcase  in  which  the  vanishing  mixed  difference  was  employed. 


Because  of  the  gradients  caused  by  surface  heating,  0  T/dxdy  will  in  gen¬ 
eral  be  positive  to  the  left  of  the  throat  and  negative  to  the  right  of 
the  throat.  This  is  caused  by  generally  decreasingly  negative  dT/dx  and 
dr/dy  with  increasing  y  and.  x,  respectively,  to  the  left  of  the  throat 
and  by  increasingly  positive  dr/dx  and  negative  dr/dy  with  increasing  y 
and  x,  respectively,  to  the  right  of  the  throat.  Therefore,  we  see  from 

Eq .  33  that  T  ,  should  tend  to  increase  with  time  more  slowly  than 
^  aniso 

T-  to  the  left  of  the  throat  and  more  rapidly  to  the  right.  Thus, 

X  SO 

except  near  the  surface  where  the  reverse  effect  dominates  as  discussed 

above,  we  should  have  T  <  T.  to  the  left  of  the  throat  and 
^  aniso  iso 

T  ,  >  T.  to  the  right.  Again,  we  see  that  these  tendencies  are 

aniso  iso  b  ^ ' 

reflected  in  the  calculated  interior  temperatures  (Tan£so  *  1508°  <  1522°  » 

T,  to  the  left  of  the  throat  and  T  .  =  1321°  >  1293°  **  T,  to  the 

iso  aniso  iso 

right).  Thus,  the  interior  anisotropic  effect  is  the  reverse  of  the  sur¬ 
face  effect  and  should  be  accounted  for  by  use  of  a  nonvanishing  mixed 


*This  subcase  illustrates  the  value  of  simulating  the  boundary  curves  y  = 
f(s,r)  accurately  in  the  analysis  because  the  only  anisotropic  effect  in  the 
solution  is  felt  in  the  value  of  Kjj,  which  depends  directly  on  the  slope 
di/dx.  To  achieve  a  similar  effect  with  alternately  horizontal  and  vertical 
segments  using  Kx  and  Ky  respectively  would  require  a  significantly  fiDer 
mesh  and  thus  considerably  more  machine  time  for  the  solution.  This  sub¬ 
case  further  sorves  as  a  partial  checkout  of  the  significance  of  the 
isolated  surface  anisotropic  effect. 


kk 


difference  analog.  We  note  in  Fig.  11  that  the  anisotropic  surface  tem¬ 
peratures  calculated  with  the  vanishing  sized  difference  are  not  influ¬ 
enced  by  the  reverse  interior  effect  md  are  therefore  higher  to  the  left 
of  tie  throat  and  lover  to  the  right  than  those  calculated  with  the  second 
order  difference  analog. 

Even  at  the  throat  plane,  both  on  the  surface  and  in  the  interior,  the 
results  obtained  are  in  agreement  with  expected  effects,  i.e. ,  the  calcu¬ 
lated  values  exhibit  slight  trends  characteristic  of  the  corresponding 
values  obtained  to  the  right  of  the  throat  because  of  the  steeper  surface 
slope  to  the  right  of  the  throat  (20  degrees  as  compared  to  10  degrees). 

Case  Co  90-degree  orientation,  =  ILgo  -  0.2  (l0-^) , 

K.  +  0.2  (10“5),  K  =  K’  =  0.  At  the  surface,  we  get 
iso  '  ' 1  xy  xy  >  b 


Therefore,  we  expect  T  .  <  T ,  at  all  points  along  the  inBide  surface, 

including  the  throat,  because  jof/9x|  <  1  everywhere  on  the  surface  (the 
maximum  inclination  from  the  horizontal  being  20  degrees).  Indeed,  calcu¬ 
lated  results  agree  with  the  expected  direction  of  deviation  (Fig.  ll). 

For  the  interior  temperatures  calculated,  there  is  only  one  case  to  analyze 

because  K  6^  T  =  0  no  matter  how  6^  T  is  calculated.  Actually  two 

cases  were  run,  one  with  vanishing  $  T  and  the  other  with  the  second 

xy 

order  calculation  given  by  Eq.  27  through  30,  to  check  the  programming  of 
the  latter.  The  results  obtained  were  identical  for  the  two  cases,  as 
expected. 
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With  the  normally  directed  heat  flux  at  the  inside  surface  applied  con¬ 
siderably  more  in  a  radial  than  in  an  axial  direction  (the  surface  being 
inclined  only  20  degrees  at  most  frcaa  the  horizontal),  the  radial  tem¬ 
perature  gradients  should  be  much  steeper  in  the  interior  than  the  axial 
and  should  also  vary  more  abruptly.  Thus,  ve  can  expect  that 


h~T  1  ST  .  S2T 


Then,  because  K  =  K.  -  0.2(l0“5),  K  =  K.  +  0.2(l0~5),  =r  =K! 

x  iso  v  '  ’  y  xso  v  7  ’  Tc  y  x 

and  K  =H£y  -0,  we  would  have 


iso 


pc§|  ~pc  HI  +  0.2(10-5)  ( 

aniso  iso 


S^T  _1  ST  S2T 

ay2  +  y  3y  ' 


3x‘ 


ISO 


and  we  would  expect  T  .  >  T.  at  positions  which  are  far  enough  into 

r  anxso  iso  v  ^ 

the  interior  to  escape  domination  by  the  reverse  anisotropic  effect  at 
the  surface.  This  is  reflected  in  the  calculated  results  (Fig.  11  ). 


Case  d.  0-degree  orientation,  K  =  K_  =  IL  +  0.2(l0~')),  K  =  Kc  «  K.  - 
- -  x  y  iso  '  y  4  iso 

0.2(10^  ),  IC  =  K'  =0.  The  reasoning  here  is  just  the  reverse  of  that 
used  for  Case  c.  Thus,  ve  get 


VKi8o^.2(l°-5)  [l  ^)2]  /  [l  +  (I)2]  <j  > 


iso 


180 


if 

mi 

£  1 

if 

mi 

>  i 

and,  therefore  we  expect  T  ,  >  T.  at  the  surface,  as  shown  in  Fig. 11. 

r  anxso  iso 


In  the  interior,  we  get 


*11  .  ‘*11.  (0 ♦*$-§)<  *$| 

aniso  iso  Sy  J  * 


iso 


3x 


iso 
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and  ve  expect  T  .  <  T.  in  the  interior,  except  near  the  surface  where 

aniso  iso  r 

the  reverse  effect  dominates. 


Case  e. 


-0.  IK! 


iso 

interior, 


45**'degree  orientation,  K  =  K  =  K.  ,  K  =  -0.1K.  , 

&  ’  x  y  iso’  xy  iso 

Similar  reas/.'iug  applies  at  the  surface  as  in  Case  a. 
however,  ve  must,  for  the  first  time,  assess  the  effect 


2 

as  well.  In  particular,  by  analyzing  the  subcase  for  which  6^  T 
to  be  identically  0,  we  isolate  the  effect  of  the  term  2K^.  3t/3x 


in  the  energy  equation  as  follows: 


K'  = 
xy 

In  the 

of  K' 
xy 

is  taken 

3T/3y 


dTs 

*371 


anise 


~  PC 


iso 


+  2K' 

xy 


ST  3T 

• 


Therefore,  because  cVf/dx,  3T/Sy,  and  are  all  negative  to  the  left  of 

the  throat  (K1  =  -0.1K!  <  0  because  K!  >  0 — see  Fig.  10) ,  we  have 

'  xy  iso  iso  6  '  * 


aniso 


and  thus  T  .  <  T.  in  the  interior  to  the  left  of  the  throat.  This 

aniso  iso 

is  borne  out  by  the  calculated  results  (T  .  =  1517°  <  1522°  =  T.  ). 

'  aniso  iso' 

Thus,  we  see  that  the  presence  of  the  term  2K1  dl'/3x  dT/3y  in  the  energy 

equation  causes  a  reversal  of  the  surface  effect  (in  contrast  to  Case  a) 

and  thus  constitutes  a  significant  anisotropic  effect.  To  the  right  of 

the  throat  we  obtain  the  consistent  result  that  in  the  interior  T  .  is 

aniso 

higher  than  T. 

6  iso 


In  the  second  subcase  of  Case  e,  in  which  the  second  order  mixed  difference 
is  employed,  the  direction  of  deviation  from  the  isotropic  results  in  the 
interior  is  the  same  as  in  Case  a  but,  as  shown  in  Fig.  11,  is  more  ex¬ 
treme  because  of  reinforcement  by  the  presence  of  the  nonzero  term 
2K*  ol/3x  dT/dy  in  the  energy  equation. 
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TASK  l.C,  MORE  THAN  ONE  CHARRING  MATERIAL 


The  major  effort  required  for  extending  the  2D-ABLATE  program  to  treat 
the  case  of  multiple  charring  ablative  thrust  chambers  was  in  programming 
rather  than  in  modification  of  the  model  equations  or  of  the  discretiza¬ 
tion  procedures.  The  latter  was  reflected  only  in  the  treatment  of  the 
generated  gas  mass  flux  across  a  material  interface;  i.e.,  in  the  simula¬ 
tion  of  interface  conditions  for  solution  of  the  continuity  equation  in 
adjacent  charring  materials.  The  method  adopted  was  to  continue  to  dis¬ 
cretize  the  continuity  equation  at  the  regular  points  only  (as  in  the 
"old  '  version  of  the  program  in  which  only  one  charring  material  could  be 
treated)  as  if  no  intervening  interface  existed.  However,  t>e  thermochemical 
properties  and  charring  temperatures  associated  with  each  ~eint  must  be 
based  on  the  properties  assigned  to  the  material  in  which  the  point  lies. 

The  bulk  of  the  programming  required  for  the  program  extension  was  thus 
of  a  logical  or  accounting  nature.  Additional  input  allocation  and  assign- 
’.ent  of  additionally  indexed  arrays  were  required  for  the  constants  and 
temperature-dependent  functions  used  to  describe  the  gas  generation  and 
ei'  ing  reactions  taking  place  within  several  rather  than  a  single  charring 

mri.'  ial.  In  Link  1,  the  selection  of  the  bounding  mesh  lines  for  the 
solution  of  the  continuity  equation  was  extended  from  a  region  lying  within 
a  singx-  material  to  a  multimeter ic i  range.  In  Links  3  and  4,  accounting 
changes  .-re  required  to  keep  track  of  the  regions  and  material  properties 
as  the  me  h  lines  cross  the  interfaces  and  to  interrelate  generated  gas 
mass  fluxes  from  one  region  to  another.  Checkout  of  the  added  Task  l.C 
program  capability  was  performed  by  comparison  with  previously  computed 
results  for  a  carbon  cloth-phenolic  thrust  chamber  backed  by  a  stainless 
steel  shell  (the  same  control  case  used  to  check  out  Tasks  l.A  and  ],B). 

Two  cases  were  run.  In  both,  cases  the  carbon  cloth  material  was  axially 
divided  into  three  separate  side-by-aide  charring  materials  (Tig. 12  ),  each 
requiring  separate  input  of  the  same  material  properties.  In  the  first  case, 
all  materials  were  treated  as  isotropic  so  that  essentially  the  control 
case  was  rerun  to  verify  that  all  the  charring  calculations  were  being 
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0  DEGREES 


peri  aed  properl}  in  each  region.  ’rs'°  .:•  lyrical  results  obtained  were 
clow.e  but  not  the  same  as  those  obtained  in  the  control  case.  The  dis¬ 
crepancy  was  caused  by  truncation  error  incurred  in  discretisation  of  the 
audit ional  interface  conditions  required  to  express  the  continuity  of  heat 
flux  from  material  to  material:  i.e.,  because  of  the  highly  nonlinear  char¬ 
acter  of  the  variation  of  hear  iiux  with  axial  distance  occurring  at  times 
in  the  neighborhood  of  the  interfaces.  This  was  particularly  true  is  the 
early  stager  of  the  run  '.rhen,  in  the  vicinity  of  an  interface,  the  tem¬ 
perature  had  risen  steeply  on  one  side  and  not  at  all  on  the  other  side. 
Despite  this  difference  in  results,  the  checkout  was  considered  to  be  suc¬ 
cessful  because  the  charring  mechanism  worked  as  intended  in  the  several 
regions. 

In  the  second  checkout  case,  the  configuration  and  materials  were  the  same 
as  in  the  first  case  but  all  three  charring  materials  were  treated  as 
anisotropic,  the  chamber  and  exit  materials  being  given  a  0-degree  orienta¬ 
tion  (0  -<=  90  degrees)  and  the  throat  materia’  a  90-degrs;e  orientation 
(0=0  degreed).  As  expected,  the  surface  temperatures  obtained  in  the 

,  0  degree  materials  were  higher  and  the  interior  temperatures  were  lower 

( 

tear  the  corresponding  temperatures  obtained  in  the  first  case,  and  the 

j  reverse  was  true  in  the  throat  material. 

1 

1 
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TASK  2,  EFFECTIVE  ABLATION  PROPER! it3 


INTRODUCTION 

The  pr^maiy  objective  of  Task  2  wss  to  define  the  effj  :tive  thermophysieal 
and  uhermochemical  properties  of  reinforced  phenolic  resin  systems  neejec 
for  ablative  design  with  the  2li-ABLATF  program.  The  two  systems  of  in¬ 
terest  in  the  program  were  phes.olic-ref raeil  and  phenol ic/carbon  cloth. 
Properties  for  these  systems  ere  used  for  a  portion  of  the  2-D  compute*, 
runs  of  Task  3* 

Task  2  also  fulfilled  the  follovJng  secondary  objectives:  (l)  operational 
checkout  of  the  2D-ABLATE  computer  program,  (2)  investigation  of  trial 
and  error  techniques  for  matching  experimental  data  with  computed  results 
and  (3)  analysis  of  parame+ric  effects  in  areas  which  overlap? 2d  with 
Task  3, 

The  physical  properties  required  in  the  2-D  program  either  as  functions 
of  te^perat— e  (l,  2,  4,  5)  or  as  constants  (3,  6,  7)  are  as  follows: 

] ,  Density-specific  heat  product 

2.  Fraction  of  resin  pyrolyzed 

3.  Heat  of  pyrolysis  of  the  resin 

4.  Enthalpy  oi  the  pyrolysis  gases 

5.  Thermal  conductivity  of  the  char  and  virgin  material  parallel 
and  perpendicular  tc  the  reinforcement 

6.  Char- reinforcement  reaction  constants 

7.  Erosion  constants 

Insofar  as  was  practical,  these  properties  were  obtained  from  the  literature. 
The  values  for  thermal  conductivity  and  the  enthalpy- temperature  relationship 
for  the  pyrolysis  gases  were  to  be  obtained  by  matching  the  calculated  temp¬ 
erature  and  char  profiles  of  2D-ABLAIE  for  assumed  property  values  to  the 
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experimental  temperature  profile  anti  depths  from  ablative-vailed  rocket 
motor  firings. 

Appropriate  data  for  the  effective  property  determinations  were  v=ry  lia¬ 
ised  because  of  the  three  basic  requirements  such  data  had  to  meet:  (1) 
relatively  constant  chamber  pressure  and  mixture  ratio,  (2)  a  sufficient 
number  of  axial  end  radial  temperature  measurements  to  define  the  tempera¬ 
ture  profile  as  a  function  of  time,  and  (3)  injector  performance  which 
gave  minimal  and  circumferential ly-uniform  surface  erosion. 

The  test  data  ultimately  selected  for  determination  of  the  properties  of 
phenolic/carbon  cloth  were  generated  by  the  Research  Division  of  Rocket- 
dyne  under  Contract  NAS7-304.  The  results  of  four  300-second  firings  with 
the  propellant  combination  OFg/HMH  at  a  nominal  chamber  pressure  of  ! 10 
psia  were  employed.  These  tests  provided  data  at  three  reinforcement 
orientations*,  0,  60,  and  84  degrees,  together  with  the  thermal  response 
of  a  nonablative  ATJ  graphite-walled  chamber. 

No  entirely  satisfactory  rocket  firing  data  were  found  to  determine  the 
properties  of  phenolic- lefrsai 1  ablatives.  Under  NASA  Contract  NAS9-150 
(Ref.  7),  effective  properties  of  phenolic-ref rasil  had  been  estimated 
by  means  of  the  previous  verison  of  2D-ABLATE  (2D-CHAR).  However, 
these  estimates  were  based  upon:  (l)  the  assumption  of  a  single  charring 

materia],  and  isotropic  properties,  and  (2)  upon  back  wall  temperatures 
which  showed  highly  damped  response. 

The  provious  estimates  of  phenolic-refrasil  properties  from  Ref.  7  were 
rechecked  by  calculating  the  temperature  history  of  the  Apollo  engine  with 
2D-ABLAIE  and  comparing  the  results  to  the  old  calculations  and  to  the 
original  experimental  data. 

In  the  following  sections,  the  recommended  properties  for  phenolic/carbon 
cloth  and  phenolic-refrasil  are  presented  graphically  and  in  a  numerical 
form  directly  applicable  for  input  to  ths  2-D  program.  Applicable  proper¬ 
ties  i or  the  hard  throat  insert  materials,  ATJ  graphite  and  silicon  carbide, 

^Orientation  is  used  here  and  throughout  the  Task  2  and  3  discussions  to  mean 
the  counterclockwise  rotation  6  as  defined  on  p.31  for  use  in  2D-ABLAIE  rather 
than  the  usual  clockwise  orientation  angle.  See  footnotes  on  pp.31  and  41 
for  the  precise  relationship  between  the  two. 
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are  also  presented.  Subsequently,  the  experimental  data  and  catching 
procedures  employed  to  define  thermal  conductivity  for  the  two  ablative 
systems  are  described. 

The  analysis  of  the  experimental  data  used  in  the  phenolic/carbon  cloth 
property  determinations  is  particularly  interesting  because  it  clearly 
shows  the  artifical  nature  of  the  usual  boundary  conditions  assumed  in 
regenerative ly  cooled  rocket  chamber  heat  transit  vaalyses.  It  also 
shows  the  significance  of  injector  effects  on  heat  transfer. 


PROPERTIES 

Product  of  Density  and  Specific  Feat 

The  recommended  values  for  the  0  product  of  phenolic-ref rasil,  phenolic/ 
carbon  cloth,  ATJ  graphite  and  silicon  carbide  are  given  in  Fig.  13.  The 
corresponding  numerical  input  for  2D-ABLATE  usage  is  presented  in  Table  2. 
The  method  of  computation  and  source  data  for  phenolic  refrasil  were  ob¬ 
tained  from  Ref.  7.  A  similar  procedure  using  manufacturer’s  listed 
values  for  virgin  density  and  the  same  degree  of  resin  pyrolysis  was  em¬ 
ployed  for  phenolic/carbon  cloth.  The  values  for  ATJ  graphite  and  silicon 
carbide  were  obtained  from  Ref.  12  and  13. 


Fraction  Resin  Pyrolyzed 

The  relationship  between  degree  of  resin  pyrolysis  and  temperature  is 
given  in  Fig.  14.  It  represents  the  recommended  input  values  for  the 
2D-ABLATE  program  for  ooth  phenolic  refrasil  and  phenolic/carbon  cloth. 

In  terms  of  the  program  input  variables 

Maximum  fraction  .esinw  pyrolyzed  «  0.45 
Maximum  temperature  of  pyrolysis  range  *»  1750  R 
Minimum  temperature  of  the  pyrolysis  range  <=  750  R 
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Figure  13.  Density  -  Specific  Heat  Product  Fir 
Ablative  Chamber  Materials 


TEMPERATURE,  R 

Figure  14  ,  Fraction  Resin  Pyrolyzed 
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TABLE  2 


RECOMMENDED  DENSITY-SPECIFIC  HEAT  PRODUCT  INPUT  TO 
2D-ABLATE  PROGRAM  FOR  ABLATIVE  WALL  AND 
HARD  THROAT  MATERIALS  / 


Phenolic/Carbon  Cloth 


Phenolic  Refrasil 


ATJ  Graphite 


Silicon  Carbide 


Temperature, 

R 

pc 

Btu/in.^-E 

350 

0.015 

360 

0.015 

1060 

0.017 

1460 

0.018 

2060 

0.0]8 

3460 

0.019 

7460 

0.019 

360 

0.0093 

660 

0.015 

900 

0.018 

1460 

0.015 

I960 

0.016 

3460 

0.018 

5460 

0.020 

570 

0.015 

560 

0:0175 

1460 

0.0195 

2060 

0.021 

3460 

0.022 

7460 

0.022 

300 

0.012 

600 

0.023 

1000 

0.029 

300G 

0.037 

6000 

0.045 

The  curve  of  Fig.  14  is  based  upon  the  data  of  Ladaeki  (Ref.  14)  for  con¬ 
stant  temperature  or  constant  heating  rate  pyrolysis  in  a  furnace.  The 
data  indicate  that  the  fraction  pyrolyzed  at  a  given  temperature  is 
practically  independent  oi  heating  rate  over  the  range  of  heating  rates 
tested. 

In  initial  computer  runs  conducted  with  a  TAP-IV,  1-dimensional  computer 
program,  an  upper  temperature  limit  for  pyrolysis  of  1600  F  (2060  R)  was 
assumed.  The  subsequent  reduction  to  1750  &  appeared  to  improve  the  fit 
of  experimental  temperature  histories  to  computer  calculations. 


Heat  of  Pyrolysis  and  Pyrolysis-Gas  Enthalpy 

I 

The  heat  of  pyrolysis  for  phenolic  resins  has  been  shown  to  vary  with 
temperatures  (Ref.  14).  However,  the  2D-ABLATE  input  format  requires 
the  assumption  of  a  constant  value.  A  standard  value  for  heat  of  pyrolysis 
of  450  Btu/lb  was  employed  in  ell  computer  runs  of  Tasks  2  and  3  in  this 
program.  Because  the  heat  of  pyrolysis  represents  only  a  small  fraction 
of  the  total  heat  absorbed  by  the  resin  and  subsequent  pyrolysis  gas  re¬ 
actions,  the  choice  of  the  value  employed  for  heat  of  pyrolysis  produces 
negligible  changes  in  calculated  char  rates  over  a  range  of  several  hun¬ 
dred  Btu/lb. 

The  2D-ABLATE  program  assumes  that  the  degree  of  pyrolysis  and  the 
composition  of  the  pyrolysis  gases  are  thermodynamically  rather  than 
kinetically  controlled  and,  that  the  pyrolysis  gases  are  in  local 
thermal  equilibrium  with  the  char.  As  a  result,  the  enthalpy  of  the 
pyrolysis  gas  from  phenolic  resin  is  defined  as  a  function  of  tempera¬ 
ture  only.  During  the  Apollo  RCS  nozzle  qualification  program  (Ref.  7), 
this  functional  dependence  was  deduced  by  two  parallel  efforts.  Develop¬ 
mental  test  data  for  the  Apollo  RCS  engine  were  matched  to  computer 
predictions  of  cuter  skin  temperatures  and  char  depth  for  aetumed  enthalpy- 
temperature  relations.  At  the  same  time,  Ladaeki 's  experime:  its  (Ref.  14) 
obtained  data  on  the  chemical  composition  of  the  pyrolysis  vapors.  An 
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enthalpy-temperature  curve  was  obtained  from  the  respective  mole  fractions , 
heats  of  formation,  and  sensible  specific  heats  of  the  gaseous  species 
observed  in  Ladacki's  experiments  (together  with  the  assumed  conversion  of 
these  species  to  carbon  monoxide  and  iydrogan  at  higher  temperature) ,  which 
substantiated  the  curve  obtained  by  matching  the  ACS  engine  firing  data.  The 
enthalpy'-tempex’atnre  curve  obtained  in  the  Apollo  qualification  stndies  is 
shown  in  Pig,  15.  Heccmmsnded  input  data  to  the  2D-ABLATE  program  are  listed 
in  Table  3*  The  higher-temperature  portion  of  the  curve,  which  is  based  entire¬ 
ly  upon  the  assumption  that  the  pyrolysis  gas  is  composed  of  carbon  monoxide  gas 
and  hydrogen  in  this  temperature  range,  is  uncertain  above  approximately  3500  R. 

TABLE  3 

EFFECTIVE  ENTHALPY  OF  PHENOLIC  RESIN  PYROLYSIS  VAPORS 


Temperature , 

R 

Enthalpy 

Btu/lb 

0 

0 

700 

0 

1110 

750 

1560 

1850 

2060 

3850 

7460 

9250 

Thermal  Conductivity 

The  values  for  thermal  conductivity  of  phenolic-refrasil  parallel-to  and 
uormal-to  the  silica  reimorcemeat  is  shown  in  Fig.  16.  The  values  shown 
were  determined  in  the  Apollo  qualification  program  (Ref.  15)  by  u  guarded 
hot-plate  method  using  samples  preconditioned  at  the  temperature. 

Experimental  data  were  obtained  at  temperatures  up  to  15 10  F,  and  the 
values  defined  for  higher  temperatures  were  estimated  theoretically  by 
assuming  internal  radiant  transport  contributions  (Ref.  7)*  The  thermal 
conductivity  values  obtained  from  the  hot  plate  measurements  were  employed 
with  a  variable  pyrolysis  gas  enthalpy  to  natch  the  outer  wall  temperatures 
of  the  SE-8  engine  during  actual  firings  to  the  computed  results  of  2D-CHAR. 
Thus,  the  thermal  conductivity  values  for  phenolic-refrasil  given  iu  Fig.  16 
and  Table  4  are  consistent  with  the  enthalpy-temperature  relations  of 
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Fig.  15  and  Table  3*  The  data  matching  done  with  2D-CESB  was  checked 
with  2D-ABLAIE  as  described  in  a  later  section  of  this  report. 

For  phenolic/carbon  cloth,  the  results  of  the  temperature-matching  proced¬ 
ures  between  the  OF^/MMH  model  motor  firing  data  and  the  calculations  of 
2D-ABLATE  predict  an  essentially  constant  value  for  thermal  conductivity 
of  8.0  x  10  5  Btu/ia.-sec-R  parallel  to  the  carbon  cloth  reinforcement 
(major  direction).  A  slightly  better  data  fit  was  obtained  by  employing 
a  moderate  temperature  dependency,  with  the  thermal  conductivity  decreas-- 
ing  from  10  x  10  at  ambient  temperatures  to  6  x  10  at  1000  F  remain¬ 
ing  constant  for  higher  temperatures.  Because  of  the  r elatively  low 
adiabatic  wall  temperatures  in  the  OFg/HMH  model  motor  (as  discussed  in 
a  later  section),  the  same  char  depth  was  obtained  with  either  defini¬ 
tion  of  thermal  conductivity.  However,  the  former  value  of  thermal  con¬ 
ductivity  predicts  a  more  conservative  char  depth  for  higher  adiabatic 
wall  temperatures  and  therefore  is  recommended  until  high  adiabatic  wall 
temperature  data  become  available. 

The  recommended  value  for  thermal  conductivity  of  phenol ic/ carbon  cloth 
normal  to  the  reinforcement  (minor  direction)  is  7*0  x  10  Btn/in.-sec-B 
between  ambient  temperatures  and  2500  B  increasing  slowly  to  1.5  x  10 
Btn/in.-sec-B  at  7000  E.  Because  of  general  similarity  between  the  con¬ 
ductivities  of  both  phenolic-refrasil  and  phenolic/carbon  cloth  normal 
to  the  reinforcement,  a  similar  temperature  dependence  is  predicted 
despite  the  lack  of  high-temperature  experimental  date  for  comparison. 

The  postulated  temperature  dependence  gives  conservative  predictions 
(greater  char  depths^  and  is  therefore  recommended.  Recommended  values 
for  the  thermal  conductivity  of  phenolic/carbon  cloth  are  given  in  Table  4. 

Beconsaended  thermal  conductivity  input  data  for  ATJ  graphite  and  silicon 
carbide  obtained  from  Ref.  12  are  also  presented  iu  Table  4. 


Char  Beinforcement  Reaction 


Provision  is  made  in  21>- ABLATE  for  reactions  between  char  and  reinforce¬ 
ment.  No  such  reactions  need  be  considered  fer  phenolic/earboa  cloth 
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char  systems.  However,  the  highly  endothermic  con"ersion  of  silica  and 
carbon  to  silicon  carbide  was  consi&ei-d  for  phenolic-refrasil  ablators. 

In  fact,  chemical  analysis  of  clJir  specimens  froa  the  SE-8  motor  walls 
(Ref.  7)  indicated  no  formation  of  silicon  carbide.  Pending  experimental 
verification  of  appreciable  carbide  formation,  no  char-reinforcement  re¬ 
action  inputs  are  recommended  for  phenolic-ref rasil. 

Erosion  Rate  Constants 

The  melting  point  of  carbon  is  above  7000  F.  Therefore,  the  logical  modes 
for  surface  erosion  of  either  ATJ  graphite  or  carbon-cloth  reinforced 
phenolic  char  are  vaporization  and  chemical  reaction  with  corrosive  com¬ 
bustion  gas  species.  Thermodynamic,  chemical  equilibrium  calculations 
and  experimental  data  (Ref.  16)  indicate  that  carbon  is  compatible  with 
HF  and  CO  but  is  severely  attacked  byHgO  and  COg  and,  to  -a  lesser  extent, 
by  Hg.  The  recommended  erosion  input  data  for  either  phenolic  carbon 
cloth  or  ATJ  graphite  are  listed  in  Table  5»  They  include  vaporization 
constants  from  Ref.  17  and  chemical  reaction  rate  constants  evaluated 
from  data  in  Ref.  18,  19,  and  20  for  reactions  with  HgO,  COg,  and  Hg. 

With  phenolic-ref rasil,  melting  of  the  silica  reinforcement  is  the  primary 
erosion  mechanism;  the  secondary  mechanism  being  the  char  vaporization. 
Recommended  erosion  input  data  for  phenolic-refrasil  are  listed  in  Table  6. 

For  silicon  carbide,  throat  erosion  occurs  primarily  by  means  of  decompo¬ 
sition  at  temperatures  above  4500  F.  This  erosion  is  treated  as  a  vapor¬ 
ization  mechanism  with  a  total  vapor  pressure  equal  to  the  sum  of  the 
partial  pressures  of  the  decomposition  products.  Recommended  erosion 
parameters  based  upon  vapor  pressure  data  in  Ref.  21  are  listed  in  Table  7. 
The  proposed  constants  have  not  been  verified  by  comparison  to  rocket 
motor  firing  data. 


Parameter 


Value 


Vaporization  Mechanism 

Diffusion  Coefficient,  in.'/sec 

Vapor  Ho 3 scalar  Weight 

Vapor  Pressure  Constant,  psia 

Molar  Latent  Heat,  Btiy 1 j  mole 

Mass  Latent  Heat,  Btu/lb 

Density,  lb/in.3 

Specific  Heat  of  Vapor,  Biu/lb-F 

HgO  Reaction  Constants 


Data  Input 
location 


0.04 

896  +  8n* 

12 

897  +  8n 

5  x  108 

898  +  8n 

514000 

899  J->  8n 

25000 

900  +  8n 

0.081 

901  +  8n 

0.44 

902  +  8n 

Molecular  Weight  of  ELO 

Free  Stream  Mole  Fraction  of  HgG 


18 

From  performance 
calculation 


Hate  Constant,  in. /sec  j 

Activation  Energy,  Btu/lb  mole  ! 

Heat  o?  Reaction,  Btu/lb  j 

Mass  Diffusion  Coefficient^  in./sec 
Product  Specific  Heat,  Btu/lb-B  j 


2,4* 

4.21  x  10* 
5000 
0.034 
0.47 


Hg  Reaction  Constants 


From  performance 
calculation 

0.095  , 

6.8  x  104 
4500 

Mass  Diffusion  Coefficient,  in. ‘/sec  0.095 
Product  Specific  Heat,  Btu/lb-R  J  1.0 


Molecular  Weight  of 

Free  Stream  Mole  Fraction  of  Eg 

Rate  Constant,  in./aec 
Activation  Energy,  Btu/lb  mole 
Heat  of  Reaction,  Btu/lb  , 


936  +  8n* 

937  +  8b 

938  +  8n 

939  +  8n 

940  +  8n 

941  +  8n 

942  +  8n 


976  +  Sn* 

977  +  8n 

978  +  8n 

979  +  8n 

980  +  8n 

981  +  8n 
582  +  8n 


COg  Reaction  Constants 


Molecular  Weight  of  C0g 

44 

1016  +  8a* 

Free  Stream  Mole  Fraction  of  COg 

From  performance 
calculation 

1017  +  8n 

Rate  Constant,  in. /sec 

0.138 

1018  +  8n 

Activation  Energy,  Btu/lb  mole 

33000 

1019  +  Sn 

Heat  of  Reaction,  Btu/lb 

4500 

1020  +  8n 

Mass  Diffusion  Coefficient,  in. "/ sec  0.0214 

1021  +  8n 

Product  Specific  Heat,  Btu/lb-R 

|  0.25 

1022  +  8n 

1 

TABLE  6 


EROSION  INPUT  DATA  FOR  PHETOLIC  REFRASIL 


ABLATIVE  WALLS 


Parameter 

V:lue 

- - - - 

Data  Input 
Location 

Melting  Mechanism 

4.6  x  Kf5 

Melt  Layer  Thermal  Conductivity,  Btu/in.-sec-R 

856  8n* 

Viscosity  Constant,  lb/in. -sec 

2.07  x  10-8 

857  +  8n 

Viscosity  Activation  Energy,  Btu/lb  mole 

Melt  Specific  Eeat,  Bta/lb-R 

2.32  x  105 

858  +  8n 

0.30 

859  +  8n 

Melt  Density,  lb/in.3 

0.08 

86p  +  8n 

Heat  of  Fusion,  Btu/lb 

100 

861  +  8n 

Melting  Temperature,  R 

5360 

862  +  8n 

Vaporization  Mechanism 

Diffusion  Coefficient,  ia.^/sec 

0.04 

896  +  8n* 

Vapor  Molecular  Weight 

12 

897  +  8n 

Vapor  Pressure  Constant,  psia 

5  x  108 

898  +•  8n 

Molar  Latent  Heat,  Btu/lbmole 

314000 

899  +  Sn 

Mass  Latent  Heat,  Btu/lb 

26000  . 

900  +  8n 

Density,  lb/in.A 

0.081 

901  *•  8n 

Specific  Heat  of  Vapor 

0.44 

902  +  8a 

TABLE  7 

EROSION  INFUT  DATA  FOR  SILICON  CARBIDE 


Parameter 

Value 

Data  Input 
Location 

Vaporization  Mechanism 

Diffusion  Coefficient,  in.  /sec 

0.01 

896  +  8n* 

•  Vapor  Molecular  Weight 

40  » 

897  +  8n 

Vapor  Pressure  Constant,  psia 

4.28  x  10* 

898  +  8n 

Molar  Latent  Heat,  Btu/lb  mole 

215,600 

899  +  8n 

Mass  Latent  Heat,  Btu/lb 

4900 

900  +  8n 

Density,  lb/in. ^ 

0.01 

901  +  8n 

Specific  Heat  of  Vapor 

0.2 

902  +  8n 

*Index  n  goes  from  0  to  one  less  than  the  number  of  regions 
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ESTIMATION  OP  EFFECTIVE  THERMAL  CONDI  "TIVITy 
Phenol ic/Carboa  Cloth 

The  effective  anisotropic  therms!  conductivity  of  carbon- cloth  reinforced 
ablatives  was  determined  by  comparing  experimental  vail  temperature  pro¬ 
files  measured  in  model  rocket  motor  firings  to  the  profiles  calculated 
by  2P-A8LATE  with  assumed  values  of  thermal  conductivity.  The  experi¬ 
mental  anta  used  for  comparison  were  ge-ierated  with  the  hardware  shown 
in  Fig.  17  .  It  consisted  of  an  ablative-walled  combustion  chamber  and 
an  ATJ  bard  throat  section  backed-up  by  porous  graphite  (Bef.  22). 

The  experimental  conditions  for  the  four  tests  employed  in  the  daia-match- 
ing  procedure  are  listed  in  Table  8.  The  four  tests  covered  reinforcement 
orientations  of  0,  60,  and  84*,  together  with  a  test  in  which  the  combus¬ 
tion  chamber  walls  consisted  of  the  nouablative  ATJ  graphite.  Considera¬ 
tion  of  three  reinforcement  orientations  permitted  the  definition  of  the 
directional  thermal  conductivity  of  the  ablatives  while  the  data  from  the 
passive  ATJ  graphite  wall  firing  established  the  hot-gas  boundary  conditions. 

TABLE  8 


BEN  CONDITIONS  FOB  OFg/MMH  FIRINGS  (EEF.  22)  SELECTED  FOB 
EVALUATION  OF  EFFECTIVE  PROPERTIES  FOR  PBTNOLlc/CARBON  CLOTH 


Test  Number 

1 

2 

3 

4 

Nominal  Chamber  Pressure,  psia 

110 

110 

110 

no 

Nominal  Mixture  Patio,  o/f 

2.0 

2.0 

2.0 

2.0 

Mainstage  Duration,  seconds 

300 

300 

300 

300 

Throat  Insert  Material 

ATJ 

ATJ 

ATJ 

ATJ 

Chamber  Insert  Material 

ATJ 

Phenolic/ 

Carbon 

Cloth 

Phenolic/ 

Carbon 

Cloth 

Phenolic/ 

Carbon 

Cloth 

Ablative  Material  Reinforcement 
Orientation  ^relative  to  radial 
coordinate),  degrees 

0 

60 

84 

*The  orientation  is  expressed  as  the  angle  between  the  reinforcement  and 
the  radial  vector  consistent  with  the  2-D  ABLATE  definitions. 
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Experimental  vail  temperature  profiles  in  the  combustion  chamber  vails 
vere  obtained  from  thermocouples  embedded  at  tbe  axial  locations  and 
depths  (from  the  hot-gas  surface)  given  in  Table  9.  The  temperature  his¬ 
tories  recorded  daring  the  four  motor  firings  are  listed  in  Table  10. 


The  overall  OP^KH H  experimental  program  from  vhich  the  data  of  tests 
1  through  4  vere  obtained  included  a  series  of  very  short  firings  in 
copper-vailed  chambers  under  conditions  of  chamber  pressure,  mixture 
ratio,  chamber  geometry,  and  injector  configuration  very  similar  to  those 
of  tests  1  through  4.  With  proper  instrumentation  and  reduction  of  cir¬ 
cumferential  and  axial  conduction,  the  transient  vail  temperature  data 
from  such  firings  permit  a  calorimetric  measurement  of  the  heat  fluxes 
to  the  copper  vails  tinder  conditions  of  high  gas-to-vall  temperature  dif¬ 
ference.  From  the  transient  heat  flux  data,  the  gas-side  heat  transfer 
coefficient  is  defined  by 


qA 

C  -  T 

aw  W 


(2-1) 


where  T  is  the  adiabatic  vail  temperature  determined  from  the  thermo- 
aw 

dynamic  gas  stagnation  temperature  by  propellant  performance  calculations 
at  the  bulk  mixture  ratio.  The  gas-side  heat  transfer  coefficient  profile 
defined  by  the  cepper-vall  actor  firings  (Ref.  16)  is  shown  in  Pig.  18 . 
The  shape  of  the  curve  in  Fig.  18  agrees  with  the  prediction  of  the  sitsp- 
lified  Bartz  equation  (Ref.  23). 


Based  upon  the  data  in  Fig.  18,  the  first  attests  to  match  experimental 
and  calculated  wall  temperature  profile  employed  a  heat  transfer  coeffic¬ 
ient  calculated  by  means  of  the  simplified  Bartz  equation  and  an  adiabatic 
wall  temperature  based  upon  the  thermodynamic  gas  stagnation  temperature 
at  a  mixture  ratio  of  2.0.  From  the  Rocketdyne  N-element  performance  pro¬ 
gram  (Ref.  24),  Tqv  was  equal  to  6200  F. 


Initial  trial  and  error  calculations  to  determine  effective  thermal  con¬ 
ductivities  were  conducted  with  TAP-4  computer  program  modified  for  ab¬ 
lative  calculations.  The  formulation,  method  of  calculation,  and  speci¬ 
fication  of  physical  properties  in  TAP-4  are  essentially  the  same  as  those 


Figure  10.  Topical  Variation  of  the  Heat  Transfer  Coefficient  With  Longitudinal  and  Cir¬ 
cumferential  Location  for  Self-Impinging  Doublet  Injector  Employed  in  Ablative 
Motor  Firings 


TABLE  9 


LOCATION  CV  THERMOCOUPLES  USED  TO  MEASURE  EXPERIMENTAL 
TEMPEEATURE  HISTORIES  USED  FOR  EVALUATION  OF 
EFFECTIVE  PROPERTIES  FOR  PHENOLIC /CARBON  CLOTH 
(Thrust  Chamber  Locations,  OFg/Amine  Tests  i  Through  4) 


Thermocouple 

Number 

Test  Number 

1  1  2 

3 

4 

Axial 

and  Radial  Coordinates, 

inches 

X 

Y 

X 

Y 

X. 

V 

X 

Y 

1 

5.23 

2.11 

2.87 

2.61 

2.87 

2.61 

2.87 

2.61 

2 

5.48 

2.94 

2.87 

2.16 

2.87 

2.16 

2.87 

2. 1  S 

3 

7.73 

2.11 

2.87 

2.93 

2.87 

2.93 

2.87 

2.93 

4 

7.98 

2.94 

2.87 

2.36 

2.87 

2.36 

2.87 

2.36 

5 

5.23 

to 

• 

* 

5.74 

2.61 

5.74 

2.61 

5.74 

2.61 

6 

5.48 

2.94 

5.74 

2.16 

5.74 

2.16 

5-74 

2.16 

1  7 

7.73 

2.11 

5.74 

2.93 

5.74 

2.93 

5-74 

2.93 

8 

7.98 

2.94 

5.74 

2.36 

5-74 

2.36 

5.74 

2.36 

9 

10.32 

2.94 

9.48 

2.46 

- 

- 

8.97 

1.90 

10 

11.32 

2.94 

9.48 

1.71 

- 

- 

8.97 

2.93 

11 

- 

- 

9.48 

2.93 

- 

- 

9.48 

2.32 

12 

- 

- 

10.32 

2.32 

- 

- 

9.48 

1.72 

13 

- 

- 

10.32 

1.57 

- 

- 

9.48 

2.93 

14 

- 

10.32 

2.93 

- 

10.32 

2.32 

15 

- 

11.38 

2.36 

- 

- 

10.32 

1.57 

16 

- 

- 

11.38 

1.76 

- 

- 

10.32 

2.93 

17 

-■ 

- 

11.38 

2.93 

- 

- 

- 

- 

$ 


I.  Test  1 


n=e, 

seconds 


Thermocouple  Kttrber 


res>erat ure,  F 


111.  Test  3 


Tiae, 
second?  ? 


40.0 

50.0 

60.0 

70.0 

80.0 

90.0 

100.0 

120.0 

140.0 

160.0 

180.0 

200.0 

220.0 

240.0 

260.0 


155 

,rr 

72  1 

215 

105 

280 

115 

345 

130 

410 

155 

470 

175 

535 

200 

665 

275 

790 

345 

920 

395 

1030 

435 

1120 

465 

1205 

490 

1280 

510 

1345 

525 

540 

|  '435 

550 

Thermocouple  Number 

"*  1  5  i  6  I  7 

Te  mpe rature,  F 


90  220 
95  330 

100  425 

110  510 


92  92 

7?Q  175 

425  400 

625  585 

790  730 

930  870 

1060  1010 
1170  1110 
1270  1205 

1560  1285 
1430  1360 

1545  1475 
1620  1580 
1680  1670 
1750  1750 
1810  1820 
1870  1890 
1945 
1995 
2045 
2075 


resx  2 


Tine, 

seconds 


7-5 

10.0 

2G.0 

30.0 

40.0 

50.0 

60.0 

70.0 

30.0 

90.0 

100.0 

120.0 

140.0 

160.0 

180.0 


Kit  HJ 


280.0  650 
3C0.0  685 


IV.  Test  4 


n  90 
75  90 

100  90 

140  90 

225  90 

295  95 

405  95 

505  100 
595  115 
665  135 
720  160 
845  210 
930  350 
1025  420 
1120  445 
1260  440 
1400  500 
1480  540 
1525  580 
1580  625 

1625  680 


Time 

seconds 


40.0 
50.0 
60.0 
70.0 
80.0 
90.0 
100.0 
120.0 
140.0 
160.0 
1S0.0  125 
200.0  150 
220.0  175 
240.0  195 
260.0  215 
280.0  230 
295.0 


25 

40  20 

65  25  20 

95  30  25 

135  35  50 

185  40  70 

235  45  85 

320  45  115 

340  55  170 

365  65  210 

445  85  250 

585  100  285 

705  115  320 

815  130  380 

905  145  460 

975  160  550 


[t. 


:  J 


1  2 


85  n  90 

85  75  90 

80  100  90 

80  140  oj 

80  215  90 

85  295  95 

95  405  95 

110  505  100 

150  595  115 

1 55  665 

195  720 

290  845 

360  950 

410  1025 

455  1120 

490  1260 

550  1400 

575  1480 

615  1525 

650  1580 

685  1625 


Thermocouple 

8  |  9 i  10 

■  1  l _ - _ 1 _ 


rerperatore,  2 


260  {  620 
795 


■urcif] 


1880  2560 


840  2500  j  2640 


1460  .  450 
151 

690  J 1555 


est  4 


TaercocGijpIe  Vradbe 


2 

J3 

L*J 

5 

6 

IB! 

8 

9  ! 

!  w 

25 

40  ! 

20 

65 

25 

95 

30 

135 

35 

185 

40 

235 

45 

320 

45 

340 

55 

365 

65 

445 

85 

585 

100 

705 

115 

S15 

130 

905 

165 

975 

l6u 

20 
25 

50  10 

70  15 

85  20 

i!5  25 

170  60 
210  85 
250  110 
285  125 
320  155 

>0  175 

460  215 

550  i  230 


10 

15 

35 

55 

75 

95 

20  135 
35  175 
55  220 
75  275 
95  335 

120  405 
125 
130 
135 


Tempe  rciure,  F 


1410  865 

0  1605  1025 

5  1770  1155 

5  I9ir  :  1280 

5  202J  I  1385 

5  2125  1485 

5  2215  1580 

5  2365  1740 

5  2495  1865 

0  2610  1995 

5  2710  2100 

5  2799  2200 

5  2855  2285 

2910  2355 

2960  2420 
2995  2480 
3025  2525 


paaloy-cd  in  2D-13LA22,  except  ties,  \%  is  ove— dixexsicxsJ  gag  s5ss 
isotreozc  is  all  parameters,  tbs  itexmscscplc  dsi-2  obtaioid  in  tie  coir- 
i^xstrom  chamber  vt*  t  ^xpeexen  ip  ns  *  S1 o  tr^r*— n^«Tg*- zxslvsis 

bscaase  •£»  variat>en  in  i&nasrefcsre  is  dsarisarted  bjr  rasaiiii  teaser— 

atere  gradients  in  this  regies. 

Snccfessri,  e  ccarpatear  rxss  ever  a  vide  range  of  radial  thesaai  ecofectivities 
did  not  produce  a  natch  vita  tie  experimental  temperature  profiles  for  an? 
reinfor-cenant  orxestatxen.  5&rtoT6r ,  vtsn  the  assumed  radial,  thamal  cts- 
^^tn-Sies  at  any  tvo  raSnforcaasnt  orientations  which  vseld  orocnce  a 
JEEtch  between,  composer— csicalersed  and  wp* r * ,al  char  dcoth*  at  these 
orientations  vere  combined  to  calculate  the  major  y—d  minor  thermal 
condg  ctiTitiea  according  to 

2  2 

2?  =  cos  6  +  ^  sin  b  (2-2) 

the  resulxsrt  definition  of  radial  conductivity  for  the  third  orientation 
did  sot  produce  a  natch  between  ccasater  calculation  exoerieaint. 


Several  computer  nms  vere  conducted  vith  2r»-i3L5rE  to  natch  the  terser— 
atoms  in  feat  1  vith  tie  Af-J  graphite  chamber  (rabies  8  and  10)  the 
calculated  surface  temperatures  in  the  throat  region  vere  close  to  the 
adiabatic  vail  temperature  (6200  P)  rather  than  the  values  of  >009  to 
>590  P  measured  vith  the  optical  pyrometer. 

Because  of  the  discrepancies  between  calculated  and  measured  temperatures, 
the  gas-side  boundary  conditions  vere  reviewed  vith  particular  attention 
given  to  possible  injector  effects  on  thrnst  voadber  heat  transfer.  The 
injector  employed  in  the  test  hardware  was  particularly  designed  for  a 
low-  tsixture  ratio  near  the  charier  walls  so  a3  to  reduce  the  concentra¬ 
tion  of  water  vapor  belov  the  corrosive  limit .for  carbon  vails.  Sampling 
of  the  injector  3pray  pattern  under  cold  flow  conditions  had  demonstrated 
that  the  mixture  ratio  near  the  vail  was  reduced  to  less  than  1.0  (Hef.  16). 
At  this  mixture  ratio,  the  thermodynamic  gas  temperature  is  >000  P.  Un¬ 
published  gaseous  mixing  experiments  at  Bonketdyne  have  shown  that  parallel 
gas  streams  have  not  approached  equilibrium  after  as  much  as  22  equivalent 
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atresn  diameters.  Therefore  r  it  is  reasonable  to  expect  that  the  combustion 
gas  is  racist.  chambers  has  a  local  iesperstare  defined  by  the  mixture 
ratio  distribsticn  si  the  injector. 


2  or  tie  injector  employed  in  the  ablative  no  tor  firings,  the  adiabatic 
vail  tesoeratare  vss  therefore  hypothesized  to  be  ach  leas  than  the 
thermodynamic  tempera  tare  for  the  balk  nix  tare  ratio  as  assumed  for  regen- 
eratire  cooling.  She  increase  is  measured  vail  temperatures  from  the 
injector  to  the  throat  for  the  data  of  feet  1  also  indicated  that  the 
zdiahaiie  vail  temperature  vas  a  fanatics  of  axial  distance. 

Bealistic  vaices  for  the  adiabatic  vail  temperature  and  the  gas-side  heat 
transfer  coefficient  in  the  test  hardware  vere  derived  by  an  iterative 
procedure,  in  adiabatic  vail  temperature  distribution  vas  first  assured. 

A  distribution  of  heat  transfer  coefficients  vas  calculated  so  as  to  re¬ 
sult  in  the  experimental  heat  flux  profile  to  copper  chamber  vails  at 
ZOO  F.  “chese  adiabatic  vail  temperature  and  film  coefficient  distribu¬ 
tions  vere  then  employed  vith  2P-ABIAfg  to  calculate  the  vail  temperature 
profile  of  Test  1.  If  a  aatch  vas  not  obtained  vith  experimental  values , 
a  nev  set  of  values  for  adiabatic  vail  temperatures  vas  postulated  and  a 
nev  iteration  vas  conducted,  The  procedure  vss  continued  until  calculated 
and  experimental  vail  temperatures  converged.  The  final  distributions  of 
adiabatic  vail  temperature  and  heat  transfer  coefficient  are  given  in  Fig.  19* 
The  gas-side  boundary  conditions  shown  in  Fig.  19  produced  the  vail  tempera¬ 
ture  histories  shove  as  dotted  curves  in  Fig.  20.  The  corresponding  experi¬ 
mental  temperature  histories  are  superimposed  thereon  as  solid  curves 
shoving  reasonably  good  agreement  between  computed  values  and  experimental 
results . 


On  the  basis  of  the  revised  gas-side  boundary  conditions,  nev  sets  of  ablative 
vail  temperature  profiles  vere  derived  vith  the  TAP-4  computer  program  for 
assumed  values  of  thermal  conductivity.  Because  of  the  relatively  low  adia¬ 
batic  vail  temperatures  involved,  the  range  of  calculated  vail  temperatures 
vas  relatively  narrov  for  carbon  cloth  ablatives.  Therefore,  it  was  consid¬ 
ered  impractical  to  define  a  temperature  dependency  for  the  thermal 
conductivities  derived.  Examples  of  the  match  between  experimental  and 
calculated  temperature  profiles  are  given  in  Fig.  21. 
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Toe  curves  represent  coapuSBth-rosnlifi.^  the  TAP-4  program  for  various 
values  of  thermal  conductivity.  Further  refinement  of  the  computed  results 
was  subsequently  accomplished  vith  the  2B-ABLATE  program.  It  is  because  of 
the  relatively  low  wall  temperatures  (both  calculated  and  experimental) 
shown  in  Fig.  21,  that  no  verification  of  the  high  temperature  thermal  con¬ 
ductivity  of  the  ablative  walls  could  be  made.  Hocket  motor  firings  giving 
verified  adiabatic  wall  temperatures  of  at  least  5000  F  under  well- 
instrumented  conditions  are  required  to  obtain  such  high  temperature 
information. 


The  comparison  between  calculated  and  experimental  profilea  for  the  three 
reinforcement  orientations  was  combined  with  that  between  the  experimental 
and  calculated  char  depths*  and  weighted  with  the  requirement  that  the 
radial  thermal  conductivities  for  these  orientations  be  related  by  Eq.  2-2. 
Based  on  all  considerations,  a  constant  value  of  8.0  x  10-^  Btu/in-sec-F 
is  recommended  for  the  major  direction  thermal  conductivity  together  with 
a  minor  direction  thermal  conductivity  which  ranges  from  7.0  x  10~^  to 
1.5  *10  5  Btu/in.-sec-F  as  given  in  Table  4. 


Phenolic  Refrasil 


The  approach  used  to  estimate  the  effective  thermal  conductivity  of 
phenol ic-refrasil  differed  markedly  from  that  employed  for  phenolic- 
carbon  cloth.  Because  the  effective  thermophysical  and  thermochemical 
properties  of  this  system  had  been  estimated  under  another  contract  with 
the  2D-CHAR  (isotropic)  computer  program  (Ref.  7),  their  validity  had  to 
be  checked  with  the  (anisotropic)  2D-ABLATE  computer  program. 

The  experimental  data  used  for  comparison  were  those  obtained  on  certain 
qualification  tests  of  the  Apollo  Command  Module  rocket  engine  and  nozzle 
extension.  The  computer  model  of  this  test  configuration  is  shown  in 
Fig.  22.  The  specific  data  employed  were  the  backwall  temperature 


*Erperiments!  char  depths  were  based  upon  photographs  of  sectioned  chamber 
walls.  Calculated  char  depths  are  based  apon  the  1000  F  isotherm. 


BACKUP  MATERIAL 


histories  during  e  given  mission  duty  cycle  at  tvo  external  points.  The 
axial  locations  (Fig.  22)  correspond  to  (l)  the  nozzle  throat,  and  (2/ 
a  point  1  inch  downstream  of  the  upstream  face  of  the  nazor  counting 
flange. 

Because  of  the  linitations  of  2B~€E£S  (mentioned  above),  the  original 
effective  properties  vere  estimated  vithcuct  consideration  of  anisotropicity 
Additionally,  the  reinforcement  orientation  in  the  chamber/ extension 
(45  degrees  relative  to  the  chamber  centerline)  differed  from  that  in  the 
nozzle  (parallel  to  the  nozzle  centerline).  Two  separate  calculations 
vere  made  to  produce  a  match  between  the  expert rental  and  predicted  back- 
vali  temperature  histories.  The  first,  in  vhich  the  reinforcement  orient¬ 
ation  in  the  entire  motor  vas  assumed  to  be  45  degrees  and  the  second,  in 
vhich  it  vas  assumed  to  be  0  degrees.  With  these  assumed  orientations, 
the  effective  thermophysical  and  thermo chemical  properties  vere  estimated 
by  separately  matching  the  throat  station  and  nozzle  station  data  (Bef. 

7) .  This  approach  vith  the  2B-CHAB  vas  required  because  that  program  vas 
limited  to  isotropic  properties  and  could  accomodate  one  charring  mater¬ 
ial  only.  The  property  variation  resulting  from  a  different  orientation 
of  an  anisotropic  material  is  in  effect  a  second  charring  material. 

Alter  resolving  the  anisotropic  properties  estimated  vith  2D-CHAB  into 
the  major  and  minor  directions,  the  properties  as  previously  estimated 
vere  submitted  on  2B-A3LATE.  The  results  are  shown  in  Fig.  23  and  24  . 
The  experimental  Apollo  data  as  veil  as  the  backvv.ll  temperature  history 
predicted  by  2D-CHAR  are  superimposed.  The  differences  are  nominal;  e.g., 
nowhere  exceeding  about  20  F  at  the  throat  station.  The  effective  thermo- 
physical  and  thermochemical  properties  thus  estimated  vere  presented  in 
Table  4.  Again  it  should  be  remembered  that  temperature  measurements  in 
th«s  range  shown  in  Fig,  24  and  25  are  not  a  verification  of  ablative 
thermal  conductivities  at  high  (>2500  F)  temperatures. 
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Figure  23  .  Throat-Station  Backwall  Temperature 
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TASK  5,  DESIGNERS  GUIDE 


The  tvo-diseasioasl  computer  pro  gran  is  ideally  suited  for  the  calculation 
of  char  rate,  erosion,  and  temperature  distribution  viihin  curved  boundaries 
mliissiterial  configurations  typical  of  real  thrust  chanter  valla  if 
all  the  ays  ten  parameters  required  for  program  input  are  known.  However , 
the  use  of  the  computer  program  to  optimize  an  ablative  vail  configuration 
within  the  constraints  of  an  overall  engine  design  requires  an  iterative' 
procedure  in  which  the  parameters  not  completely  fixed  by  mission  require¬ 
ments  arc  systematically  varied.  Because  of  the  large  number  of  possible 
parameters  io  be  considered,  a  systematic  ablative  design  requires  that  the 
approximate  effects  of  the  most  important  parameters  be  known  in  advance; 
otherwise  the  number  of  iterations  becomes  excessive. 

Subsequent  sections  of  this  report  provide  the  designer  with  information 
relating  the  useful  life  of  ablative  chamber  walls  to  the  more  important 
operating  variables  and  material  properties.  Discussio  is  restricted  to 
the  following  thermal  responses:  (l)  char  rate,  (2)  thermal  penetrr  .ion, 
and  (3)  surface  erosion  and  surface  temperature.  Other  factors  in  mechan¬ 
ical  design  which  are  related  to  thermal  history  such  as  material  strength 
and  thermal  stress  are  nnt,  considered.  Discussion  is  also  restricted  to 
the  phenolic  resin  systems  and  to  the  coianon  reinforcements,  refrasil  and 
carbon  cloth. 


SIGNIFICANT  PARAMETERS 

The  variables  known  to  be  important  to  ablative  chamber  design  can  be 
divided  into  three  categories: 

1.  Combustion  gas  properties  and  operating  conditions 

2.  Chamber  geometry 

3.  Wall  materia'  properties 


Ice  prop-ill  act  ecuEbinsties  sad  aixtare  ratio  defifc;  the  gas  teroa  rasare 
and  such  properties  as  specific  beat,  molecular  weight,  sad  viseosl'Sy 
which  is  taro  define  the  beat  transfer  coefficient  gad  the  coseestratian 
of  rcrrcsiTe  gaseous  species  (EgO,  ^)>  The  fcdiahstic  ■wall 

tespercmre  and  the  effective  gas  preparties  are  specified  by  the  aixtcre 
ratio  near  the  wall.  The  chamber  operating  conditions  which  affect 
ablative  charier  design  are  duty  cycle ,  ehanber  preecore,  and  thrast 
level.  Duty  cycle  specifies  not  only  the  total  firing  duration  bat  also 
the  available  cooldown  periods  between  successive  firings.  To  the  extent 
that  sensible  heat  stored  in  a  char  layer  or  refractory  throat  insert  can 
be  lost  ay  radiation  between  firing  periods,  the  length  of  the  poises 
and  the  interposed  cooldown  periods  can  reduce  char  depth  *na  surface 
erosion. 

Chamber  pressure  determines  the  heat  transfer  coefficient  and  the  concen¬ 
tration  of  corrosive  species  in  the  eosbnstion  gases.  With  other  conditions 
fixed,  the  magnitude  of  the  heat  transfer  coefficient  determines  hov  closely 
the  surface  temperature  approaches  the  adiabatic  wall  temperature.  Over- 
vide  ranges  of  heat  transfer  coefficient,  variations  in  A"?  between  the 
surface  and  adiabatic  wall  temperatures  produce  relatively  snail  percen¬ 
tage  changes  in  the  AT  between  surface  and  charring  reaction  zone  and 
therefore  relatively  minor  differences  in  char  rate.  However,  because  of 
the  Arrhenius  type  of  temperature  dependence  for  surface  reactions,  the 
differences  in  surface  temperature  produced  by  changes  in  heat  transfer 
coefficient  can  produce  significant  differences  in  surface  erosion.  Thrust 
level  combined  with  chamber  pres  rare  and  contraction  ratio  determines  the 
throat  and  chamber  dimensions  and  thereby  the  two-dimensional  nature  of 
the  axisymaetric  system. 

The  obvious  variables  of  chamber  geometry  are  the  expansion  and  contraction 
ratios,  the  contraction  and  expansion  l/2-angles  and  the  arrangement  of 
hard  or  soit  throat  inserts  with  multilayer  ablative  walls  of  different 
reinforcements  and  reinforcement  orientations.  Th.*  effects  of  these  vari¬ 
ables  are  tvo-dimenbional  in  nature  and  many  of  the  boundary  conditions 
are  difficult  to  specify  for  a  real  system. 
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It  part  icalzriy  difficali  to  specify  lie  beat  transfer  eoefticient  for 
his A  csatrar-tico  ratio's  arc  high  eetriraciica  AS-argles,  amd  ike  resraliarsi 
erotica  rates  art  cgrrespcsdSrgly  Ereeriain.  sit  ezpzzsiea  ratio  deter- 
aiaty  the  radiative  Tiev  fsciar  bsiveen  tie  ikrcai  regiea  ar>3  fret  space 
and  osriiaaly  ctfices  tit  max  9 min,  ienggraiarc  of  ssrsalziive  threat  inserts- 

The  most  iaspsriani  vail  mzierizl  psrareters  are  iierrgal  coarsctiYity,  r*32n 
cGsteoi,  eheaicai  reaciicms  between  char,  pyrolysis  gases  and  reioforcenant. 
and  sr~  face  erosics  resistance-  The  items!  ecn&ciiviiy  of  the  cfear  eoctrels 
the  diffusion  of  krai  frets  tie  ecscbasiics  g2«es  to  tie  charring  reaction 
zone.  High  thermal  ecndreiivity  produces  £  high  charring  rate  hat  also 
produces  a  lover  sarfaee  bfisneraiure  and  therefore  a  lever  surface  erosion. 

The  thermal  ccodactivify  o?  the  virgin  ablative  determines  the  thermal 
penetration  through  the  virgin  ablative-  Hesin  content  determines  the  rate 
of  char  front  advance  as  do  the  energy  absorbing  processes  occurring  in  the 
char  layer  between  the  char  and  the  percolating  pyrolysis  gases-  Although 
the  reaction  be tveen  refrasil  (silica)  reinforcement  and  char  has  been 
proposed  as  a  high  best  sink,  Sockeid yna  experience  on  the  Apollo  program 
indicates  that  this  reaction  does  not  ©eeur. 

The  resistance  to  erosion  of  ablative  or  passive  refractory  isaieriala  is 
dependent  upon  a  high  melting  point,  lev  vapor  pressure,  and  c crpetabi  1  iiy 
vitfc  boundary  layer  gas  species.  The  optical  c  orb  instil  of  lev  erosion, 
lev  char  rate,  and  lov  thermal  penetration  vill  usually  require  a  rslti- 
layer  combination  of  ablative  and  passive  materials. 

Although  the  beat  transfer  in  an  ablative  systen  is,  in  fact,  ivo  dioen- 
sional,  tbe  radial  temperature  gradients  through  n<fst  ebariber  vjlls  are 
each  larger  than  tbe  axial  gradients.  £or  convenience  in  prelioinary 
design,  axial  conduction  can  be  ignored  and  a  "first  cut"  design  can  he 
based  on  radial  heat  transfer  only.  Similarly,  it  is  most  convenient  to 
evaluate  the  effects  of  the  various  physical,  cbeaieal,  end  transport 
properties  of  the  vail  tsaterial  and  combustion  gas  on  a  one— dimensional 


hssiso  2a  the  follcviEg  section  of  this  report,  the  thermal  responses  of 
ablative  vail 3  to  a.  one-diisgasiosal  racial  heat  flux  are  presented  graph¬ 
ically  and  analyzed  for  the  following  parameters: 

Vail  the  real  conductivity 
Adiabatic  vail  tecierafere 
Essin  fraction 
Heat  transfer  coefficient 
Chamber  radios 
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Definitions 

Erperirsents  (Ref.  ll)  have  shovn  that  charring  of  phenolic  resins  occurs  at 
a  rapid  hot  finite  rate  over  the  approximate  temperature  range  300  F  to 
1600  F.  Therefore,  a  narrow  zone  of  partially  charred  resin  rather  than 
a  distinct  char-virgin  interface  exists  in  ablative  walls  both  during  and 
after  burn  time.  However,  because  of  the  thermodynamic  description  of 
the  char  reaction  in  the  computer  program  it  is  convenient  to  define  a 
char  front  for  purposes  of  comparing  various  parametric  effects.  In  the 
following  sections  of  the  report,  char  rate  is  defined  as  the  rate  of 
advanc-.  of  the  1000  F  isotherm.  Similarly,  char  depth  is  the  distance 
from  the  hat-gas  surface  to  the  1000  F  isotherm. 

Transfer  of  heat  from  ablative  walls  to  outer  portions  of  the  vehicle  is 
directly  related  to  the  temperature  profile  in  the  virgin  material.  This 
can  be  loosely  expressed  as  thermal  penetration  through  the  wall.  In  the 
following  sections,  thermal  penetration  is  defined  as  the  location  of  the 
100  F  isotherm  with  respect  to  the  surface  for  a  system  initially  at  70  F. 

Surface  erosion  is  the  increase  (assumed  circumferentially  uniform)  in 
the  diameter  of  the  chamber  at  a  given  location. 


Surface  temperature  is  the  temperature  of  the  ablative  wall  exposed  to 
the  combustion  gas.  This  response  is  directly  related  to  surface  erosion 
and  in  some  cases  is  used  in  place  of  erosion  for  evaluating  parametric 
effects . 


Thermal  Conductivity 


The  effect  of  ablative  wall  thermal  conductivity  on  char  depth  is  shown 
in  Fig.  25  through  27  for  adiabatic  wall  temperatures  of  3500,  5000, 
and  7000  F  respectively.  All  curves  are  based  upon  a  chamber  radius  of 
1.90  inch  and  a  wall  thickness  of  1.20  inch.  The  3500  F  curves  are  based 
upon  a  heat  transfer  coefficient  of  G.OOO65  Btu/in^-sec-F  and  a  resin  content 
of  0.30;  the  50Q0  and  7000  curves  are  based  upon  a  heat  transfer  coefficient 
of  O0OOO5  Btu/in  -scc~F  and  a  resin  content  of  0.25*  All  three  graphs  show 
that  char  rate  is  a  strong  function  of  thermal  conductivity.  If  a  modified 
Biot  modulus  is  defined  as: 

Bi  -  hg  (3-1) 

kc 

where  A  y  is  the  char  depth,  a  simple  expression  can  be  written  for  the 
c 

effect  of  thermal  conductivity  on  char  rate  which  applies  for  Biot  moduli 
greater  than  about  10.  This  simple  expression  which  can  be  employed  to 
interpolate  between  curves  is 

■J5" 


(Ayc)2 


The  curves  of  Fig.  25  through  27  are  for  heat  transfer  coefficients  in 
the  range  of  0.0005  to  O.GOO65  Btu/hr-ft^-F  typical  for  a  combustion 
chamber  at  chamber  pressures  from  70  to  150  psia.  The  relationship  between 
char  depth  and  thermal  conductivity  for  very  high  and  very  low  values  of 
heat  transfer  coefficient  i&  given  in  Fig.  28  .  For  high  values  of 
film  coefficient,  Eq.  3-1  and  3-2  still  apply  for  interpolation.  However, 
char  depth  is  relatively  independent  of  thermal  conductivity  at  very  low 
values  of  film  coefficient  (Fig,  28). 
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CONDITIONS 
hg  =  0.00065 
Fore  =  0.30 


1.6  X  jO*5 


k  =  9  X  10 
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Figure  25  .  Effect  of  Thermal  Conductivity  on  Char  Front  Advance 
1,20  — . -  . ?atic  Vail  Temperature  of 

- - - -  I  k  =  2.1  X  10"^ 

CONDITIONS  / 

1.00  -  - /  I 

r  hg  ■  0.0005  / 

FRES  “  0,25  / 


k  -  7-0  X  10 
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Figure  26  .  Effect  of  Thermal  Conductivity  on  Char  Front  Advance 
For  an  Adiabatic  Wall  Temperature  of  5000  F 
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The  effect  of  thermal  conductivity  on  thermal  penetration  is  shown  in 
Pig.  29.  As  weald  be  expected,  thermal  penetration  depends  upon  the 
conductivity  of  both  the  char  and  the  virgin  materials.  The  results  shewn 
in  Fig.  29  can  be  extended  to  other  adiabatic  wall  temperatures  and  other 
char-  and  virgin-region  thermal  conductivities  by  means  of  the  relation: 


In  Eq.  3-3,  A y is  the  thermal  penetration  depth  (the  100  F  isotherm) 

while  Ay  is  the  char  depth  and  k  is  fche  thermal  conductivity  of  the 

virgin  material.  The  char  depth  can  he  obtained  from  Fig.  25  through  27 

while  a  reference  Ay  is  obtained  from  Fig.  29. 

? 


'Hie  effect  of  thermal  conductivity  upon  the  gas-side  surface  temperature 
is  shown  in  Fig.  30  for  an  adiabatic  wail  temperature  of  5000  F  and  a 
heat  transfer  coefficxent  of  0.0005  Btu/in  -sec-F.  The  effect  of  increasing 
thermal  conductivity  is  to  reduce  the  rate  at  which  the  surface  temperature 
approaches  the  adiabatic  wall  temperature  or,  in  other  words,  to  lower  the 
surface  temperature  at  any  time.  Interpolation  between  the  curves  of 
Fig.  30  and  extrapolation  to  other  adiabatic  wall  temperatures  can  be 
made  with  fair  accuracy  by  the  relationship: 


where  the  effective  temperature  of  the  heat  sink  is  taken  to  be  1200  F. 
Eq.  3-4  nay  be  used  ^hen  the  modified  Biot  number  defined  by  Eq.  3-1  is 
greater  than  10. 

In  general,  all  three  responses,  char  rate,  thermal  penetration,  and 
surface  temperature  are  related  to  the  square  root  of  the  thermal  conduc¬ 
tivity.  These  thermal  conductivity  effects  with  ablation  are  therefore 
very  similar  to  the  relationships  for  ordinary  transient  conduction.  In 
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using  the  curves  of  Fig.  25  through  3 0  .  a  mean  thermal  conductivity 
for  the  char  layer  should  be  estimated. 

For  lev  thermal  conductivity  material.,  in  vfaich  the  char  rate  ie  relatively 
low  and  the  radius  vector  does  not  change  drastically  with  time,  the  char 
depth  is  knovn  to  increase  approximately  as  the  square  root  of  time.  This 
is  true  of  the  curves  of  Mg.  25  through  27  for  turn  times  greater  than 
approximately  100  seconds. 


Adiabatic  Wall  Temperature 

Bates  of  char  front  advance  at  T  of  3500,  5000,  6000,  and  7000  F  are 

SV 

compared  in  Fig.  31  •  The  curves  apply  for  a  heat  transfer  coefficient 
of  0.0005  Btu/in2-sec-F,  a  resin  content  of  0.25  and  a  thermal  conductivity 
assumed  to  vary  linearly  from  2  x  10~^  to  3  x  10~^  Btu/in-sec-F  between 
2500  and  7000  F.  As  shown  (Fig.  31  )»  the  char  depth  for  a  burn  time  of 
150  seconds  is  almost  50  percent  greater  for  an  adiabatic  wall  temperature 
of  7000  F  than  for  a  temperature  of  3500  F. 

Even  on  0  one-dimensional  basis,  the  effect  of  adiabatic  wall  temperature 
on  char  rate  is  considerably  more  difficult  to  generalize  than  the  effect 
of  thermal  conductivity.  The  view  factor  of  the  surface  to  low-temperature 
regions  (free  space)  determines  the  degree  by  which  radiation  can  reduce 
the  heat  flux  into  tl.e  wall  and  thereby  reduce  the  char  rate.  Similarly, 
surface  reactions  can  absorb  significant  amounts  of  heat  to  reduce  char 
rate,  but  the  resultant  erosion  will  remove  insulating  char  and  thereby 
increase  char  rate.  The  curves  of  Mg.  31  are  for  negligible  radiation 
and  surface  reactions. 


An  approximate  method  for  calculating  the  effect  of  adiabatic  wall  temper¬ 
ature  upon  char  depth  is  to  apply  the  formula 


CHAR  DEPTH,  INCH 


where  the  enthalpy  parameter  lb  is  obtained  from  Fig.  J2  .  Che  enthalpy 
parameter  is  defined  by 


aw 


-  1200 

JH 


(5-6) 


where  AH  is  the  total  heat  absorbed  by  resin  pyrolysis,  gas  cr^eki’.ig, 
and  sensible  heat  absorption  by  the  pyrolysis  gases.  Equation  3-5  applies 
for  modified  Biot  numbers  greater  than  10  and  for  negligible  radiation 
and  surface  erosion. 

The  effect  of  adiabatic  wall  temperature  upon  thermal  penetration  is  shown 
in  Fig.  53  .  Coaqiarison  of  Fig.  33  and  31  indicates  that  the  effect  of 
adiabatic  wall  temperature  on  thermal  penetration  can  be  equated  to  the 
differences  produced  in  char  rate.  Under  steady  firing  conditions,  the 
difference  between  the  thermal  penetration  front  and  the  char  front  is 
effectively  independent  of  adiabatic  wall  temperature  and  is  a  function 
only  of  virgin  thermal  conductivity  and  time. 

The  relation  between  surface  temperature  and  adiabatic  vail  temperature 
is  shown  in  Fig.  34 .  The  surface  temperature  is  shown  to  gradually 
approach  the  adiabatic  wall  temperature.  The  curves  of  Fig.  34  can  be 
effectively  normalized  into  a  single  carve  by  plotting  the  ratio 


l  ^ 


L-i 


T  -  T 
aw  s 


T  -  1200 
aw 


vs  time  as  shown  in  Fig.  35.  The  results  shown  in  Fig.  35  provide 
the  justification  for  employing  Eq.  3-4  to  account  for  thermal  conductivity 
and  adiabatic  wall  temperature  changes  from  values  employed  in  Fig.  30. 


Hesin  Fraction 


The  effect  of  variable  phenolic  resin  content  upon  the  char  rate  of  abla¬ 
tive  walls  under  conditions  of  high  thermal  conductivity  and  high  heat 


SURFACE  TEMPERATURE ,  F  THERMAL  PENETRATION*  INCH 


7000  {•' 


transfer  coefficient  is  shown  in  rig.  yS  .  Yhs-  tvo  resin  fractions 
0.25  and  0.45,  represent  practical  iisit?  for  this  variable.  In  rig.  -57 
the  char  rates  for  the  tvo  resin  fractions  are  compared  ande?  conditions  of 
lav  therral  conductivity  and  lev  heat  transfer  coefficient.  As  expected, 
the  larger  heat  capacity  of  the  higb-resin-content  ablative  eaterial  re¬ 
sults  in  a  lover  char  rate.  The  effect  of  resin  content  oa  char  rate  if 
not  amenable  to  generalization  under  high  heat  transfer  conditions  because 
of  the  coupled  effect  of  chamber  radius  on  overall  heat  capacity  and  radial 
char  conductance,  (k/Ay  ).  Under  conditions  of  lover  heat  fluz,  an-?  partic¬ 
ularly  for  conditions  in  vhich  the  modified  Biot  ncsber  defined  by  &j.  5-1 
is  greater  than  10,  the  effect  of  resin  content  on  char  rate  can  be  approx¬ 
imated  by 


2  _  /FB  1 
^  yc^  1  *1  7 R  2 


(3—7)' 


The  effects  of  resin  content  on  thermal  penetration  are  shown  in  Fig.  53 
For  high  conductivity  ablatives,  the  effect  of  resin  content  on  therma.. 
penetration  is  negligible.  For  low  conductivity  ablatives,  higher  resin 
content  reduces  thermal  penetration  by  the  sane  degree  to  which  it  reduces 
the  char  rate.  Under  steady  firing  conditions,  the  difference  between 
the  thermal  penetration  front  and  the  char  front  is  independent  of  resin 
fraction  and  depends  only  on  the  virgin  thermal  conductivity  and  the  burn 
time. 


The  surface  temperature  histories  of  ablative  walls  with  resin  contents 
of  0.25  and  0,45  for  an  adiabatic  wall  temperature  of  5000  F,  a  heat 
transfer  coefficient  of  0.0005  Btu/in  -sec-F  and  an  average  char  thermal 
conductivity  of  7  x  10  Btu/in-sec-F  are  compared  in  Fig.  39  .  The 
resin  fraction  is  shown  to  have  a  negligible  effect  on  surface  temperature. 
A  similar  result  is  obtained  over  wide  ranges  of  thermal  conductivity, 
adiabatic  wall  temperature,  and  gas-side  heat  transfer  coefficient.  How¬ 
ever,  close  to  the  point  of  complete  depletion  of  resin,  the  surface 
temperature  deflects  upwards.  The  lover  the  initial  resin  content,  the 
earlier  is  the  inflection  of  the  associated  temperature  history. 
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Figure  39  .  Effect  of  Resin  Fraction  or  Surface 
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97 


Heat  Transfer  Coefficient 


The  effect  of  heat  transfer  coefficient  upon  char  rate  through  high 
conductivity  and  low  conductivity  ablatives  is  shown  in  Fig,  40  and  41, 
respectively.  Through  a  high  conductivity  vail,  the  effect  of  an  increasing 
heat  transfer  coefficient  is  to  significantly  increase  the  char  rate.  In 
fact,  for  lov  values  of  the  modified  Biot  number  defined  by  Eq.  3-1,  the 
increase  in  char  rate  is  more  than  proportional  to  the  increase  in  heat 
transfer  coefficient.  For  low  conductivity  ablatives  (Fig.  41  ),  the 
effect  of  heat  transfer  coefficient  becomes  negligible  for  heat  transfer 

.  ft 

coefficients  above  0.001  Btu/in  -sec-F.  In  terms  of  chamber  operating 
variables,  this  means  that  char  rate  is  essentially  independent  of  chamber 
pressure  for  chamber  pressures  above  approximately  200  psia. 

The  effect-  of  heat  transfer  coefficient  on  the  surface  temperature  of  a 
low  conductivity  ablative  is  shown  in  Fig.  42  .  For  low  conductivity 
ablatives,  ana  with  negligible  endothermic  surface  reactions  and  surface 
radiation,  the  surface  temperature  closely  approaches  the  adiabatic  wall 
temperature  for  high  or  moderate  heat  transfer  coefficients.  In  fact,  for 
chamber  pressures  above  200  psia,  the. ablative  surface  temperature  for  low 
conductivity  ablatives*  will  approximate  the  adiabatic  wall  temperature 
within  100  F  in  the  combustion  chamber.  Ho  simple  general  statement  can 
be  made  with  regard  to  high  thermal  conductivity  ablatives. 

The  effect  of  heat  transfer  coefficient  upon  thermal  penetration  through  a 
low  conductivity  ablative  is  shown  in  Fig.  43.  As  with  the  adiabatic  wall 
temperature  and  the  resin  content,  the  effect  of  heat  transfer  coefficient 
is  to  shift  the  thermal  penetration  front  by  a  distance  equivalent  to  the 
shift  in  char  front. 

C harnber  Radius 


The  effect  of  chamber  radius  on  char  rate  is  shown  in  Fig.  44 «  For  low 
values  of  ablative-wall  thermal  conductivity  and  heat  transfer  coefficient, 
there  is  essentially  no  difference  in  the  char  depth  between  chamber  radii 

"^PhanolTcTrefraa i 1  or  phenolic- graphite  at  orientations  greater  than  75 
degrees. 
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Figure  44.  Effect  of  Chamber  Radius  on  Char  Front  Advance 


of  1.90  and  10.0  inches.  For  high  values  of  ablative  wall  thermal  conduc¬ 
tivity  and  heat  transfer  coefficient;  the  char  depth  for  a  chamber  radius 
of  1.90  inches  is  less  than  for  a  chamber  radius  of  10.0  inches.  The  dif¬ 
ference  represents  the  competing  effects  of  the  greater  heat  sink  per  unit 
area  of  hot-gas  surface  and  the  greater  average  normal  area  to  conduction 
heat  transfer  per  unit  hot-gas  surface  for  the  sr-all  chamber  radius.  Because 
the  correction  for  chamber  radius  is  small,  the  approximation  resulting 
from  solution  of  the  implicit  relation 
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(3-8) 


may  be  employed  for  low  thermal  conductivity  ablatives.  To  correct  the 

parametric  curves  of  this  manual  for  chamber  radius,  a  value  of  1.90  inches 

should  be  employed  for  R  . 

ox 


The  effects  of  chamber  radius  on  surface  temperature  aad  thermal  penetration 
are  negligible  for  ablative  materials  and  need  not  be  shown.  This  is  not 
necessarily  the  case  for  passive  refractory  surfaces,. 


NOZZLE  THROAT  EROSION 


Surface  erosion  in  ablative  rocket  engine  vails  is  important  primarily  in 
tbe  nozzle  throat  region  because  the  throat  area  increase  produced  by 
recession  may  reduce  chamber  pressure  which  in  turn  effects  thrust,  specific 
impulse,  and  propellant  flovrates.  Some  estimation  of  throat  lifetime  is 
therefore  desirable  before  a  preliminary  nultimaterial,  rocket-engine  wall 
configuration  should  be  specified  for  analysis  by  2B-ABLATE. 

Experimental  results  (fief.  25  )  of  firings  with  the  propellant  combination 
of  NT0/5Q/S  "  50J&  UDHH  (flame  temperature  =  5180  F)  indicate  that 
refrasil-reinforced  ablative  materials  are  poor  choices  for  throat  walls 
for  chamber  pressures  of  100  psia  or  higher.  In  environments  which  do  not 
contain  water  vapor,  phenolic/carbou-cloth  soft  throats  might  be  considered 
for  low  chamber  pressures.  However,  hard  graphite  throats  would  be  pre¬ 
ferable  from  strength  considerations. 

The  consEon  materials  for  hard  throat  inserts  are  silicon  carbide  for  pro¬ 
pellant  combinations  such  as  NT0/5t$  NgH^  -  50/£  UDMH  which  contain  appre¬ 
ciable  amounts  of  water  vapor  and  graphite  for  fluorinated  oxidizer 
systems  which  contain  little  or  no  water  vapor.  In  noneorrosive  gas 
streams,  graphite  is  able  to  withstand  higher  temperatures  than  silicon 
carbide  which  vaporizes  because  of  decomposition  at  temperatures  in  the 
neighborhood  of  4700  F  (Ref.  2l) » 

The  erosion  mechanisms  for  ATJ  graphite  provided  for  in  2D-ABLA3E  in¬ 
clude  vaporization  and  chemical  reaction.  Vaporization  is  a  function  of 
adiabatic  wall  temperature  together  with  heat  and  mass  transfer  coeffi¬ 
cients  which  jointly  determine  the  surface  temperature  (vapor  pressure) 
and  rates  at  whii>  heat  diffuses  to  the  surface  and  carbon  vapo-s  diffuse 
into  the  main  gas  stream.  Chemical  reaction  is  determined  by  adiabatic 
wall  temperature,  heat  transfer  coefficients  and  free-stream  concentra¬ 
tions  of  corrosive  species.  To  a  lesser  extent,  the  thermal  conduction 
of  the  throat  backup  material  and  the  view  factor  to  the  nozzle  exit  ar^ 
also  important. 
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Representative  erosion  rates  for  ATJ  graphite  throat  inserts  by  hydrogen 
and  by  water  vapor  are  given  in  Pigs.  45  and  46  ,  respectively,  as  func¬ 
tions  of  adiabatic  vail  temperature  and  .free-stream  mole  fraction  of  the 
reactant  gas  species.  The  erosion  rates  -were  calculated  by  a  Rocketdyne 
one-disensicnal  computer  program,  THAB,  employing  the  erosion  mechanisms 
of  2B-ABLATE  and  the  recommended  erosion  properties  from  Tables  5  through 
?.  The  calculations  assiane  a  2-inch  throat  diameter,  an  expansion  ratio 
of  40,  and  a  phenolic-refrasil  overwrap.  The  data  apply  directly  for  a 
chamber  pressure  of  500  psia.  Estimations  for  other  chamber  pressures 
may  be  obtained  from  the  graphs  by  employing  an  effective  mole  fraction, 
Xeff»  defined  by 


eff 


-  (— )  X 

"  l500'  p 


where  X  is  the  free-stream  mole  fraction  at  the  chamber  pressure  of 

Pc 

interest. 


Representative  erosion  rates  for  silicon  carbide  throat  inserts  are  given 
in  Pig.  47*  The  method  of  calculation  and  assumed  operating  conditions 
are  those  described  in  the  previous  paragraph  for  ATJ  graphite.  Because 
the  erosion  of  silicon  is  assumed  to  be  strictly  by  vaporization,  the 
data  in  Pig.  47  are  independent  of  chamber  pressure  except  for  the 
effect  of  chamber  pressure  on  heat  transfer  coefficient  and  adiabatic 
wall  temperature. 


THROAT  RECESSION,  INCH/SECOND 


ADIABATIC  WAIL  TEMPERATURE ,  R 

Figure  46.  Recession  ol'  ATJ  Graphite  Throat  Inserts  in  a 
Hydrogen  Environment  at  500  psia 
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TWO-DIMENSIONAL  ABLATION 


Analysis  based  upon  one-dimensional  beat  transfer  permits  estimations  of 
the  required  thicknesses  of  ablative  vails,  durability  of  throat  r.nserts, 
and  insuletive  overwrap  for  an  ablative  engine  at  given  operating  condi¬ 
tions.  From  these  estimations  and  from  the  chamber  dimensions  specified 
by  chamber  pressure,  thrust,  and  expansion  ratio,  a  preliminary  ablative 
wall  configuration  can  be  drafted.  The  (Lar  rate  and  temperature  distri¬ 
bution  in  the  actual  chamber  walls  will  differ  from  one-dimensional  esti¬ 
mations  for  the  following  reasons: 

1.  The  boundary  conditions  along  the  combnstion  chamber  and  nozzle 
walls  vary  from  the  injector  to  the  nozzle  exit.  These  boundary 
conditions  include  adiabatic  wall  temperature,  heat  transfer  co¬ 
efficient,  and  radiation  view  factors. 

2.  The  boundaries  in  the  throat  region  are  not  normal  to  the  radial 
heat  flux  vector  of  a  one-dimensional  analysis. 

3.  The  differences  in  thermal  conductivity  between  typical  throat 
inserts  and  ablative  materials  amplify  differences  in  radial  tem¬ 
perature  profile  produced  by  axial  variation  in  boundary  condi¬ 
tions. 

4.  The  anisotropic  thermal  conductivity  of  reinforced  ablatives  will 
amplify  axial  conduction  produced  by  differences  in  radial  tem¬ 
perature  profile. 

To  determine  the  importance  of  the  two-dimensional  effects  in  ablative 
chamber  wells  under  steady  firing  conditions,  seven  computer  runs  were 
conducted  with  2D-ABLATE  for  representative  thrust  chamber  configurations 
and  operating  conditions.  The  conditions  investigated  are  listed  in 
Table  11. 

Two  propellant  combinations,  ClF^/NgH^  and  NgO^/N^-UDMHC 5-0-50)  are  con¬ 
sidered.  Carbon  cloth/phenolic  ablative  walls  in  the  combustion  chamber 
mi  nozzle  together  with  an  ATJ  graphite  throat  insert  are  used  with  the 
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propellant  combination  because  of  the  compatibility  of  carbon, 
vith  3?  and  also  because  of  the  high  melting  paint  of  the  reinforcement. 
With  the  Xo0 (50-50)  propellant  combination,  phenolic-refrasil 
is  used  for  the  combustion  chamber  and  nozzle  vails;  either  phenolic-re¬ 
frasil  or  a  silicon  carbide  insert  is  used  at  the  throat.  The  operating 
conditions  (Table  ll)  include  both  high  and  low  chamber  pressure  (i.e., 
high  and  low  heat  transfer  coefficient)  and  high  and  low  mixture  ratios 
(i.e.,  high  and  low  adiabatic  wall  temperature). 


Thrust  Chamber  Geometry 

The  ablative  motor  configurations  considered  in  the  two-dimensional  com¬ 
puter  runs  are  shown  in  Fig.  48  through  51*  In  all  configurations, 
internal  wall  materials  in  the  combustion  chamber  and  throat  regions  are 
surrounded  by  an  overwrap  of  phenolic-refrasil  whose  reinforcement  is 
oriented  for  low  radial  thermal  conductivity.  For  convenience  in  computer 
programming,  the  boundaries  and  wall  material  interfaces  cousist  of 
straight-line  segments.  Although  these  geometries  are  idealized  and  im¬ 
practical  from  mechanical  strength  and  molding  considerations,  the  general 
results  of  the  thermal  analysis  of  the  simpler  geometries  are  applicable 
to  more  complex  boundaries. 

The  mixture  ratios  given  iri  Table  11  are  assumed  to  be  the  mixture  ratios 
close  to  the  chamber  wall,  so  that  the  flame  temperature  for  the  subject 
propellant  vombination  at  the  given  mixture  ratio  is  close  to  the  adiabatic 
wall  temperature  once  combustion  is  complete.  The  heat  transfer  coeffi¬ 
cient  in  the  throat  region  and  in  the  nozzle  is  defined  from  graphical 
correlations  given  in  Appendix  D.  Radiation  to  the  nozzle  exit  is  based 
upon  a  graphical  correlation  of  view  factors  given  in  Appendix  D. 

On  the  basis  of  the  analysis  of  the  OFg/MMH  ablative  chamber  in  Task  II, 
the  following  assumptions  were  made  for  the  two-dimensional  computations: 

1.  The  full  length  of  the  combustion  chamber  is  required  for  combus¬ 
tion;  i.e.,  the  local  adiabatic  wall  temperature  rises  linearly 
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Figure  48  .  Thrust  Chamber  Configuration  for  Computer  Run 
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from  a  value  of  2460  H  (2000  F)  at  tbe  injector  face  to  the  final 
adiabatic  wall  temperature  d» fined  by  the  mixture  ratio  in  the 
convergent  section  of  the  nf  zzle. 

2.  The  heat  transfer  coefficient  varies  locally  in  the  combustion 
chamber,  so  that  when  combined  with  the  A  T  between  the  true 
local  adiabatic  wall  temperature  and  a  wall  at  500  F  (regenerative- 
cooling  wail  temperature),,  the  same  heat  flux  is  predicted  as  is 
obtained  by  a  calculation  using  the  simplified  Bartz  equation 
(Ref.  23)  for  heat  transfer  coefficient  and  a  Al  based  upon 
the  thermodynamic  gas  temperature. 

The  adiabatic  wall  temperatures  after  complete  combustion  and  the  heat 
transfer  coefficients  at  the  nozzle  throat  are  listed  for  the  seven 
computer  runs  in  Table  12. 

TABLE  12 

ADIAEATIC  WALL  TEMPERATURES  AND  THROAT  HEAT  TRANSFER 
COEFFICIENTS  FOR  TWO-DIMENSIONAL  COMPUTER  RUNS 


Computer 

Run 

TADW 

(Mixture  Ratio), 

R 

Throat  Heat 
Transfer  Coefficient, 
Btu/in.2-sec-R 

1 

6780 

0.0018 

2 

6780 

0.00029 

3 

5000 

0.00031 

4 

4600 

0.0056 

5 

4500 

0.00038 

6 

5400 

0,00030 

7 

6000 

0.0051 
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Results 


The  final  char  front,  gas-side  surface  temperature  profile,  and  outside 
skin  temperature  profile  after  a  300-second  bum  time  are  graphically 
illustrated  in  Fig.  52  for  tbe  chambers  lined  with  the  high-conductivity, 
phenolic  carbon/cloth  ablative.  Iu  all  cases,  the  char  front  has  penetrated 
completely  through  the  carbon  cloth  ablative  into  the  phenolic-refrasil 
overwrap.  Around  the  chamber  and  throat  regions,  where  the  overwrap  of 
the  low-conductivity  75-degree-oriented  phenolic-refrasil  is  1.0  inch  thick, 
the  outside  skin  temperature  is  near  or  below  100  F.  In  the  nozzle,  where 
the  overwrap  tapers  to  zero  thickness,  the  skin  temperature  climbs  to  tem¬ 
peratures  over  1000  F.  The  higher  skin  temperatures  in  the  nozzle  also  re¬ 
flect  the  change  in  orientation  between  the  overwrap  reinforcement  and  the 
nearest  hot-gas  surface  due  to  the  expansion  half-angle  of  the  nozzle. 

In  all  cases  shown  in  Fig.  52,  the  inner  surface  temperatures  rise  in  the 
combustion  chamber  because  of  the  increase  in  adiabatic  wail  temperature. 

In  the  nozzle  region,  where  the  adiabatic  wall  temperature  is  constant,  the 
inner  surface  temperatures  fall  because  of  the  decrease  in  heat  transfer 
coefficient.  The  difference  in  the  inner  surface  temperatures  shown  in 
Fig.  52  between  computer  Runs  1  and  2,  which  have  equivalent  adiabatic 
wall  temperature  distributions,  shows  the  effect  of  chamber  pressure  (heat 
transfer  coefficient)  upon  surface  temperature. 

Char  front  and  inside  and  outside  surface  temperatures  for  phenolic-refrasil 
lined  chambers  after  a  300-second  burn  time  with  NgO^/NgH^-UDMH  (50-50)  at 
a  chamber  pressure  of  1000  psia  and  mixture  ratios  of  1.0  and  2.0  are  graph¬ 
ically  illustrated  in  Fig.  53*  Because  of  the  lower  thermal  conductivity 
of  phenolic-refrasil,  the  char  depth  is  not  as  high  ae  with  the  phenolic 
carbon-cloth  systems.  Near  the  injector  end,  the  caar  front  does  not  pene¬ 
trate  through  the  inner  chamber  liner.  The  lower  thermal  conductivity 
also  reduces  the  outer  skin  temperature. 

For  high  heat  transfer  coefficients,  the  char  depth  can  be  calculated  by 
one-dimensional  relations.  In  Fig.  54  and  55,  the  char  penet,  etion 
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calculated  by  2D-ABLATE  in  computer  Run  1  fsr  phenol? c/carbon  cloth  at 
a  chamber  pressure  of  500  psia  is  compared  to  the  char  penetration  cal¬ 
culated  by  the  one-diaeosionui  graphical  relations  of  Fig,  25  through 
44  .  In  the  one-dimensional  calculation,  local  values  of  adiabatic  wall 
temperature  and  heat  transfer  coefficient  were  used  together  with  a  radial 
thermal  conductivity  calculated  from: 

K  =  Kn  cos"  9  *■  Kt  ain'^  6 
y  v  5 


-  *• 


i 

i 


! 
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and  the  recommended  values  of  major  and  minor  axis  thermal  conductivity 
from  Task  2,  Char  rate  in  the  overwrap  layer  was  calculated  by  combining 
the  thermal  resistance  of  the  totally  charred  inner  layer  with  the  gas- 
side  resistance  to  yield  an  overall  heat  transfer  coefficient  defined  by: 


1 

U 


A  x. 


(3-9) 


derived  in  Appendix  E,  As  shown  in  Fig.  54  and  55,  the  agreement  between 
the  one-dimensional  and  two-dimensional  calculations  of  char  depth  is 
rather  good  for  the  operating  conditions  investigated . * 

Comparison  of  the  results  of  one-  and  two-ditaensi onai  calculations  of  com¬ 
bustion  chamber  end  nozzle  char  rates  for  chamber  walls  of  phenol ic-refrasil 
at  a  chamber  pressure  of  1000  psia  is  given  in  Fig.  56  through  59  .  The 
agreement  is  not  os  good  as  with  the  phene lic/carbon  cloth  but  is  certainly 
satisfactory  for  preliminary  design  purposes,  considering  the  uncertainties 
in  the  specification  of  adiabatic  wall  temperature  and  heat  transfer  coef¬ 
ficient  in  real  rocket  engines. 

Operation  at  low  chamber  pressures  (low  heat  transfer  coefficients)  intro¬ 
duces  differences  between  one-  and  two-dimensional  calculations  of  char 
rote  (Fig.  60  and  6l).  In  this  case,  the  char  rates  predicted  by  2D- 
ABLATE  are  significantly  higher  through  phenolic/carbon  cloth  inner  cham¬ 
ber  liners  than  those  predicted  by  one-dimensional  calculations  at  a  cham¬ 
ber  pressure  of  50  psia.  The  difference  shown  in  Fig.  60  and  6l  results 
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1'igure  jb.  Comparison  of  Char  Advance  Through  Phenolic  Refraail 
Combustion  Chamber  WallB  by  Two-  and  One-Dimensional 
Computer  Programs  at  7-0  Inches  From  Injector 
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Figure  l'fj.  Comparison  of  Char  Advance  Through  Phenolic  Refrasil 
Combustion  Chamber  Wall.:  By  ‘it/o-  and  One-Dimensional 
Computer  Programs  at  8.5  Inches  Fr,.a  Injector 
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CONDITIONS: 
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Figure  jG.  Cczp? rison  uf  Char  Advance  in  Combustion  Chamber 
By  Two-  and  One-Dimensional  Computer  Programs  at 
8.5  Inches  From  In, lector 
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CONDITIONS: 
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Figure  59.  Comparison  of  Char  Advance  in  Combustion  Chamber 
By  Two-  and  One-Dimensional  Computer  Programs  at 
7-0  Inches  From  Injector 
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Figure  60  . 
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Comparison  of  Char  Advance  in  Combustion  Chamber 
By  Two-  and  One-Dimensional  Computer  Programs  at 
3.6  Inches  From  Injector  —  Run  2,  Low  Chamber  Pressure 
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Figure  Cl.  Comparison  of  Char  Advance  in  Combustion  Chamber 

By  Two-  and  One-Dimensional  Technique  at  3.6  Inches 
From  Injector  —  Run  3,  Low  Chamber  Pressure  — - 
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from  the  combined  effects  of  the  axial  gradient  in  adiabatic  wall  tempera¬ 
ture  through  the  combustion  chamber  and  the  high  internal  conductivity  in 
the  phenolic  carbon  cloth  wall  relative  to  the  surface  heat  transfer  coef¬ 
ficient  at  low  chamber  pressures.  The  net  effect  is  conduction  of  extra 
heat*  from  the  throat  region  (where  the  adiabatic  wall  temperature  and  beat 
transfer  coefficient  are  higher)  to  the  combustion  chamber  j'alls  and,  thereby, 
production  of  additional  charring,  above  that  produced  Lv  radial  heat  trans¬ 
fer.  The  anisotropic  conductivity  of  reinforced  ^blatrves,  particularly 
favor  this  type  of  conduction  when  the  fibers  are  oriented  r*  appreciable 
angles  from  the  radial  chamber  coordinate. 

At  low  chamber  pressures,  the  char  rate  through  soft  throat  sections  is 
calculated  to  be  higher  by  2D-ABIATE  than  is  predicted  from  one-dimensional 
calculations  (Fig.  62  and  63).  In  this  case,  the  higher  char  rates  are 
apparently  due  to  the  larger  surface  area  for  convective  heat  transfer  in 
the  throat  region  relative  to  the  area  defined  for  one-dimensional  radial 
heat  transfer,  i.e.,  the  throat  acts  like  a  fin.  At  high  chamber  pressures, 
where  the  surface  resistance  to  heat  transfer  is  negligible  lompared  to  the 
internal  resistance  in  the  walls,  the  fin  effect  of  the  throat  is  insig¬ 
nificant. 

Gas-side  surface  temperatures  at  relatively  high  and  relatively  low  chamber 
pressures  are  presented  in  Fig.  64  and  65  •  At  high  chamber  pressures, 
the  normalized  difference  between  adiabatic  wal]  temperature  and  sn'iace 
temperature,  T^^  -  ^gA^p^  -  1200, drops  below  0.10  in  less  than  100  seconds 
of  burn  time  except  in  the  expansion  region  of  the  nozzle.  This  means  that 
surface  temperatures  will  typically  approach  the  adiabatic  wa*I  temperature 
to  within  F  unless  lowered  by  surface  reactions.  Because  surface  re¬ 
actions  will  involve  some  degree  of  erosion,  a  conservative  preliminary 
estimation  of  the  onset  of  throat  erosion  should  be  made  on  the  basis  of 
erosion  relations  presented  in  Fig.  45  through  47  at  various  free-stream 
concentrations  of  corrosive  species  in  the  combustion  gas  streams. 


*In  any  conduction  process,  hi _vi  resistance  at  the  boundaries  relative 
to  the  internal  resistance  ^ars"  the  internal  temperature  profiles. 
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Figure  Gh  .  Approach  oi  Surface  Temperature  to  Adiabatic  Wall 
.  Temperature  for  Relatively  High  Chamber  Pressure 


Figure  .  Approach  oi  Surface  Temperature  to  Adiabatic  Wall 
Temperature  for  Relatively  Low  Chamber  Pressure 
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According  to  Fig.  65,  the  gas-side  surface  temperatures  will  be  appreciably 
belov  che  adiabatic  vail  temperature  in  all  portions  of  the  chamber  at  low 
chamber  pressures  (<  50  psia; .  A  preliminary  estimation  of  throat  ero- 
sion  based  upon  the  adiabatic  wall  temperature  will  therefore  be  con¬ 
servative  at  these  operating  conditions. 


Intermittent  Duty  Cycles 


To  examine  the  effects  of  a  multiple  start  duty  cycle,  Computer  Hun  8 
was  made  wild  the  thrust  chamber  eonfiguraf on  and  propellant  combina¬ 
tion  of  Computer  Run  1  (Table  11  and  Fig.  48)  but  with  a  burn  schedule 
as  shown  below: 


Time,  seconds 


Duty 


0  to  200 
200  to  2200 
2200  to  225G 
2250  to  2550 
2550  to  2580 


Fire  at  500-psia  chamber  pressure 
Soakback 

Fire  at  100-psia  chamber  pressure 
Soakback 

Fire  at  500-psia  chamber  pressure 


The  computer  run  was  terminated  during  the  final  burn  period  when  the 
time  limit  on  the  IBM-560  computer  was  exceeded. 


At  the  end  of  the  first  burn  period,  the  char  front  (1000  F  isotherm)  had 
penetrated  approximately  0.20  inch  ir.oo  the  phenolic  refrasil  overwrap 
around  the  combustion  chamber  and  approximately  0.30  inch  into  the  over- 
wrap  in  the  throat  region.  In  the  nozzle  region  charring  was  complete 
at  area  ratios  greater  than  3-0. 

At  the  end  of  the  2000-second  first  soakback  period,  the  char  front  ex¬ 
tended  approximately  0.50  inch  into  the  refrasil  overwrap  arouna  the 
combustion  chamber  with  essentially  the  same  penetration  in  the  throat 
region.  In  the  nozzle,  the  wall  was  completely  charred  through  for  area 
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ratios  between  2.8  and  8.  In  the  combustion  chamber,  charring  continued 
for  practically  the  entire  soakback  period.  The  entire  second  burn  period 
was  required  to  reheat  the  char  layer  in  the  combustion  chamber  and  throat 
regions  with  no  additional  charring  occurring  during  this  period.  At  the 
end  of  the  second  soakback  period  (300  second) ,  the  char  front  had  pene¬ 
trated,  approximately  0.58  inc'h  into  the  refrasil  overwrap  in  the  combustion 
chamber  and  0.65  inch  in  the  throat  region.  Charring  was  still  in  progress 
at  the  end  of  the  soakback  period. 

During  the  first  soakback  period,  the  additional  char  penetration  in  the 
combustion  chamber  was  approximately  45  percent  of  the  total  penetration 
through  both  the  inner  carbon  cloth  and  outer  refrasil  ablative  materials 
which  took  place  during  the  first  burn  period;  but  it  was  2~l/2  times  the 
penetration  through  the  refrasil  overwrap  alone.  During  the  first  500  . 
seconds  of  soakback,  the  char  ra*e  in  the  refrasil  was  almost  equal  to 
that  at  the  end  of  the  burn  period.  In  the  throat  region,  the  additional 
char  penetration  through  the  refrasil  during  the  first  soak  period  was 
approximately  70  percent  of  the  penetration  through  the  refrasil  during 
the  first  burn  period.  The  reason  for  the  high  additional  char  penetration 
during  soakback  for  the  conditions  of  Computer  Run  1  is  the  high  thermal 
conductivity  of  both  the  inner  carbon  cloth  chamber  liner  and  the  ATJ- 
graphite  throat  insert.  As  a  consequence  of  their  high  conductivities, 
these  portions  of  the  chamber  wall  attain  a  relatively  uniform  temperature 
close  to  the  adiabatic  wall  temperature  and  can  store  a  significant  amount 
of  heat.  Because  of  the  low  conductivity  of  phenolic  refrasil,  it  requires 
long  soakback  periods  for  this  heat  to  be  dissipated.  The  high  thermal 
conductivity  of  the  inner  wall  materials  together  with  the  low  conductivity 
of  the  phenolic  refrasil  also  produce  the  long  heat-up  period  before 
charring  resumes  upon  ablative-motor  restart. 

At  the  end  of  the  first  bum  period,  the  outer  skin  temperature  of  the 
phenolic  refrasil  had  risen  only  from  70  to  80  F  in  the  chamber  and 
throat  regions,  but  had  risen  to  2760  F  at  an  area  ratio  of  5.  At  the 
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end  of  the  first  soakback  period,  the  skin  temperature  in  the  throat  region 
had  risen  to  520  F,  but  the  nozzle  skin  temperature  had  dropped  to  1975  F. 
During  the  subsequent  burn  and  soakback  periods,  the  chamber  and  throat 
region  skin  temperatures  were  almost  constant,  but  the  nozzle  tempera¬ 
tures  fluctuated  with  heating  load. 

The  results  of  the  multiple  restart  computer  experiment  described  above 
are  specific  to  the  engine  configuration  and  mission  cycle  specified. 
However,  they  reveal  the  general  qualitative  trends  to  be  encountered 
with  this  type  of  operation. 


The  duty  cycle  of  the  Apollo  Command  Module  rocket  engine  which  provided 
the  data  for  evaluating  the  thermal  conductivity  of  phenolic-refrasil 
under  Task  2  of  this  study  furnishes  an  example  of  a  combined  pulse  mode 
and  multiple  restart  operation.  The  portion  of  the  duty  cycle  of  the 
Apollo  engine  examined  and  the  calculated  char  penetration  through  the 
combustion  chamber  walls  during  the  periods  of  pulse,  full  burn,  and 
soakback  operation  are  summarized  in  Table  13.  Charring  occurs  during 
both  burn  and  soakback  periods  with  the  actual  burn  time  for  the  entire 
duty  cycle  being  74.7  seconds.  The  adiabatic  wall  temperature  and  heat 
transfer  coefficient  in  the  Apollo  engine  are  4000  F  and  0.00048  Btu/in.^- 
sec-F  in  the  combustion  chamber  just  ahead  of  the  nozzle  throat,  while 
the  radial  char  thermal  conductivity  for  the  phenolic-refrasil  walls 
(45-degree  reinforcement  orientation)  is  approximately  1.2  x  10”'* 
Btu/in-sec-F.  From  one-dimensional  calculation  for  the  actual  burn 
time  (Figs.  25  and  31  )>  the  char  depth  is  approximately  0.23  inch, 
while  rigorous  calculation  by  means  of  2D-ABLATE  results  in  a  char 
depth  of  0.57  inches. 

The  general  conclusion  to  be  drawn  from  the  intermittent  duty  cycle  modes 
of  operation  examined  is  that  the  char  rates  are  signif icantly  different 
from  the  rates  under  steady  firing  conditions  and  require  calculation  by 
numerical  methods  which  consider  the  specific  geometrical  arrangement 
and  type  of  burn  cycle  of  the  system. 
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Start  of  Interval, 
Seconds 


End  of  Interval, 
Seconds 


Fraction  Burn  j  Char  Depth  at 
Daring  Interval  End  of  Interval 


.0 

180 

0.08 

0.19 

180 

350 

0.000? 

0.19 

350 

363.4 

]  .00 

0.22 

363-4 

768 

0.0725 

0.41 

768 

1052 

0.0 

0.44 

1052 

1132 

0.089 

0.44 

1132 

1324 

0.0 

0.44 

1324 

1440 

0.008 

0,44 

1440 

1516 

0.123 

0.48 

1516 

i570 

0.00 

0.53  . 

1570 

1760 

0.0148 

0.57 

S0XHAHT  OF  DESIffi  PHOCH3UHES 


To  effectively  employ  2D-ABL4TE  and  the  parametric  relations  in  this 
report,  the  following  general  procedures  are  recoanended: 

J.  Convert  engine  specifications  into  heat  transfer  boundary 

conditions.  The  engine  specifications  are  propellant  eombi- 
nation,  mixture  ratio,  chamber  pressure,  thrust,  expansion 
ratio,  duty  cycle,  and  proposed  injector.  Fron  these  speci¬ 
fications,  the  approximate  chamber  vail  contours  are  defined 
before  a  heat  transfer  design  begins.  The  necessary  heat 
transfer  boundary  conditions  are  adiabatic  wall  temperature 
profile,  heat  transfer  coefficient  profile,  and  the  internal 
radiation  view  factors. 

The  adiabatic  wall  temperature  is  determined  from  the  flame  ~  * 
temperature  of  the  propellant  combination  at  the  operating 
nixtore  ratio.  Flame  temperature  is  obtained  from  any  avail¬ 
able  propellant  performance  program  such  a3  the  Hocketdyne 
N-eleneni  program  (Bef.  24). 

Local  heat  transfer  coefficients  along  the  chamber  axis  can  be 
obtained  from  the  graphical  relations  of  Appendix  D  based  upon 
the  simplified  Bartz  equation  (Bef.  23)  ssing  combustion  g gs 
properties  from  propellant  performance  data.  Badiation  view 
factors  for  the  throat  and  nozzle  regions  are  obtainable  from 
graphical  correlations  such  as  those  given  in  Appendix  D. 

2.  Choose  ablative  wall  materials  on  the  basis  of  adiabatic  wall 
temperature  and  compatability  with  combustion  gases,  carbon- 
cloth  phenolics  for  high- temperature  gases  with  minimal  water 
content,  phenolic-refrasil  for  low  flame  temperature,  and 
water-containing  combustion  gases.  Reinforcement  orientation 
should  be  30  degrees  or  more  from  chamber  axis  to  prevent 
de lamination. 
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Pron  the  average  thermal  conductivity  of  the  ablati'e  vali  chosen, 
the  adiabatic  wall  iesperature,  and  gas-side  heat  transfer  coef¬ 
ficient,  estimate  the  expected  char  depth  fros  the  curves  of 
Pig.  25  through  44  .  Interpolate  between  corves  by  scans  of 
Eq.  5-1,  5-2,  3-5,  3-6,  3-7,  and  3-S.  Extrapolate  to  longer 
bam  times  by  assnsing  that  char  depth  varies  as  the  square 
root  of  tine  for  long  bum  time-  Hake  an  initial  allowance  of 
50  percent  sore  char  depth  for  intermittent  dety  cycle. 

If  char  depth  is  considered  excessive,  try  overwrap  of  phenolic 
refraHl  with  orientation  nearly  parallel  to  chamber  axis  for 
lew  char  rate.  Calculate  char  rate  in  the  o-yerwr^p  layer  from 
Bq.  3-9,  starting  the  calculation  at  the  tine  charring  of 
inner  vail  layer  is  complete- 

Por  chamber  preesores  above  50  psia  and  long  duty  cycles, 
consider  hard  throats:  graphite  for  i 1 nominated  oxidizers, 
silicon  carbide  for  combustion  gases  with  appreciable  oxygen 
or  water  vapor.  Estimate  throat  erosion  from  Fig.  45  through 
4?  •  If  erosion  is  excessive  by  order  of  magnitude,  chamber 
pressure  ecd/or  mixture  ratio  will  have  to  be  changed.  If 
erosion  is  marginal,  two-dimensional  calculation  may  change 
prediction. 

Estimate  thermal  penetration  depth  from  Eq.  3-3  and  Fig.  29. 

X?  limitations  on  the  outer  skin  temperature  exist,  sake  initial 
estimation  of  additional  overwrap  required  over  and  above  char 
depth  so  that  thermal  penetration  does  not  reach  outer  skin. 

Lay  out  exact  coatyurs  of  wall  segments  based  upon  stress  and 
fabrication  considerations. 

Test  configuration  with  2B-ABLATE.  If  erosion  is  excessive, 
chamber  pressure  and/or  mixture  ratio  must  be  changed.  Either 
hulk  mixture  ratio  or  mixture  ratio  near  the  wall  (by  injector 
modification)  may  be  changed.  Adjust  thicknesses  of  inner  and 
outer  ablative  layera  from  char  and  thermal  penetration  results. 
Consider  outer  insulation  when  necessary. 
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'  ANALYZE  WITH  2-D  ABLATE 


MODIFY  THERHAL  DESIGN  I 


MODIFIED  DESIGN 


ANALYZE  WITH  2-D  ABLATE 


SUCCESSIVE  ITERATIONS 


FINAL  DESIGN 


Figure  66.  Design  Procedure  for  Ablative  Rocket  Engines 


RECOMMENDATIONS  FOR  FUTURE  EFFORTS 


The  main  area  of  uncertainty  in  the  design  of  ablative ly  protected  liquid 
engines  is  in  the  prediction  of  surface  recession  rates.  Theoretical  models 
have  been  drawn  up  for  surface  degradation  based  upon  three  primary  mechanisms: 

1.  Melting  of  the  surface  leading  to  mechanical  removal  of  char 
through  shear  and/or  pressure  gradient 

2.  Chemical  reaction  with  the  combustion  gas 

3.  Sublimation 

The  present  2-D  program  includes  provisions  for  all  these  mechanisms.  The 

theoretical  models  on  which  the  calculations  are  based  have  been  shown  to  be 

in  reasonably  good  agreement  with  simple  experiments.  However,  because  all 

/  -A  E^kTs 

models  contain  Arrhenius-type  expontial  terms  ),  they  are  very 

sensitive  to  c*  efficiency  and  to  mixture  ratio  which  determine  gas 
temperature.  Moreover,  because  the  local  mixture  ratio  affects  the 
surface  reaction  (rather  than  the  overall  mixture  ratio),  the  reported 
bulk  mixture  ratios  given  in  the  usual  engine  test  reports  are  not 
necessarily  an  accurate  measure  of  the  ablative  surface  environment. 
Consequently,  most  of  the  data  from  small  engine  firings  are  relatively 
useless  for  checking  the  validity  of  computer  program  calculations. 

A  particularly  disturbing  observation  concerning  average  surface  recession 
rates  as  usually  reported  (total  change  in  cross-sectional  area  divided  by 
firing  time)  is  the  fact  that  "after"  photographs  often  show  scalloping 
attributed  to  oxidizer  streaking,  whereas  the  computer  program  assumes 
uniform  recession.  It  is  obvious  that  available  surface  recession  date, 
can  only  be  applied  by  a  person  thoroughly  familiar  with  the  injector, 
chamber,  and  instrumentation  employed  in  the  test. 

Several  combinations  of  analytical  and  experimental  investigations  could 
advance  the  state-of-the-art  knowledge  of  surface  recession.  SUch  inves¬ 
tigations  are  recommended  to  include  the  following: 
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1.  Critical  re-examination  of  available  ablative  motor  test  data;  e.g., 
establish  contact  with  the  original  investigators,  define  the 
additional  complications  in  the  reported  test  data  caused  by 
injector  effects  and  process  these  runs  through  the  2-P  program 

to  determine  whether  the  overall  growth  of  the  cross-sectional 
aretTcan  be  adjusted  to  a  uniform  recession  model  by  a  simple 
accomodation  coefficient.  The  restrictions  upon  the  use  of  such 
a  coefficient  would  be  established  possibly  in  terms  of  an  in¬ 
jector  spray  pattern  rating. 

2.  Other  analytical  attempts  to  describe  nonuniform  erosion  with 
the  help  of  the  2-D  program. 

3-  Generation  of  additional  data  in  critical  regions  using  hard¬ 
ware  specifically  designed  to  produce  a  uniform  mixture  ratio 
and  a  well-developed  flow  profile  throughout  the  ablative  portion 
of  the  nozzle.  Such  firings  would  check  the  validity  of  the 
uniform  recession  under  the  idealized  conditions  assumed. 

4.  Variation  of  the  wall  mixture  ratio  in  the  test  apparatus  of  3 
above  through  selective  hole  enlargement  in  the  injector.  The 
resultant  effects  on  throat  erosion  would  be  measured  in  a  few 
firings.  The  results  would  be  compared  to  the  predictions  of 
previous  hot-gas  mixing  experiments  obtained  at  Eocketdyne  and 
elsewhere  to  determine  the  possible  variation  in  mixture  ratio 
effects  in  an  ablative  chamber  design. 

5.  Extension  of  the  2D-ABIATE  to  cover  heat  transfer  and  erosion 
mechanisms  heretofore  not  included;  e.g.,  effect  of  solid  particles 
in  the  gas  stream,  radiation  from  gas  to  wall,  variation  of  the 
gas-side  heat  transfer  coefficient  with  time,  boundary  layer  cooling. 

6.  Extension  of  the  2D-ABIATE  to  include  calculation  of  view  factor; 
thermal  stress  analysis;  convective  heating,  charring  and  erosion 
at  any  boundary  with  a  more  flexible  mesh  procedure  to  handle 
alternative  configurations. 

It  should  be  noted  that  all  recommended  tasks  are  independent  of  v&ch  other; 
i.e.,  any  task  or  portion  of  a  task  may  be  executed  with  little  effect  on 
the  information  to  be  obtained  in  the  other  tasks. 
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NOMENCLATURE 


areas  of  thrust  chamber  inside  subsurfaces,  i=l,  . « . ,  K 

input  constant  characteristics  of  the  blocking  phenomenon 

specific  heat  of  vail  material 

free  stream  specific  heat  at  -constant  pressure 

fraction  of  mass  converted  to  gas  or  fraction  remaining 
as  solid  due  to  gas  generation  reaction  r  in  charring 
material 


=  view  factors  from  inside  subsurface  i  to  subsurface  k 

=  radial  position  of  receding  hot  gas  boundary 

=  axial  and  radial  components  of  generated  gas  mass  flax 
in  charring  material 

=  components  of  generated  gae  mass  flux  in  minor  and  major 
directions  of  conductivity  for  anisotropic  charring  material 

=  esbs  flux  of  gaseous  species  j  at  exposed  inside  surface 

=  enthalpy  of  gaseous  species  j  at  exposed  inside  surface 

=  enthalpy  of  gases  generated  in  charring  material 

=  heat  of  decomposition  mode  r  at  eroding  vail  surface 

=  heat  transfer  coefficient 

=  basic  cbnvective  heat  transfer  coefficient  prior  to 
modification  to  accourt  for  blocking 

=  thermal  conductivity  o:f  vail  material 

=  conductivity  in  the  normal  direction  at  the  surface  of 
an  anisotropic  material  (see  Eq.  23) 

=  conductivities  in  the  axial,  radial,  and  "mixed"  directions 
for  an  anisotropic  material  (see  Eq.  22) 

=  conductivities  in  the  minor  and  major  direction  for  an 
anisotropic  wall  material 
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n,  s 


=-•  outward  normal  and  counterclockwise  tangential  directions 


u  . 

J 

«r 


=  known  constant  for  gaseous  species  j  at  exposed  inside 
surface  (see  Eg.  9) 

=  reference  heat  of  reaction  for  gas-generation  reaction  r 
in  charring  material 

=  heat  flux 

-  rati°<  Vx'Vv 

=  temperature  of  mil  material 


T  (x)  =  adiabatic  wall  temperature 

T  (x,y,T)  =  maximum  value  of  T  achieved  by  time  T  at  point  (x,y)  in 

x  wall  material 


x,y 

X 

Ax,  Ay 
fif(x) 

6 

c 

p(T) 

v ,  $ 

or 

T 

AT 

e 

u 

Pr 


=  minimum  pyrolysis  temperature  of  charring  material 
=  normal  velocity  of  gas  or  liquid  at  exposed  inside  surface 

=  axial  and  radial  coordinates 
=  mole  fraction 

=  axial  and  radial  distance  increments 

=  radiation  view  factor  from  wall  surface  to  outside 
environment 

=  spatial  difference  operator 
=  emissivity  of  naterial  surface 
=  spatial  density 

=  major  and  minor  directions  of  conductivity  for  an  aniso¬ 
tropic  material 

=  Stef an-Bo ltzmann  constant 

=  time 

-  time  increment  used  for  time  step  calculations 

angular  displacement  of  £  and  fj  directions  from  x  and  y 
directions,  respectively,  for  an  anisotropic,  csterial 

=  overall  heat  transfer  coefficient 

=  Frandil  number 


vfiSf  •'  - - 
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Scb3crists 


A,  3,  C,  D  =  points  adjacent  to  interior  point  0  in  finite  differ¬ 
ence  analog  of  energy  eoasti^n 

cosy  =  doe  to  ccirrecticn 

efi  =  effective 


sit?  =  ace  to  environssiial  heating 


0 

r 


*  interior  point  (x^,  y^.} 

=  index  ragging  over  erosion  nodes  at  the  exposed  inside 
surface  or  over  gas  generation  reactions  in  a  charring 
cate rial 


rad  =  due  to  radiation 

rerad  =  dee  to  radiative  exchange  at  the  exposed  inside  surface 

v  =  virgin 

c  '  -  char 


Units 

When  not  otherwise  specified,  t._j  following  units  apply: 

h  =  Btu/ in,  -sec-P 

a 

k  =  ht u/in.-sec-P 
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For  the  checkout  case  described  earlier  is  this  report,  the  »eelin£iion 

of  the  inside  surface  frees  the  horizontal  did  sot  exceed  2D  degrees  so 

that  the  assumption  ei  so  occnltaiion  or  shading  {required  for  the  use 

of  Eq.  A-l)  is  reasonable.  Using  Eq.  A-l,  ve  can  calculate  the  desired 

array  F.  .,  i,j  »  1,...,  »,  as  foilovs: 
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3be  derivation  of  Eg.  .4-3  3rd  A-4  requires  use  of  as  interia  array  of 
^iev  factors  free  the  inside  surface  areas  to  the  base  areas  A^  , 

j 

i  «  1, — ,  X,  j  =  1, —  N  -  1,  given  by  ihs  following: 
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(A-5) 


fba  procedure:  vill  be  indicated  for  several  values  of  i  enH  k.  The  rest 
is  done  aiailarly.  To  get  F^  ye  start  by  uriting  tie  identities  resist¬ 
ing  surfaces  yithia  the  first  disk: 


F,  ,  +  F,  +  F, 

1,1  1,0^  1,0^ 


S’1  *  \’a2  =  1: 


P  .  +  F  =  1 . 

V1  V°t 

proa  Eq.  A-7  and  A-S  ana  the  identity 


A  F  =  A  F  , 
r  ?,8  s  s,r 


(A-6) 


(A-7) 


(a-s) 


(A-9) 


f  - 


ye  obtain  the  folloving  relationships: 


pi,a.  <A-10> 

Fl,<^  ’  \1,1  ~  Fo2.V  '  **  ’  ^°2  ”  *al\’e^  ^  h'  <A-U) 

Substitution  o f  F,  and  F.  in  Eq.  A-6  and  rearranging  vill  yield 
*■»  **1  *■  '“2 

Eq.  A-3  for  the  case  of  i  =  1. 
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Sinilarly,  to  get  F.  „  ve  start  by  writing  one  identity  based  upon  analysis 

—  5^ 

of  the  second  disk: 


'y*  *  V  - 1  <A-! 

as  veil  a?  two  others  obtained  by  treating  the  first  two  disks  as  a  single 
disk  (by  removing  side  ) : 

*  Fl,l  *  ?1,2  +  P1,CC  =  lf  ^A"1 

FCC,ftj  *  ?Gyl  +  F0^,2  *  1#  ^A_i 

By  combining  Eq.  A-6  and  A-ll  with  A-13,  and  Eq.  A-12  with  A-14,  we  obtsin 
the  following  expressions: 


1,2  +  Fl,OL  *  \>°2  ^  Ai’ 


s  +.  f  s  F 

°5»i  «r»V 

Fron  Eq.  A-16  and  A— 9  ve  get 


?l,OL  "  ^3^  ?0^,C^  ”  AOq?<^,OL^  ’  Al* 
and  finally,  fron  Eq.  A-15  and  A-17,  we  obtain 

Fl,2  “  [S  ^  ”  "S  *Fqv°2  ~  FaL»°3^  1  Al* 

Eq,  A-18  is  a  special  case  of  Eq.  A-4  for  i  =  1  and  k  =  1. 


In  a  giailar  fashion,  ve  can  derive  the  following: 


?1.3  +  Fl,°*  “  '  \  \><*J  /  *1 ' 


F  +  F  =  P 

V“l  \,1 


^  VV*4  “l  W 


"1* 


\A-XL) 


and,  finally,  the  folloving  special  case  of  Eq.  A-4  for  i  =  1  and  k  -  2: 

Fl»3  =  IS  ^?a2,CC5  ~  ?a2 ’°k  ~  A°1  ~  Al*  ^4-22^ 

The  remainder  of  the  ^  can  be  obtained  in  a  similar  fashion  (or  by 
induction). 


-4PPE®rX  B 


DEVELOPMENT  OF  THE  ANISOTROPIC  ENERGY  EQUATION 


To  derive  the  fore  of  the  anisotropic  energy  equation  as  giyen  by  Eq.  21 
and  the  expression  given  by  Eq.  25  for  the  normally  directed  conductivity 
at  a  boundary  of  an  anisotropic  material,  a  transformation  through  a  rota¬ 
tion  6  is  required  from  gradients  in  the  x  and  y  directions  to  gradients 
in  the  Z  and  7]  directions  and  vice  versa.  Geometrical  considerations 
yield  the  following  expressions: 


3  o  3  -  O  3 

■gr  =  cos  8  -g-  +  s  in 8  -gp 

(9-1) 

(B-2) 

fz  -  If  -  —  k 

(B-3) 

-  »in0  ||  +  cos 0  . 

(B-4) 

If  we  further  stipulate  that  the  two  coordinate  systems  possess  a  conaon 
origin,  then*  from  Eq.  B-l  through  B-4,  ve  can  also  write 

x  =  §  co bB  -  7)  sin9,  (B-5) 

j  =  Z  sin9  +  7]  cos6.  (B~6) 

To  derive  Eq.  21,  we  start  with  the  energy  equation  expressed  in  the  co¬ 
ordinates  Z  and  7j  and  than  iranefora  it  to  x  and  y  coordinates  using 
Eq.  B-l  through  B-6.  This  would  be  straightforward  were  it  not  for  the 
need  to  include  the  radial  effect  in  the  Z  £ad  7]  system,  which  makes  it 
more  difficult  to  get  started.  If  x  and  y  were  cartesian  coordinates 
rather  than  cylindrical ,  the  energy  equation  would  take  the  following 
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form  in  the  £  and  V  system  (vhere  for  convenience  of  notation  ve  omit  the 
subscript  "eff"  on  C): 


9T  3  ,  dT\  B  j  Bt,  dS  /  ot  „  3T  . 

=  (K£Sf'  -  dT  ^G4^|  *  Gi?S5?  3 


(B-7) 


A  strictly  geometric  argument  entailing  Eq.  B-6  above  would  show  that  the 
radial  effect  can  be  introduced  into  Eq,  B-7  by  adding  the  following  term 
to  the  right  hand  side: 


K 

_ L 

sin  8  + 


T^e  <8iE0lf  +  cose^), 


(B-3) 


where  the  expression  for  K  ,  the  radial  component  of  conductivity,  is  to 
be  determined  as  a  function  of  K^,  K^.,  and  0  during  the  transformation  of 
the  rest  of  the  equation  to  x  and  y  coordinates.  Thus,  ve  transform 
Eq-  B— 7 ,  augmented  by  the  radial  term  B~8,  to  x~y  coordinates  by  sub¬ 
stituting  from  Eq.  B-i  through  B-6,  as  follows: 

PC  ||  =  (cosQ  ^  +  sin0  |~)  [K|  (cos0  ||  -t-  sin0  |~)]  + 

(«sin0  |-  +  cos9  ^ {-sinS  ||  +  cosS  ||)]  +  ~  - 

H  [g^  (cos0  §  +  sin3  ||)  +  G?.  (-sin0  ||  +  coS0  ||)] 

=  [(*g  cos20  +  sin2©)  ||  ]  +  ^  [(K^  eiu2B  +  cos 20) 


5ta s 

y  By 


^  ”  H  [(G£  0080  "  G77  sin0)  §  +  (G|  sin8  +  G7?  c080)  lyj  + 


If  i>£  -  V  sin0  0080  ID +  b  [<K£  -  V  sin0  cos6 13 


|r(KxS)+i|(yKy|) 

V  J&.  +  9  3T 

zy  BxBy  dT  Bx  By’ 


dH  (r  3T  „  3Tx 
dT  (Gx  dT  +  Gy  By'  * 


(B-9) 


15? 


(B-10) 
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J55253ST 


1»  i  tt 


2^ 

4L 

«$• 


Sever*!  iiaite  difference  analogs  of  tie  arises  teggerstare  derivative, 
cPy/gx £y,  -sere  expressed  in  Eq.  25  rfcrcngh  35*  varying  ires  first  to 
second  order  aceeraxe  depending  npaa  the  geccsstry.  Is  xis  derivaiioxs 
biiosf,  ti»  ariied  difference  o~  at  a  regnlar  point  0  (see  Fig.  8  for 
labeling  of  the  points)  given  by  £q.  25  is  obtained  from  Taylor  series 
expressed  abeni  tbs  point  0  for  xhe  resrperainres  at  the  points  S,  *,  G, 
and  E-  Gn  the  other  hard,  the  rtaaicisg  analogs,  given  by  Bo.  27  ihroegh 
56,  vere  derived  fresa  Taylor  series  for  the  tessseratare  gradients  at  the 
points  A,  3,  C,  and  35.  The  latter  approach  proved  to  be  no  re  applicable 
to  the  nesh-bctsdaiy  configurations  explored  in  the  25-AELATF  progras, 
and  tics  Eq.  2o  Iras  cot  used  in  the  pregraa  except  throega  chance  vhen 
Fq.  2/  or  28  happened  to  reduce  to  Eq.  25  (vhieh  is  a  special  case  of 
Ba.  26). 


To  derive  Eq.  2a,  ve  express  the  texpera tares  at  the  points  E„  F,  S,  and 
3  (Fig-  ?a)  in  Tarlcr  series  abotrfc  a  regular  point  0.  explicitly  tncleding 
terns  np  through  fosrth  order.  Only  the  expression  for  will  be  exhibited , 
as  follows,  tha  resainder  being  easily  expressed  similarly  (the  subscripts 
on  T  refer  to  partial  derivatives  at  the  point  8): 


T  = 
E 


Tn  + 


4 Ax,?  ft 
A  x,0 


&7BTy,o 


Ax'- 


xx, 0 


Lx 


A 


^bVo 


Ay 


B 


yy,o 


Aza 

— _dL  t 
6  xrx,0 


Ax* 

-  T 

2  ^B-xxyjO 


syy,0 


AyB 


T. 


yyy»o 


Axa  Axaa  Axa  AyB  AyB 

24  *xxx3:,0+  5  yB^xxsy,0  +  2  2  ^xxyy,0  +  ^  XA  *1)  Txyyy,0  + 


24 


T  a 

yyyy.o 


(c-i) 
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CtaBOsnitig  1^,  Ip  za5  1^  ve  obtain 

3^  -  %  *  ?G  -  ^  =  (Ax^  -  Axc)  (Ay3  -r  Ay^)  [t_  0  - 


*XA~**£ 


a4  -  4ia  Aic  *  Alc 


‘ayy.O 


?  _  * 
xszy.O 


?!_,  n  + 

^xxy.O 


(C-2) 


*Z£-*XC  ~'V3  -  A'VI 
2  2 


*4  '  AyD  *  *4 


xxyy.O  ~ 


,0  +*'€ J 


\  ~  \  +  fG 


'ZA  +  “XC'  '^'B  +  ^yD' 


=  TW  n  +  0  (*x4  -  Axr>  +  0  (Ay*  -  AyJ  + 


0  (Ax|  -  Axa  Axc  +  Ax2  )  +  0  [(Axa  -  Axc)  (AyB  -  Ay^]  * 


0  (Ay^  -  AyB  Ay^AyJ  ) 


(C-3) 


O 

The  left  side  of  Eq.  C-3  is  the  mixed  difference  analog,  6^  Tq,  as  ex¬ 
pressed  by  Eq.  26,  and  the  right  side  indicates  how  closely  the  difference 
analog  approximates  the  mixed  derivative,  T  _.  If  Ax  »  Axn,  Ax,  « 
Axp,  Ayg  »  Ayp,  or  Ay^  «  Ay^,  Eq.  C-3  is  essentially  first  order 
accurate.  If,  on  the  other  hand,  Ax  =  Ax„  =  Ax  and  Ay  =  Ay  =  Ay, 

A  t  d  D  q 

Eq.  C-3  reduces  to  Eq.  25  and  the  truncation  error  reduces  to  0  (A  XJ  + 

0  (Ay2). 


The  derivation  of  Eq.  27  and  28  is  considerably  simpler,  requiring  fewer 
terms  of  the  Taylor  series  and  yet  yields,  in  general,  a  more  accurate 
formula.  Here  we  express  temperature  gradients  instead  of  temperatures 
at  the  points  A,  B,  C,  and  D  in  Taylor  series  about  the  point  0,  as 
follows: 


Ty,A  “  7y,0  +  AzA  Txy ,0  +  2  Txxy 


,0  +  T"  Txxxy,0  +-‘*  fc"4) 
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A  —  4* 

S3 

*4 

CXIX,0  ~ 

2 

Jl  * 

xxy,0 

-  ~r 

a-vbV,o  - 

a4 

2 

*? 

*3yy,0 

Av5 
‘r  “ IT 

AV^,0  + 

2 

T 

xyy,Q 

'  6 

From  Eq.  C— 4  and  C-5-.  vs  obtain 


■4  Va-(6=€  -  \,a  -  l4  tv 


5a  Cxc  aV 


£x.  &xp 

- 7 - T _ _  n  T  n  +  0  (£  x ,  L  x_) 

o  xxx y,0  xy,0  'AC' 


and,  similarly,  from  C-6  and  C-7,  we  have 


0  (&yb  ^yD)  (c-9 

Eq.  C-8  is  applicable  to  the  2D-ABLATE  program,  as  discussed  between 
Eq.  28  and  29  in  the  text,  whenever  both  the  points  A  and  C  are  regular 
and,  similarly,  Eq.  C-9  is  applicable  when  points  B  and  D  are  regular. 

To  obtain  Eq.  27  and  28  from  Eq.  C-8  and  C-9>  respectively,  we  replace 

Ty,0»  Ty,A>  Ty,C’  Tx,0’  Tx,B’  and  Tx,D  ia  ^  C_S  and  °"9  by  8econd  order 
accurate  centered  finite  difference  analogs.  In  doing  this,  however,  ex¬ 
amination  of  the  error  term  involved  in  each  would  seem  to  indicate  that 
Eq.  28,  although  second  order  accurate  in  the  x-direction,  is  in  all  cases 
only  first  order  accurate  in  the  y-direction,  and  that  Eq.  29  similarly,  is 
only  first  order  accurate  in  the  x-direction.  That  this  cannot  be  true 
and,  indeed,  leads  to  contradiction  is  demonstrated  by  analyzing  the  case 
in  which  all  the  points  diagonally  adjacent  to  the  point  0  are  regular 
(Pig.  9a).  In  that  case,  Eq,  28  and  29  reduce  to  the  same  difference 
equation.  But  then  the  asseition  above  concerning  the  error  would  lead 


to  the  conclusion  tint  the  resulting  difference  equation  both  is  sad  is 
not  second  order  accurate,  vhitb  is  absurd.  As  a  natter  of  fact,  by 
deriving  Eq.  28  or  29  for  the  case  of  Fig.  9a  frees  Taylor  series  for  the 
tenpsratures  at  the  points  A,  3,  C,  D,  Z.  F,  G,  and  H,  the  error  can  be 
shave  to  he  0  (£xa£xc)  *  0  This  is  an  explicit  demon¬ 

stration  of  second  order  accuracy  and  proves  that  Eq.  28  and  29  are  sore 
accurate  than  Eq.  26.  Vfcen  the  points  diagonally  adjacent  to  the  point 
0  are  not  regular  (in  -which  case,  Eq.  26  vonld  not  even  be  appJ '  „ab2e), 
it  can  be  shown  that  Eq.  28  is  between  first  and  second  order  accurate 
in  the  y-direchion  (in  a  sense  sirsilar  to  Bq.  C-5)  and  Bq.  29  is  the  sane 
in  the  x-direction.  For  the  point  configuration  depicted  in  Fig.  C-l, 
for  example,  the  error  in  Eq.  28  is  0  (Ait  A  x^,)  ~  0  (Ay^  -  A  y^,  )  + 

0  (Ayg  iA  y_j) ,  which  yields  essentially  first  order  accuracy  in  the  y- 
direction  for  Ay  «  A  y  and  second  order  accuracy  for  Ay  ~  A  y_. 

As  stated  above,  either  Eq.  28  or  29  does  not  apply  when  one  or  sore  of 
the  points  A,  B,  C,  and  D  is  irregular  (Fig.  £r).  If,  for  example, 
point  C  is  irregular  we  would  obtain  the  following  expression  (instead 
of  Eq.  C-8  above)  directly  from  Eq.  C-4; 


+  0  (Axa) 


(c-io) 


As  with  Eq.  C-8,  when  we  substitute  the  second  order  accurate  difference 
analog  and  5  Tq  for  the  first  derivatives  in  Eq.  C-10,  the  resulting 
difference  analog  (given  as  part  of  Eq.  29)  is  between  first  and  second 
order  accurate  in  the  y-direction.  The  rest  of  the  expressions  given  by 
Eq.  29  and  30  are  similarly  obtained. 
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I 

% 


*  Ay 


*-  — - 9 


B 


o  Ax, 


Ay* 


S2SSS  TgT??} 


S?  ^77  ggfc?  -1 


*53  SOIiZHS  713?  £5£PiSS  10 
gjgSSES  iSTBSSSSiS? 


JLt  trie  threat  el  a  roc&st  gross!*,  the  sinpiiJiea  camelaticc  of  Baris 
(Stf.  25)  ear  be  --cites  as  follsvs: 


(B-l 


v2sre  fc 


is 


reel 


2 

fer  eoefriei eat  (3sa/in.  -sec- ?)  at  ire  iarcat 


T?**  throat  feaeig,  5*,  *rd  -£&»  esasber  pressure,  ?  ,  are  defined  hj  She 

c 

*g^rjrnp  operating  cgaditicss.  Eae  specific  beat  of  ccffibssticEi  gas, 

C  ,  and  tbs  characteristic  exhaust  velocity  of  the  gas,  c*,  are  given 
p 

by  propellant  pcrforrsrce  data  st  ti2  operating  csndiiions  (Sef.  24). 


Propellant  pexfossEBce  data  also  provide  tie  cccbastior*  gas  ierperstare, 

^  y  sraragg  anl ftenlar  yoj g^i .  fmw 

vMch  tie  gas  viscosity,  jx,  tie  Bcsaatl  rather,  Rr,  and  tie  reference 
property  gsx^Beter,  <7,  are  defined. 


For  rest  applications,  tie  dimensional  groups  cf  £□.  B-l  can  be  obtained 

vith  sufficient  accuracy  freae  ample  graphical  ccrrelations.  Figure  B-l 
0  2 

shn-tfs  (l/B*)  *  as  a  function  of  P  and  the  rocket  engine  throat,  F. 
Figure  B-2  shews  (P^/c*)  *  as  a  function  of  (P ^/c*).  In  Pig.  B-5,  the 

viscosity  factor,  *  ,  is  plotted  rs  the  product,  H?  ,  while  Fig.  B-4 

0  £  § 
shovs  the  relation  between  (a/Pr  *  )  and  y. 


For  ablative  rotor  calculations, 
« 

related  to  h  by 


h^  at  the  axial  location  x  can  be 

s 


^  4  -r4J-* 


Figure  D-5.  Variation  of  View  Factor  Between  Nozzle  Surface  Element 
and  Nozzle  Exit  as  a  Function  of  the  Ratio  Between  Nozzle 
Exit  Radius  and  Local  Radius 
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vies  p„F_,  A3  Lj  tii®  beat  aesorbssS  pss  aait  to  less  of  tt~>  c-r^sis  charred 

ard  T  is  tfee  effective  cisr  fresi.  Viik  this  ac-^ei,  U  is  gives  sy 
c  ^  * 


1  J:  +  Ie 
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Equation  E-2  can  be  integrated  to  give 


is  doainaiiiog  so  that  the 


Fs?  long  firing  tine 3,  the  tern  ■  ■  „  ~A  g 

c*v  R  & 


••jne-direasienal  interpolation  equations  (Eq.  3-2,  5-6,  and  3-7)  are 
obtained.  The  best  agreement  betveen  the  one-disensional  calculations 


ac3  Eq.  K-4  occurs  vhen  5^  is  defined  as  120(1  F. 


For  charring  through  railtiple  layers,  Bq»  F-4  becomes 


2  ATeff  *2 


VA  *  H 


C 


where 


in  Eq,  £~5«  the  tine  begins  when  the  inner  layer  of  thickness  Y,  ins; 
cftffipletely  charred  through. 
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The  two-dimensional  ablative  heat  transfer  computer  program  generated  Tinder 
contract  AF04(6ll)-9714  was  refined  and  extended  to  handle  anisotropic  mate¬ 
rials,  more  than  one  charring  material  and  reradiation  at  the  heated  surface. 
Results  of  the  two-dimensional  program  were  compared  to  experimental  data  to 
determine  effective  values  of  material  properties  used  in  the  analysis  to 
simulate  gas  generation  and  cracking  reactions.  The  resulting  effective 
properties  for  two  ablative  materials  were  used  in  a  parametric  study  gener¬ 
ating  basic  information  for  the  design  of  ablative  systems  in  liquid  rocket 
thrust  chambers.  Graphs  and  charts  showing  the  variation  of  thermal 
penetration,  char  depth  and  surface  erosion  are  included. 
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